Fatsllfra 13 GGG 334

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 93

List des matières

Chapitre 1 : Les voiles en béton armé sous chargement sismique


Comportement et principe de dimensionnement.

1.1. Introduction .................................................................................................................. 19

1.2. Pourquoi utiliser des murs en cisaillement ........................................................................ 19

1.3. Les caractéristiques essentielles influençant le comportement des murs voiles BA : ...... 20

1.4. Mode de fonctionnement des murs voiles ......................................................................... 20

1.4.1. Mode de fonctionnements des murs voiles élancés ................................................... 20

1.4.1.1. Endommagements en flexion .............................................................................. 20

1.4.1.2. Endommagement en flexion/cisaillement ........................................................... 20

1.4.1.3. Endommagement en cisaillement........................................................................ 20

.................................................................................................................................................. 21

1.4.2. Mode de fonctionnement des murs voiles courts ....................................................... 21

1.5. Le principe de dimensionnement dans les codes de calcul : ............................................. 21

1.6. Principe de dimensionnement des murs voile ................................................................... 23

1.6.1. Dimensionnement en capacité .................................................................................... 23

1.6.2. Dimensionnement par création des rotules plastique ................................................. 23

1.7. Les niveaux de modélisation ............................................................................................. 24

1.8. Les échelles de modélisation ......................................................................................... 24

1.8.1. Approche globale ................................................................................................... 24

1.8.2. Approche macroscopique ....................................................................................... 25

1.8.3. Approche semi-globale .......................................................................................... 25

1.8.3.1. Approche multifibre ........................................................................................ 26

1.8.3.2. Approche multicouche .................................................................................... 26

1.8.3.2. Approche multicouche .................................................................................... 27

2.1. Introduction ....................................................................................................................... 30

2.2. Le processus de fissuration dans les éléments en Béton armé .......................................... 30

8
2.2.1. Formation des fissures ................................................................................................ 30

2.2.2. Ouverture des fissures ................................................................................................ 30

2.2.2.1. Cas de la fissuration accidentelle ........................................................................ 30

2.2.2.2. Cas de la fissuration systématique : ........................................................................ 31

2.3. Contexte règlementaire et méthodes associées ................................................................. 33

2.3.1. Aperçu rapide sur les approches codifiées de limitation du dommage ...................... 33

2.3.1.1. Les codes basés sur des approches en force ........................................................ 33

2.3.1.2. Les codes basés sur des approches en déplacement ............................................ 34

2.3.1.3. FEMA 356 ........................................................................................................... 34

2.3.1.4 AIJ ........................................................................................................................ 34

2.3.1.5 ASCE/SEI............................................................................................................. 35

2.3.1.6 ATC-40................................................................................................................. 35

2.3.1.7 Comparaison des différents codes et commentaires ............................................ 36

2.3.2. Codification des ouvertures des fissures .................................................................... 37

2.3.2.1. Eurocode 2........................................................................................................... 37

2.3.2.2. Le code ACI (American Concrete Institute) ....................................................... 39

2.3.2.3. CEB FIP 1990 [CEB93] : .................................................................................... 40

2.3.2.4. Norme SIA 262 ................................................................................................... 41

2.3.2.5. Fib model concrete (MC2010) ............................................................................ 42

2.3.3. Les approches expérimentales pour l’estimation des ouvertures des fissures ........... 43

2.3.3.1. Approche Chi et Kirstein..................................................................................... 43

2.3.3.2. Approche Etienne Gallitre ................................................................................... 43

2.3.4. La codification des efforts de cisaillement ................................................................. 44

2.3.4.1. Code BAEL ......................................................................................................... 44

2.3.4.2. PS 92 ................................................................................................................... 45

a. Mécanisme de flexion............................................................................................... 46

b. Mécanisme d'effort tranchant ................................................................................... 47

c. Résistance limite associée ........................................................................................ 47

2.3.4.3. E.C.2 .................................................................................................................... 47

9
a. Résistance au cisaillement sans ferraillage réparti qui correspond à la contribution
du béton : ...................................................................................................................... 47

b- La résistance à l’effort tranchant avec des armatures transversales réparties (non


inclinées) est la suivante :............................................................................................. 48

2.3.4.4. E.C.8 .................................................................................................................... 48

2.3.4.5. ACI 349-01 .......................................................................................................... 49

2.4.3.6. JEAG40611987 ................................................................................................... 49

2.4. Conclusion......................................................................................................................... 50

3.1. Introduction ....................................................................................................................... 52

3.2. Présentation du programme SAFE .................................................................................... 52

3.2.1. Caractéristiques de la maquette .................................................................................. 52

3.2.2. Condition aux limites ................................................................................................. 54

3.2.3. Essai de pseudo dynamique........................................................................................ 55

3.4. Les modèles utilisés pour la simulation ............................................................................ 57

3.4.1. Le modèle de Laborderie............................................................................................ 57

3.4.2. Le modèle multicouche .............................................................................................. 58

3.4.3. Le modèle fichant (micro-iso) ................................................................................... 59

3.4.4. Loi de comportement de l’acier ................................................................................. 59

3.3.5 Les paramètres utilisés pour la loi élasto-plastique sont : ........................................... 59

3.5. Histoire de chargement...................................................................................................... 60

3.6. Résultats et discussions ..................................................................................................... 61

3.6.1. Résultat push-over (Statique monotone) .................................................................... 61

3.6.1.1. Voile T10............................................................................................................. 61

3.6.1.2. Voile T12............................................................................................................. 62

3.6.2. Résultat cyclique ........................................................................................................ 62

3.6.2.1. Approches multi couches .................................................................................... 62

3.6.2.2 Approche multi fibres : ........................................................................................ 63

3.6.3 Comparaison approche macroscopique Vs. Expérimentation..................................... 65

3.6.3.1 Courbe effort-déplacement (statique-monotone) ..................................................... 65

10
3.6.3.2. Courbe effort déplacement (cyclique (T10)) ....................................................... 65

3.5.3.3 Courbe effort déplacement (cyclique (T12)) :...................................................... 68

3.6.3.4. Discussions des résultats globaux : ..................................................................... 71

3.7. Résultats en termes comportement local (endommagement-fissuration) ......................... 72

3.7.1. Cas d’un chargement push-over ................................................................................. 72

3.7.2. Cas cyclique ............................................................................................................... 74

3.7.3. Discussion des résultats locaux (cartes d’endommagement) ..................................... 79

3.8. Conclusion......................................................................................................................... 80

4.1. Introduction ....................................................................................................................... 82

4.2. L’évolution des ouvertures de fissures OUVFISS ........................................................ 82

4.3. Les résultats du comportement local (Ouverture des fissures) ......................................... 83

4.3.1. Cas statique ................................................................................................................ 83

4.3.2. Cas cyclique ............................................................................................................... 87

4.4. Comparaison des ouvertures des fissures avec la réglementation..................................... 91

4.5. Comparaison les ouvertures des fissures avec expérimental ............................................ 93

4.5.1. Discussion .................................................................................................................. 93

4.6. Comparaison effort-tranchant avec la réglementation internationale............................ 94

4.6.1. Cas statique ................................................................................................................ 95

4.6.2. Cas cyclique ............................................................................................................... 95

4.6.3. Commentaire .............................................................................................................. 96

4.7. Conclusion......................................................................................................................... 96

Conclusion et perspective......................................................................................................... 97

11
List des illustrations

FIGURE1. 1:DOMMAGES SUBIS PAR LES VOILES (A) RUPTURE DE SECTION A LA BASE DU VOILE (SEISME DE
BOUMERDES DU 21/05/2003) ; (C) ENDOMMAGEMENT AUX NIVEAUX DE NŒUD POTEAUX FORTE POUTRE
FAIBLE. (B-D) : RUPTURE EN CISAILLEMENT DES VOILES (SEISME DE L'AQUILA 06/04/2009. ............... 19
FIGURE1. 2::MODES DE RUPTURE DE VOILES ELANCES (DAVIDOVICI ET AL, 1985)...................................... 21
FIGURE1. 3: MODES DE RUPTURE DE VOILES COURTS (DAVIDOVICI ET AL. 1985). ...................................... 22
FIGURE1. 4:: COMPORTEMENT DUCTILE D’UN VOILE ELANCE. RUPTURE DUE A LA FLEXION. .......................... 22
FIGURE1. 5: COMPORTEMENT NON DUCTILE D’UN VOILE. RUPTURE DUE A L’EFFORT TRANCHANT.................. 23
FIGURE1. 6::PRINCIPE DE DIMENSIONNEMENT SELON LES REGLEMENTS (EC8 VS PS92)[XUAN HUY NGUYEN
,2006]. .............................................................................................................................. 23
FIGURE1. 7: MODELISATION SMART2008 : TROIS TYPES D'APPROCHES (A) APPROCHE GLOBALE ; (B-C)
APPROCHE LOCALE ; (D) APPROCHE SEMI-LOCAL (LERMITTE 2010). ................................................ 24
FIGURE1. 8:PRINCIPE DE FONCTIONNEMENT DE L’ELEMENT POUTRE MULTIFIBRE [CEDRIC G.2006]. ............. 26
FIGURE1. 9::PRINCIPE DE FONCTIONNEMENT DE L’ELEMENT POUTRE MULTICOUCHE. [MOHAMMED IBRAHIM
2003]. ............................................................................................................................... 27

FIGURE2. 1:PRESENTATION DE LA DEFORMATION D’ACIER (HAUT) ET DE BETON (BAS) EN FONCTION DE LA


LONGUEUR DE LA BARRE DANS LE CAS DE FISSURATION ACCIDENTELLE. [ETIENNE GALLITRE, 2008]. ....... 31
FIGURE2. 2: REPRESENTATION DES DEFORMATIONS DE L’ACIER ET DU BETON AVEC L’AXE DE LA BARRE D’ACIER EN
ABSCISSE DANS LE CAS DE LA FISSURATION SYSTEMATIQUE [ETIENNE GALLITRE, 2008]. ....................... 32
FIGURE2. 3: CONTRAINTE D’ACIER A LA FONCTION DE LEUR DEFORMATION MOYEN [ETIENNE GALLITRE, 2008].
.......................................................................................................................................... 32
FIGURE2. 4: ELEMENT SOLLICITE EN TRACTION. ................................................................................... 39

FIGURE3. 1:PLAN DE FERRAILLAGE (T10). .......................................................................................... 53


FIGURE3. 2:PLAN DE FERRAILLAGE (T10). .......................................................................................... 54
FIGURE3. 3: PLAN DE FERRAILLAGE (T12). .......................................................................................... 54
FIGURE3. 4: ILLUSTRATION DE PRINCIPE DU DISPOSITIF DE CHARGEMENT. ................................................. 55
FIGURE3. 5: DIFFERENTES APPROCHES UTILISEES DANS LA MODELISATION (A)MULTIFIBRE (B) MULTICOUCHE (C)
APPROCHE MACROSCOPIQUE.................................................................................................... 57
- FIGURE3. 6:REPONSE UNI AXIALE DU MODELE [LA BORDERIE, 1991]. [PANAGIOTIS.K, 2006]. ............. 57
FIGURE3. 7:LOI DE COMPORTEMENT DE L’ACIER ELASTO-PLASTIQUE. [D. COMBESCURE, 2001]. .............. 59
FIGURE3. 8:: DEPLACEMENT –TEMPS (T10)........................................................................................ 60
FIGURE3. 9: : DEPLACEMENT –TEMPS (T12). ...................................................................................... 60
FIGURE3. 12:LA REPONSE GLOBAL (EXPERIMENTAL VS NUMERIQUE (CYCLIQUE)) POUR VOILE (A) T10 ET (B) T12-
MULTI COUCHES. ................................................................................................................... 63
FIGURE3. 13:LA REPONSE GLOBAL (EXPERIMENTAL AVEC NUMERIQUE (CYCLIQUE)) POUR VOILE T10 (A) ET T12
(B)-MULTI FIBRES................................................................................................................... 64
FIGURE3. 14:COURBE DE CAPACITE MACRO AVEC SEMI-GLOBAL (MULTICOUCHE). ...................................... 65
FIGURE3. 15:COMPARAISON DIFFERENT CYCLES NUMERIQUES AVEC EXPERIMENTALES SAFE (T10) PAR
APPROCHE MACROSCOPIQUE.................................................................................................... 68
FIGURE3. 16: COMPARAISON DIFFERENT CYCLES NUMERIQUES AVEC EXPERIMENTALES SAFE (T12) PAR
APPROCHE MACROSCOPIQUE.................................................................................................... 71
FIGURE3. 17:EVOLUTION ENDOMMAGEMENT PAR TRACTION (CISAILLEMENT) VOILE T10. ........................... 73
FIGURE3. 18: EVOLUTION ENDOMMAGEMENT PAR TRACTION VOILE T12.................................................. 74
FIGURE3. 19:EVOLUTION D’ENDOMMAGEMENT POUR DIFFERENT CYCLE T10 ET COMPARAISON AVEC
L’EXPERIMENTAL.................................................................................................................... 77

12
FIGURE3. 20: EVOLUTION D’ENDOMMAGEMENT POUR DIFFERENT CYCLE T12 ET COMPARAISON AVEC
L’EXPERIMENTAL.................................................................................................................... 79

FIGURE4. 1:ILLUSTRATION SUR EVOLUTION DES OUVERTURES DES FISSURES POUR VOILE T10. ...................... 85
FIGURE4. 2:ILLUSTRATION SUR EVOLUTION DES OUVERTURES DES FISSURES POUR VOILE T12. ...................... 87
FIGURE4. 3:LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 0.78MM. ......................... 87
FIGURE4. 4:LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 5.43MM. ......................... 88
FIGURE4. 5: LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 11.57MM. ...................... 89
FIGURE4. 6:6 LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 13.50MM. .................... 89
FIGURE4. 7:LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 0.25 MM. ........................ 89
FIGURE4. 8:LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 0.85MM. ......................... 90
FIGURE4. 9: LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 2.60MM. ........................ 90
FIGURE4. 10:LA CARTE DE L’OUVERTURE DE FISSURES POUR UN DEPLACEMENT DE 7.8MM. ......................... 91

13
Liste des tableaux
TABLEAU2. 1: TABLEAU INDIQUANT LES LIMITES DES DOMMAGES SELON LE CODE DE CALCUL FEMA 365....... 34
TABLEAU2. 2: INDIQUE LES LIMITES DES DOMMAGES SELON LA REGLEMENTATION JAPONAISE. ...................... 34
TABLEAU2. 3: TABLEAU DONNE LES LIMITES DRIFT ENDOMMAGEMENT..................................................... 35
TABLEAU2. 4: DRIFT LIMITE DONNE PAR ATC -40................................................................................ 36
TABLEAU2. 5: SYNTHESE DES DRIFTS LIMITES DONNES DANS LES CODES. ................................................... 36
TABLEAU2. 6:DIAMETRE MAXIMAL Φ_S DES BARRES POUR LA MAITRISE DE LA FISSURATION [S.MULTON-
2012]. ............................................................................................................................... 38
TABLEAU2. 7: DIAMETRE MAXIMAL Φ_S DES BARRES POUR LA MAITRISE DE LA FISSURATION [S.MULTON-
2012]. ............................................................................................................................... 38

TABLEAU3. 1: DESCRIPTION DE LA MAQUETTE SAFE. ........................................................................... 53


TABLEAU3. 2:CARACTERISTIQUE GEOMETRIQUE ET MECANIQUE DES MAQUETTES SAFE. ............................. 56
TABLEAU3. 3: TABLEAU DES PARAMETRES DU MODELE UNILATERAL UTILISES DANS LA SIMULATION. [D.
COMBESCURE, 2001]. ....................................................................................................... 58
TABLEAU3. 4:LES PARAMETRES DE MODELE MULTICOUCHE. ................................................................... 58
TABLEAU3. 5: LES PARAMETRES DE FICHANT UTILISE DANS LA MODELISATION ........................................... 59
TABLEAU3. 6:TABLEAU DES PARAMETRES DU MODELE ELASTO-PLASTIQUE UTILISES DANS LA SIMULATION. ...... 59

TABLEAU4. 1: IDENTIFICATION PARAMETRE DE MODELE (OUVFISS) UTILISEE DANS LA SIMULATION. ............ 83


TABLEAU4. 2: COMPARAISON ENTRE DIFFERENT REGLEMENTATIONS AVEC LES APPROCHES DE MODELISATION. 93
TABLEAU4. 3: COMPARAISON LES OUVERTURES DES FISSURES (MC2010, MACRO, MULTICOUCHE ET
EXPERIMENTALE) POUR T10 ET T12. ......................................................................................... 93
TABLEAU4. 4: COMPARAISON LES EFFORTS DE CISAILLEMENT STATIQUE POUR DIFFERENT CODE DE CALCUL. .... 95
TABLEAU4. 5:COMPARAISON LES EFFORTS DE CISAILLEMENT DYNAMIQUE POUR DIFFERENT CODE DE CALCUL. . 96

14
Notations et symboles utilisés
W : ouverture de fissure.
𝑙𝑠 : Longueur de glissement.
𝜀𝑠 : Déformation de l’acier suivant l’axe de la barre.
𝜀𝑐 : Déformation du béton.
𝜀𝑠𝑟1 : Déformation de l’acier au point de non glissement sous l’action des forces provoquant la
fissuration du béton.
𝜀𝑠𝑟2 : Déformation de l’acier au droit de la fissure.
εsm : Déformation moyenne des aciers sur la longueur de glissement.
εcm : Déformation moyenne du béton sur la longueur de glissement.
𝜀𝑐1 : Déformation du béton au point de glissement nul.
𝜀𝑠1 : Déformation de l’acier dans l’hypothèse d’absence de fissuration.
𝜀𝑠2 : Déformation de l’acier au droit de la fissure.
Δ𝜀𝑠𝑟 : variation de déformation entre le point de non glissement et la section fissurée.
𝜎𝑠2 : Représente la contrainte de l’acier au droit de la fissure.
𝜎𝑠𝑟 : Représente la contrainte de l’acier juste avant la fissuration.
𝑆𝑟,𝑚𝑎𝑥 : Espacement maximal des fissures.
𝜀𝑠𝑚 : Allongement moyen des armatures, sous la combinaison d’actions considérée, tenant
compte de la contribution du béton tendu.
𝜀𝑐𝑚 : Allongement moyen du béton entre les fissures.
𝛽𝑡𝑚 : Facteur permettant de passer de la déformation maximale de l’acier à la déformation
moyenne.
𝜎𝑠 : Contrainte de l’acier dépende de l’espacement maximal ou de diamétre maximale 𝛷𝑠 des
barres.
𝜎𝑠2 : Contrainte de l’acier au droit de la fissure.
𝜎𝑠𝑟 : Contrainte de l’acier à l’instant où la fissure est provoquée.
𝜌𝑝,𝑒𝑓𝑓 : Ratio d’armature calculée par rapport à une section dite «Effective » du béton.
𝐴𝑐,𝑒𝑓𝑓 : est l’aire de béton entourant l’armature tendue sur une hauteur (ℎ𝑐,𝑒𝑓𝑓 ).
𝛷𝑠 : Diamètre des barres.
c : Enrobage des armatures longitudinales.
𝑓𝑠 : Contrainte dans les aciers.
𝑓𝑦 : Limite élastique des aciers.
𝐸𝑠 : Module élasticité d’acier.
𝜎𝑠 : Contrainte d’acier au voisinage de la fissuration.
𝑓𝑐𝑡 : définit comme le produit entre taux armatures et la contrainte au voisinage de fissure.
𝑁𝑥𝑦 , 𝑁𝑥𝑥 , 𝑁𝑦𝑦 : Les contraintes principales en MPA.
𝜏𝑢 : Contrainte de cisaillement.
𝜎𝑐𝑚 : Contrainte moyenne de compression.
𝜎𝑛 : Contrainte normal.
𝑓𝑐𝑘 : Résistance du béton en compression à 28 jours.
𝑓𝑐𝑑 : Résistance de dimensionnement du béton en compression en MPa.
𝑉𝑅𝑑𝑖 : Résistance au glissement en MPa.
𝑓𝑐𝑡𝑚
: Rapport entre la contrainte caractéristique de béton et la contrainte de cisaillement.
𝜏
𝑏𝑚𝑠
𝑉𝑢 : Valeur de calcul de l’effort tranchant.

15
Introduction générale

La simulation du comportement linéaire et non-linéaire des structures en béton armé


soumises à des chargements sismiques est un problème d’actualité. Des séismes plus au
moins puissants (Kobe (Japon), Izmit (Turquie) ou Athènes (Grèce), …) ont montré que
même les constructions conformes aux codes parasismiques de dernière génération ne sont
pas à l’abri des catastrophes naturelles.
Les voiles sont des structures verticales à deux dimensions. Ils présentent généralement une
grande résistance et rigidité vis-à-vis des forces horizontales dans leurs plans. Ils sont très
présents dans les constructions. Ils jouent un rôle très important pour la sécurité sous
chargement sismique.
La dégradation que peuvent subir des murs voiles en béton armé est toujours susceptible
d’apparaitre. Plusieurs études de recherche ont été réalisées pour évaluer l’endommagent du
béton en termes d’ouverture de fissures. A l’heure actuelle, les codes de conception et de
dimensionnement proposent des formules empiriques, pour que l’ingénieur puisse obtenir des
informations sur les mécanismes de rupture des éléments et l’évolution des ouvertures des
fissures.
Les murs voiles soumis à un chargement sismique dans leurs plans. L’estimation des
ouvertures de fissures causées par le chargement de cisaillement induit par le séisme n’est pas
une tâche facile comparée aux cas de traction ou de compression simple. Les règlements
parasismiques ont généralement tendance à surestimer les ouvertures des fissures en utilisant
des coefficients de sécurité amplifiés et aussi, en négligeant quelques paramètres tels que les
effets d’échelle et les effets de vitesse. C’est la raison pour laquelle une coordination
internationale est donc nécessaire pour la compréhension et la prédiction du comportement
sismique des murs voiles en béton armé, cette dernière a pour objectif d’améliorer les codes
de calcul, et pour vérifier si les règles de conception estiment bien les paramètres locaux.
Les codes de conception cherchent à limiter la participation des ouvertures des fissures
uniquement par les contraintes des aciers en excluant les contraintes de béton, ce qui est vrai
uniquement dans le cas des sollicitations normales d’une poutre soumise à la flexion. Or dans
les murs voiles ; les sollicitations engendrées par le séisme sont très complexes.
Plusieurs travaux expérimentaux et numériques traitent le comportement des voiles sous
différentes sollicitations ont été réalisées. Parmi ces programmes ; le programme CASSBA
(Conception et Analyse Sismique des Structures en Béton Armé), CAMUS c’est la continuité
du programme CASSBA (conception et analyse des murs sous séisme), dans ce dernier, on a
testé des voiles élancés faiblement armé sur une table vibrante dans le laboratoire LNEC
(laboratoire national for civil engineerings) a Lisbonne, le but de ce programme est de
démontrer la légitimité de la conception eurocode8 et PS92 et aussi la prédiction du
comportement sismique des structures en béton armé. Le programme SAFE (Structure Armé
Faiblement Elancé) a été élaboré en 2000 en collaboration entre deux sociétés EDF et
COGEMA dans le laboratoire ELSA,JRC, ISPRA par des essais pseudo dynamique, l’objectif
de ce programme est de voir le comportement des structures courtes fortement armées qui
sont utilisées dans le dimensionnement des ouvrages nucléaires.
L’objectif de notre étude est de confronter la modélisation numérique à l’expérimentation et
aux règlements internationaux traitant le dimensionnement des murs voiles sous chargement
sismique. Des informations globales (Courbe de capacité) et d’autres, locales (Ouverture de

16
fissures) sont évaluées numériquement et comparées à l’expérimentation et aux codes de
conception.
Le présent travail comporte quatre parties :

Le premier chapitre est dédié à une étude bibliographique sur le comportement des murs voile
ssous sollicitation sismique et à la présentation des différentes échelles de modélisation. Le
deuxième chapitre expose quelques codes de conception et leurs méthodes de calcul des
efforts de cisaillement et l’évaluation des ouvertures des fissures.
Le troisième chapitre traite une modélisation de la maquette SAFE (Structure Armée
Faiblement Elancé) par différente approche avec le logiciel CAST3M, ces simulations ont
pour but d’étudier la capacité et les performances des modèles élasto-plastique en
dommageable à reproduire le comportement structurel des murs voiles en béton armé.
Le quatrième chapitre est consacré à l’évaluation des ouvertures des fissures par la procédure
OUVFISS dans la partie post traitement et d’évaluer les efforts de cisaillement, on comparant
avec la réglementation internationale. Finalement, des conclusions et des perspectives de ces
travaux sont présentées.

17
Chapitre 01

Les voiles en béton armé sous chargement sismique :


Comportement et principe de dimensionnement.

18
1.1. Introduction

L’utilisation des voiles en béton armé comme des éléments porteurs pour les
constructions parasismiques devient de plus en plus fréquente. La raison est que les murs de
contreventement, outre leur rôle porteur vis-à-vis des charges verticales et sont particularité
de prendre les efforts horizontal, ont une grande capacité de dissipation énergie lors d’un
événement sismique.
Les structures des murs en cisaillement (voile) ont été pratiquées depuis 1960pour les
immeubles moyens et grands hauteur et pour les ouvrages spéciaux (les ouvrages nucléaires)
dans de nombreux pays qui sont exposés au risque sismique.
Pour rendre les voiles comme des éléments de résistances et comme des éléments de
dissipation d’énergie dans les risques majeurs (séisme), il est primordial de connaitre leur
comportement linéaire et non linéaire en termes d’endommagement et d’ouverture des
fissures sous différentes excitations (sollicitations).
La figure 1.1 illustre quelque photo de modes de rupture des voiles en béton armé sous
sollicitation sismique.

Figure1. 1:Dommages subis par les voiles (a) Rupture de section à la base du
voile (Séisme de Boumerdes du 21/05/2003) ; (c) endommagement aux niveaux de
nœud poteaux forte poutre faible. (b-d) : Rupture en cisaillement des voiles
(Séisme de L'Aquila 06/04/2009.

1.2. Pourquoi utiliser des murs en cisaillement

Il y a plusieurs raisons d’utiliser les murs voiles comme des éléments résistants aux
efforts sismiques. Parmi ces raisons :
 Rendre la structure plus rigide.
 Diminuer les phénomènes d’instabilité.
 Rendre le comportement de la structure plus fiable par rapport aux structures auto
Stable.
 L’utilisation des voiles réduit l’influence des éléments non porteurs sur le
comportement de la structure, influence que nous ne maîtrisons pas aujourd’hui.

19
 Un voile fissuré garde une partie de sa résistance, ce qui n’est pas général le cas d’un
poteau.
 Réduire les déplacements horizontaux lors des séismes qui abaissent la conséquence
psychologique sur les habitants de haut bâtiment, et aussi réduire le cout de réparation
des dommages des éléments non porteurs.

1.3. Les caractéristiques essentielles influençant le comportement des murs voiles BA :

 Élancement qui peut être défini comme le rapport entre la hauteur et largeur du voile
(H/L).
 Les caractéristiques géométriques des éléments.
 Le pourcentage et la disposition des armatures (ex armatures horizontales empêchent
l’ouverture des fissures).
 Intensité de l’effort normal et de l’effort de cisaillement induit par le séisme.
 Type de chargement : statique ou cyclique, monotone.

1.4. Mode de fonctionnement des murs voiles

Les règlements parasismiques dans le monde classifient les voiles selon leur
élancement on distingue 2 types des voiles :
Les voiles élancés : qui ont un élancement de l’ordre supérieur de 1,5.
Les voiles courts : qui ont un élancement de l’ordre inférieur de 1,5.
Pour les deux types, les voiles ont un comportement différent et un mode de rupture
Différent vis-à-vis des sollicitations.

1.4.1. Mode de fonctionnements des murs voiles élancés

Les modes de ruine des voiles élancés sont représentés sur la figure 1.2. Chaque mode
de rupture et discuté brièvement dans la partie suivant :

1.4.1.1. Endommagements en flexion

a. Mode F1 : rupture par plastification des armatures verticales tendue et endommagement de


béton par compression. C’est le mécanisme le plus fiable qui correspond à la formation des
rotules plastique dans la partie inférieure du voile avec une importante dissipation d’énergie.
Ce mode de ruine apparait dans les voiles plus élancés soumis à des efforts de cisaillement
modérés et un effort de compression faible.
b. Mode F2 : rupture par écrasement du béton, ce mode de ruine est rencontré pour les voiles
assez fortement armé, sollicité par des efforts normaux plus importants.
c. Mode F3 : la ruine s’effectue par la rupture fragile des armatures verticales tendues c’est
un mode qui est rencontré dans les voiles faiblement armés aux extrémités. Pour améliorer la
notion de ductilité, il faut concentrer les armatures dans l’extrémité.

1.4.1.2. Endommagement en flexion/cisaillement

Mode F/T : ce mode de ruine consiste à la plastification des armatures verticales de flexion
et armatures transversales, c’est ce qui se produit pour les voiles qui sont sollicités par des
efforts de flexion plus au moins importants ou avec insuffisance des armatures horizontales.

1.4.1.3. Endommagement en cisaillement

20
Mode T : la rupture se fait par écrasement de béton de l’âme, ce mode de rupture est
rencontré dans le cas des efforts de cisaillement plus important.

Figure1. 2::Modes de rupture de voiles élancés (Davidovici et al, 1985).

1.4.2. Mode de fonctionnement des murs voiles courts


Dans les voiles courts, on distingue 3 modes de ruine :
Mode T1 : Endommagement par glissement (Slinding shear) : La rupture par glissement à
l’encastrement est obtenue par plastification progressive des armatures verticales. Ce type de
ruine engendre une perte importante sur la rigidité et de la dissipation énergie.
Mode T2 : Endommagement par plastification des aciers horizontaux ou verticaux (tension-
failure) : Rupture par effort tranchant et rupture des armatures le long de fissure diagonale
c’est un cas aussi fréquemment rencontré.
Mode T3 : Endommagement par ruine par compression du béton de la diagonale
(compression-failure) : Rupture par écrasement du béton de l’âme, à la base des bielles
transmettant les efforts de compression. C’est un mode de ruine caractéristique des voiles
fortement armés, surtout s’ils sont associés à des raidisseurs sur leur bord.

1.5. Le principe de dimensionnement dans les codes de calcul :


La plupart des codes de construction parasismique dans le monde essaient de
dimensionner des bâtiments résistant aux séismes de façon à ce que l’énergie sismique doive
être absorbée et dissipée par des déformations post-élastique (plastique,
endommagement, …).

21
Figure1. 3: Modes de rupture de voiles courts (Davidovici et al. 1985).
Le calcul d’un voile court en flexion ne peut plus être basé sur l’hypothèse de la planéité des
sections. Si l’on veut assimiler son comportement à celui d’une poutre, il faut enrichir la
cinématique de la poutre en la dotant d’une distribution non-linéaire des déformations. Dans
la plupart des cas le calcul est effectué en utilisant des éléments plaques (Panagiotis Kotronis,
2006).
Les codes réglementaires parasismiques préconisent des coefficients du comportement q
inférieurs pour des structures à voiles par rapport aux structures à portiques, malgré
l’influence bénéfique des voiles sur leur comportement. La raison en est la volonté d’éviter
des ruptures fragiles. Alors que les voiles pour lesquels la flexion et prépondérante présentent
en général une grande capacité de dissipation d’énergie (Figure 1.4), les voiles cisaillés ne
sont pas suffisamment ductiles et les courbes efforts tranchant/cisaillement sont pincées
(Figure 1.5). Pour ces voiles, il n’y a pas de mode fondamentalement ductile, à moins de
dispositions d’armatures tout à fait spécifiques.
Aujourd’hui nous savons pourtant comment dimensionner une structure pour éloigner la
possibilité d’une rupture par effort tranchant. La tendance actuelle est donc pour une
augmentation des coefficients de comportement dans les codes réglementaires.

Figure1. 4:: Comportement ductile d’un voile élancé. Rupture due à la flexion.

22
Figure1. 5: Comportement non ductile d’un voile. Rupture due à l’effort tranchant.

1.6. Principe de dimensionnement des murs voile


Dans tous les règlements parasismiques, deux principes de dimensionnement sont
retenus :

1.6.1. Dimensionnement en capacité


Ce principe, utilisé dans l’eurocode 8, qui prévoit la formation d’un mode de rupture
par création des rotules plastique à la partie inférieure du voile qui est supposée encastré, une
grande capacité de ductilité qui s’organise dans cette région, pour assurer la dissipation
d’énergie il faut éviter la rupture fragile par effort tranchant en assurant de la plastification de
la zone critique et augmenter la ductilité par présence des cadres disposés autour des
armatures verticales.

1.6.2. Dimensionnement par création des rotules plastique


Il existe d’autre conception de dimensionnement sismique par exemple des
réglementations française PS92. Le fonctionnement de voile doit assurer la dissipation
d’énergie par endommagement de béton (ouverture des fissures) et plastification des aciers
en transformation d’énergie par soulèvement des masses, ces voiles devraient avoir un
fonctionnement multi-fusible résultant en une ductilité répartie et une fissuration distribuée
sur toute la hauteur.

Figure1. 6::principe de dimensionnement selon les règlements (EC8 vs PS92)[Xuan


Huy Nguyen ,2006].

23
1.7. Les niveaux de modélisation
La modélisation des structures passe d’abord par le choix d’un niveau de modélisation
et d’une échelle de discrétisation. Nous pouvons distinguer quatre niveaux de modélisation
[MESTAT et al 1995] :
1.7.1. Le niveau géologique : il vise à traiter une structure dans un environnement naturel par
référence à des données géologiques. Les principales caractéristiques de l’ouvrage sont prises
en compte même si sa géométrie peut être simplifiée. Le but de ce type de modélisation et de
déterminer les déformations du sol et celles de l’ouvrage. Ce niveau de modélisation est bien
adapté à la simulation des ouvrages géotechniques comme les ouvrages des soutènements, les
barrages-poids ….
1.7.2. Le niveau global : il vise à traiter une structure dans son ensemble, comme un bâtiment
modélisé par des éléments poutres ou plaques et/ou coques par exemple (structure benchmark,
structure Nice).
1.7.3. Le niveau semi - local : il correspond à l’étude d’un élément de structure. Pour un
bâtiment, c’est le comportement des planchers soumis à des charges localisées qui peut être
visé par exemple, les voiles sollicités aux cisaillements.
1.7.4. Le niveau local : il correspond à l’étude d’une partie détaillée de la structure. Pour un
bâtiment, il s’agit des parties dont les dimensions sont faibles comparés aux éléments.

1.8. Les échelles de modélisation


La modélisation numérique est un outil incontournable dans le domaine de l’analyse et
de l'étude des structures/ouvrages de génie civil. Elle contribue à la compréhension du
comportement mécanique des matériaux utilisés dans les constructions afin de maîtriser les
mécanismes entraînant leur dégradation et permettre ainsi la prédiction de leur comportement
sous diverses sollicitations.
Plusieurs approches peuvent être adoptées afin de modéliser le comportement d’un élément
de structure soumise à des chargements mécaniques : on peut citer trois approches : approche
globale, approche semi-globale et approche macroscopique et ils existent d’autres approches
de modélisation voir la figure 1.7.1

Figure1. 7: Modélisation SMART2008 : Trois types d'approches (a) Approche globale


; (b-c) Approche locale ; (d) Approche semi-local (Lermitte 2010).

1.8.1. Approche globale


Le principe de cette approche consiste à intégrer les lois de comportement global qui
sont des relations entre les variables généralisées. Les déformations généralisées (εx, εy, εz,
γx,γy, γz) et les contraintes associées (Mx, My, Mz, Ty, Tz ou Nx). Ce niveau bien adapté à
l’étude des structures complexe sous chargement sismique, la formulation des lois de
comportement s’appuie sur la théorie classique de la plasticité avec introduction de seuil de
plasticité f (N, V, M)=0 ou la mécanique d’endommagement par association d’une variable

24
interne qui présente la ruine de matériaux. Ce type de modélisation bien avancé est intégrée
dans les nouveaux codes de calcul par exemple EC8.
Concernant la modélisation des éléments de structure, l’état de l’art des modèles globaux
développés pour le calcul des voiles est moins avancé que celui des modèles construits pour le
calcul des poutres et des poteaux. En effet, dans le cas des voiles en béton armé, les
déformations inélastiques dues à l’effort tranchant peuvent avoir une influence notable sur la
réponse globale. Une première approche pour prendre en compte ces effets réside dans
l’enrichissement des modèles de flexion. Ainsi, la loi de flexion de Takeda a été modifiée par
(Roufaiel and Meyer (1987)) afin de décrire le pincement des boucles d’hystérésis associée à
la fissuration diagonale par cisaillement. [Khuong LENGUYEN, 2015]
Le comportement de matériaux ainsi que les mécanismes de formation d’ouverture et
refermeture des fissures ne peuvent pas être considérés.
Parmi les avantages d’approche globale

 Elle est adaptée à la modélisation des structures minces de type coque (dalle et voile).
 elle est peu coûteuse en temps CPU et en taille mémoire.
 elle peut être combinée à l’emploi d’autres éléments de structure (coques
multicouches, …).

Les inconvénients d’approche globale


 Elle ne permet pas de représenter finement les phénomènes mécaniques et les
réponses locales.
 la mise en œuvre de modèles de comportement globaux nécessite une identification
des paramètres pouvant être difficile.

1.8.2. Approche macroscopique


Ce type de modélisation basé sur la théorie de thermodynamique et l’approche locale
très utile pour la validation du comportement expérimental d’un élément de structure et pour
confirmer les modèles globaux. Les lois de comportement des matériaux sont entièrement
décrites par des variables locales (contrainte –déformation), ce type de modélisation est en
général indépendant de la géométrie de la structure. Ce type de modélisation donne des
informations locales sur l’état de la fissuration et de l’endommagement. L’analyse est
relativement lourde, prendre un temps plus important à cause de l’élaboration de modèle
théorique et implémentation dans le code de calcul par le modélisateur.
Parmi les avantages et les inconvénients d’approche locale.
 elle permet de représenter finement des géométries complexes du type nœuds
d’ossatures ou zones que l’on cherche à modéliser avec précision (y compris
l’ensemble du ferraillage longitudinal et transverse).
 elle permet de représenter tous les types de chargement.
 elle permet d’accéder aux grandeurs globales et locales.
 elle est plus coûteuse en temps CPU et en taille mémoire que les modélisations semi-
globales et globales.

1.8.3. Approche semi-globale


Cette approche est un compromis entre les 2 approches précédentes, cette dernière
permet d’utiliser le modèle local du comportement du béton et acier, dans le cadre de la
cinématique simplifiée qui permet un calcul des variables statiques locales (contrainte) qui
sont intégrées pour les variables statiques généralisées, et donne une idée sur les variables
internes tel que l’endommagement du béton et les ouvertures des fissures.
La modélisation des poutres a donné lieu au développement d’éléments multicouches par La
Borderie (1991) et Merabet (1990) pour des analyses bidimensionnelles et d’éléments
multifibres par (Kotronis and Mazars 2005; Mazars et al. 2006) pour des analyses
tridimensionnelles. L’avantage important de ces approches réside dans le couplage implicite

25
des efforts de flexion et de l’effort normal. Cette approche est la mieux adaptée pour la
modélisation des systèmes de poutres et poteaux. Dans ce cadre, elle a permis d’obtenir
d’excellents résultats, en statique comme en dynamique. Dans le cas des voiles relativement
élancés, l’approche semi-locale (semi-globale) peut apporter des résultats intéressants, à
condition que la perturbation apportée par l’effort tranchant ne soit pas très importante. Dans
le cas des voiles faiblement élancés, Kotronis and Mazars (2005) ont proposé un modèle
d'endommagement continu et différente stratégie numériques simplifiés pour simuler le
comportement de ces murs soumis à des mouvements sismiques. Pour la modélisation 2D des
murs en béton armé dont le comportement est contrôlé principalement par la flexion, un
élément de poutre d’Euler est adopté. Pour les problèmes 3D, un élément de poutre de
Timoshenko multifibres ayant des fonctions d'interpolation d'ordre supérieure a été
développé. Enfin, pour décrire le comportement des murs à faible élancement, les auteurs
proposent un modèle béton armé équivalent (ERC model). La comparaison des résultats
numériques avec les résultats expérimentaux des murs testés sur une table vibrante montre les
avantages, mais aussi les limites de l'approche adoptée [Khuong LENGUYEN, 2015]

1.8.3.1. Approche multifibre


L’approche multifibres s’appuie sur une cinématique du type poutre Timoshenko ou
Euler-Bernoulli. Un élément de structure présente donc deux niveaux de discrétisation. Le
premier niveau est réalisé à travers des éléments poutre et, le second niveau consiste à
discrétiser chaque section selon un certain nombre de fibres. Ainsi, au premier niveau, il est
possible de définir les déformations généralisées (effort normal, efforts tranchants, moments
fléchissant), le lien entre ces deux quantités étant fait au moyen de termes de rigidité lorsque
le comportement est élastique. Elle permet aussi de nous donner idée claire sur les
déformations locales dans une structure en termes d’endommagement et ouverture des
fissures.

Figure1. 8:Principe de fonctionnement de l’élément poutre multifibre [Cédric G.2006].

Les avantages et inconvénients de modélisation multifibre :


 Elle est adaptée à la modélisation des éléments de type poutre.
 Moins coûteuses dans le temps exécution et la mémoire de stockage.
 Elle permet de réduire la taille de problème par rapport à un modèle en élément massif
(approche globale).
 Ce type modélisation il est adaptable pour des structures simples.
 C'est un modèle robuste pour définir le comportement local de béton.

Inconvénient de cette approche ne permet pas de représenté le ferraillage transversal de la


structure.

1.8.3.2. Approche multicouche


Une plaque multicouche est un ensemble des éléments finis basé sur la théorie des
plaques. Ce niveau de modélisation consiste à découper par des couches de béton armé
suivant la troisième dimension (épaisseur h).

26
Contrairement à la modélisation multifibre qui associe les lois de comportement pour chaque
fibre, ce type de modélisation associe les lois de comportement local en terme (contrainte
déformation) pour différents matériaux.

Dans le cadre d’une modélisation multicouche, il y a deux niveaux d’intégrations. Il y a la


modélisation dite «plane» qui sera représentée par des cellules (de support géométrique plan
avec trois points de Hammer sur la surface) et une modélisation suivant l’épaisseur
(perpendiculaire au support géométrique, discrétisée en couches).

1.8.3.2. Approche multicouche

Figure1. 9::Principe de fonctionnement de l’élément poutre multicouche. [Mohammed


Ibrahim 2003].
Parmi les avantages d’une modélisation multicouche :
 C’est une solution très fiable et robuste pour l’analyse non linéaire des structures pour
différent chargement.
 Elle est adaptée à la modélisation des structures minces de type (voile, dalle,
fondation).
 La réduction des problèmes (passage global vers un passage semi-global) en termes de
modèle (passage à échelle massif vers une structure simple).
 Rapidité dans le temps d’exécution.
 Le seul inconvénient elle ne permet pas de représente le ferraillage transversal de
manière fiable.

27
1.9. Conclusion

Le comportement des murs voiles soumis à des efforts tranchants est un problème
actualité dans le monde entier. Pour comprendre le comportement linéaire et non linéaire des
voiles soumis à des sollicitations sévères plusieurs programmes ont été élaborés pour la
prédiction de comportement des voiles au chargement sismique et améliorer les codes de
calcul sismique (développement des formules pour les ouvertures des fissures et effort de
cisaillement statique et dynamique). Le chapitre suivant est dédié à une exposition des
différents règlements traitant les ouvertures de fissures et les efforts de cisaillement dans un
mur voile soumis à un chargement sismique.

28
Chapitre 02

Codification des ouvertures des fissures et des efforts de cisaillement selon


la réglementation.

29
2.1. Introduction

La fissuration dans les ouvrages en béton armé est un processus très complexe que
nous maîtrisons très peu d’aujourd’hui. En effet la plupart des réglementations limitent les
contraintes des aciers sans tenir de la contribution de béton pour la résistance aux différentes
sollicitations.
Les formules actuelles qui sont élaborées par les codes de calcul pour estimation des
ouvertures des fissures sont exprimées en fonction des contraintes des aciers aux droits des
fissures, la formule est bien adaptée pour des éléments qui sont soumis aux sollicitations en
flexion ou effort normal tel que les poutres et les poteaux, or ceci n’est plus le cas pour les
voiles soumis aux cisaillements (cisaillement statique et dynamique) ou le béton contribue
dans la résistance aux efforts tranchants.
Dans ce chapitre on va présenter brièvement les différents codes qui traitent le calcul d’effort
tranchant et l’évaluation des dommages en termes d’ouverture des fissures.

2.2. Le processus de fissuration dans les éléments en Béton armé

Dans une structure en béton armé, les fissures peuvent apparaitre très tôt. Elles
prennent tout d’abord la forme des microfissures (due à la composition de béton) qui se
propagent lorsque le chargement devient important et que la limite en traction est dépassée,
des macro-fissures peuvent apparaitre et se développer jusqu’à la rupture de la structure.

2.2.1. Formation des fissures

Pour qu’il y ait une fissure, il faut un glissement de l’acier dans le béton sur une
longueur qui s’appelle longueur de glissement 𝑙𝑠 , ce glissement dû à la différence entre la
déformation de l’acier et le béton dans une longueur, l’ensemble de ce terme est variable le
long de la barre et la durée de chargement, en fonction du temps.

2.2.2. Ouverture des fissures

L’étude de la fissuration du béton a débuté par celle des tirants du fait de la simplicité
apparente de ces éléments travaillant dans une seule direction. Ce sont les résultats de ces
premières études qui ont donné lieu aux premières formulations codifiées dans la
réglementation française sur la fissuration, c'est-à-dire celles figurant dans les règles CCBA
68. La formule proposée pour ouverture des fissures c’est la différence entre déformations
acier et le béton multiplié par la longueur de glissement.

W = 𝒍𝒔 *(𝜺𝒔𝒎 -𝜺𝒄𝒎 ) (2.1)

Avec :
𝑙𝑠 : Longueur de glissement.
𝜀𝑠𝑚 : Déformation moyenne des aciers sur la longueur de glissement.
𝜀𝑐𝑚 : Déformation moyenne du béton sur la longueur de glissement.
Ces deux grandeurs physiquement simples sont en fait non accessibles par un calcul
direct et d’autant moins qu’il faut distinguer deux cas très différents, celui de la fissuration
(accidentelle) et celui de la fissuration (systématique), ces deux notions explicitées ci-dessous.

2.2.2.1. Cas de la fissuration accidentelle

Ce cas donne lieu au schéma de déformation de l’acier en fonction de son abscisse sur
la barre pouvant être représenté sur la figure (2.1) :

30
Figure2. 1:présentation de la déformation d’acier (haut) et de béton (bas) en fonction
de la longueur de la barre dans le cas de fissuration accidentelle. [Etienne Gallitre,
2008].
La fissuration accidentelle présente une particularité par apport à la fissuration
systématique, un point de glissement nul c.-à-d.la déformation de béton égale à celle de
l’acier.
Avec :
εsr2: Déformation de l’acier au droit de la fissure (grandeur accessible dans un calcul béton
armé) conventionnel dans le cas des poutres ou des tirants.
εc1 : Déformation du béton au point de glissement nul : grandeur accessible avec certaines
précautions (avant fissuration systématique), puisque dans ce cas il s’agit de la déformation
du béton armé considéré en dehors de la zone fissurée.
εsr1 : Déformation de l’acier au point de glissement nul.
εcr1: Déformation du béton au point de glissement nul.
Δεsr: Différentiel maximal de déformation de l’acier le long de la barre (différence de deux
grandeurs accessibles).
La grandeur [βt,m , Δεsr], permettant de passer de la déformation maximale à la déformation
moyenne est accessible par l’intermédiaire de ce coefficient βtm égal à un facteur
d’intégration de la déformation de l’acier le long de la barre. Ce coefficient est voisin de 0,5,
valant 0,6 compte tenu du caractère convexe de la courbe représentant εs le long de la barre

Dans le cas des poutres ou des tirants, où on a accès à 𝜀𝑠2 , puisque l’on considère que tout
l’effort de traction passe par les seuls aciers, on a
𝜺𝒔𝒎 = 𝜺𝒔𝟐 - βt,m (εsr2 - εsr1) (2.2)

2.2.2.2. Cas de la fissuration systématique :

Lorsque la fissuration est systématique, on obtient le schéma donné par la (figure


2.2) :

31
Figure2. 2: Représentation des déformations de l’acier et du béton avec l’axe de la
barre d’acier en abscisse dans le cas de la fissuration systématique [Etienne Gallitre,
2008].

Toutes les fissures sont apparues et les déformations de la structure se traduisent


essentiellement par une augmentation de la déformation des aciers et de l’ouverture des
fissures ; la déformation du béton entre les fissures est quant à elle plafonnée par sa capacité
de déformation (les pics de déformation du béton ne sont pas identiques).

Lorsque les déformations vont encore augmenter, l’allongement des aciers croit jusqu’à leur
limite élastique et au-delà : la raideur correspondante va se rapprocher de la raideur des aciers
seuls. Mais ces déformations « moyennes » des aciers comportent une ordonnée à l’origine
non nulle correspondant aux restrictions de déformations dans les zones de glissement nul.
Il est d’usage de représenter l’ensemble de ces phénomènes sur le diagramme de la figure
(2.3) visualisant la déformation moyenne des aciers qui constitue la donnée de base pour
estimer l’ouverture de fissure.

Figure2. 3: contrainte d’acier à la fonction de leur déformation moyen [Etienne


Gallitre, 2008].

32
Avec :
𝜎𝑠2 : représente la contrainte de l’acier au droit de la fissure.
𝜎𝑠𝑟 : représente la contrainte de l’acier juste avant la fissuration.

Ce phénomène est appelé “stiffening effect“ : au fur et à mesure que les déformations
augmentent, l’effet relatif du béton autour des aciers diminue, car le glissement augmente et
l’effet restrictif du béton diminue.

Le stade initial I correspond à l’absence de fissures. Le stade II intermédiaire correspond à la


phase d’apparition de la première fissure jusqu’au début de la fissuration systématique. Les
fissures apparaissent les unes après les autres, avec des longueurs de glissement dont
l’extrémité à glissement nulle (entre les fissures) connaît des points où la déformation de
l’acier est égale à celle du béton. Le stade III correspond à la fissuration systématique avant
plastification des aciers (domaine de fonctionnement qui va nous intéresser). La totalité des
fissures est apparue et l’augmentation des déformations s’accompagne d’une augmentation de
la taille des fissures. De plus, au point de glissement nul la déformation des aciers n’est pas
nécessairement identique à celle du béton.

Le stade IV est celui de la plastification des aciers qui s’accompagne d’une diminution de la
raideur.

La distinction entre ces différents stades est fondamentale pour accéder aux grandeurs qui
nous occupent (principalement ls et εsm), Le passage du stade I au stade II s’opère par un
transfert de la totalité de l’effort sur la section d’acier seul. La possibilité de ce transfert est
conditionnée par la condition (très théorique) de non-fragilité. [Etienne Gallitre, 2008].

2.3. Contexte règlementaire et méthodes associées

2.3.1. Aperçu rapide sur les approches codifiées de limitation du dommage

2.3.1.1. Les codes basés sur des approches en force

Les approches en force correspondent à la culture des ingénieurs en général, elles


constituent donc la base de la plupart des codes à vocation normative. Mais en analysant la
réglementation BAEL, on observe tout de même que cette dernière en introduisant la notion
d’états limites a donné un critère en déformation pour l’état limite ultime (critère à 1 %
d’allongement des aciers par exemple, soit 5 fois la limite élastique conventionnelle). Cette
notion a bien entendu été reprise par les Eurocodes sur une même base, mais avec des critères
définis à partir de la capacité ultime des matériaux eux-mêmes. Il s’agit toutefois de critères
très locaux constituant en fait une sophistication de l’approche en force, compte tenu de la
nature de la sollicitation prise en compte.

On reviendra un peu plus loin dans l’analyse plus détaillée de chacun des règlements sur les
états limites de service qui constituent aussi une façon indirecte de limiter le dommage, qui
trouve une traduction en termes d’ouverture de fissure relativement explicite dans l’EC2.
[Etienne Gallitre, 2008].

33
2.3.1.2. Les codes basés sur des approches en déplacement

Se sont développés depuis une dizaine d’années, des codes basés sur des approches en
déplacement. Ces approches ne sont pas seulement utiles pour un objectif de non-
effondrement comme on pourrait le penser, certaines d’entre elles contiennent dans leurs
énoncés également une relation un peu plus explicite entre le niveau de déformation atteint et
le dommage attendu. Nous allons donc ci-après présenter très rapidement les critères
fonctionnels en déplacement indiqués dans les différents codes.

2.3.1.3. FEMA 356

Le code American FEMA 356 identifie le degré de dommage pour les murs voile en
fonction d’un coefficient qui s’appelle DRIFT (le déplacement différentiel entre les étages) le
tableau (2.1) présente les degrés de dommage des structures en voile en béton armé.

Niveau de Niveau Niveau non- Occupation Dommages


dommage opérabilité effondrement immédiate limités
La limite de Non défini Б <2 % Б<0.5% Б<1%
déplacement
entre les
étages(%)

Tableau2. 1: Tableau indiquant les limites des dommages selon le code de calcul FEMA 365.
2.3.1.4 AIJ

Le code parasismique Japonais définit les classes de dommage selon les ouvertures des
fissures et limitation des contraintes. Le tableau (2.2) donne des valeurs pour différentes
classes qui sont définies selon la réglementation japonaise :

Exigence Terme anglais Dommage Cracking Strain criteria


donné pour (except cracking description (%)
traduire description)
l’exigence
Non- « Safety » «no decay in No limited 1-2
effondrement lateral and axial
capacity »
Dommages « Restorability » « no crush W < 2 mm 0.8-1
limités of core »
Dommages très « Resparability » « no crush of W < 1 mm 0.4-0.6
limités cover »
Réversible « Serviceability » « no yielding » W < 0.2 mm 0.12-0.33

Tableau2. 2: indique les limites des dommages selon la réglementation japonaise.

34
Par ailleurs, les documents japonais analysés étaient également accompagnés d’une étude très
documentée constituant les « backgrounds documents » de ce code, lesquels comprenaient des
essais d’éléments de murs à l’échelle 1 et des études non linéaires associées c.-à-d. simulation
d’une structure réelle. [Etienne Gallitre, 2008].

2.3.1.5 ASCE/SEI

Le document (ASCE/SEI 43-05) donne les déformations admissibles en fonction de ces états
limites : (voir tableau 2.3) :

Limite Définition des limites déplacement


dommage entre les
étages (%)
A déformations permanentes importantes, proches de 0.75
l’effondrement.
B déformations permanentes modérées. 0.6
C déformations permanentes limitées. 0.4
D comportement essentiellement élastique. 0.4

Tableau2. 3: tableau donne les limites drift endommagement.

Ces états limites ne décrivent pas la fissuration associée, sachant que ces définitions sont
également valables pour les tuyauteries et les équipements. Ce document ouvre la voie des
calculs non linéaires, qu’ils soient linéaires équivalents (sur la base des coefficients
forfaitaires) ou explicites. Les analyses pseudo-statiques en poussée progressive (push-over)
sont autorisées pour les structures où il existe véritablement un mode fondamental. La notion
de comportement, essentiellement élastique, pose question pour un ouvrage en béton armé
dont on vient de voir dans l’exploitation de SAFE que la fissuration, source de perte de
rigidité, intervient très tôt et bien avant un drift de 0,4%. Nous interprétons cette notion
comme un comportement élastique du matériau acier seul, justifiant ainsi le mot
«essentiellement ». Nous reviendrons bien entendu sur cette question ouverte, qui dans le
domaine nucléaire, pourrait être interprétée comme le domaine de validité des modèles de
calculs permettant d’obtenir les spectres de planchers.

Pour compléter, il faut bien citer dans le présent sous-chapitre l’ATC-40 (compte tenu de la «
popularité » de ce texte provenant du récent développement de l’approche en déplacement
aux Etats Unis. Cette popularité provient du fait qu’il s’agit du seul document à caractère
normatif qui donne les outils véritablement opérationnels et simples pour mener jusqu’à son
terme une étude pseudo statique non linéaire avec correction d’amortissement. Cette méthode
concerne avant tout l’étude d’ouvrages «plastifiés ». En effet, ce dernier document ne trouve
son intérêt que si les courbes d’hystérésis reconstruites à partir de l’étude pseudo statique
exhibent une plastification significative permettant de corriger l’amortissement et donc de
réduire la sollicitation sismique. Or, cette plastification « facile » à obtenir pour une ossature à
partir de modèles à fibres est beaucoup plus difficile à exhiber pour un bureau d’études
industrielles qui ne dispose pas de modèle de comportement validé de voile en cisaillement
(par exemple de type bi linéaire voire tri linéaire).

2.3.1.6 ATC-40

Ce code présente les critères fonctionnels suivants pour les voiles fonctionnant en
cisaillement :

35
Type

d’éléments Eléments primaires


Eléments Secondaires

Exigence Occupation
Sécurité Stabilité Sécurité Stabilité
fonctionnelle Immédiate

déplacement
entre les
étages 0.4 0.6 0.75 0.75 1.5
(%)

Tableau2. 4: drift limite donné par ATC -40.


Les éléments primaires sont ceux qui sont nécessaires dans la résistance aux efforts
horizontaux, par opposition aux secondaires qui ne sont pas nécessaires.

2.3.1.7 Comparaison des différents codes et commentaires

Les différents critères sont comparés dans le tableau ci-dessous


Type

d’éléments Eléments primaires


Eléments Secondaires

Exigence Occupation
Sécurité Stabilité Sécurité Stabilité
fonctionnelle Immédiate

déplacement
entre les
étages 0.4 0.6 0.75 0.75 1.5
(%)

(*) Ce document ne donne pas directement les valeurs, mais les formulations permettent de
les obtenir.
(**) Suivant interprétation d’une synthèse du document japonais traduit en anglais et donnant
des exemples.

Tableau2. 5: synthèse des drifts limites donnés dans les codes.


D’une manière générale, l’écart entre FEMA/AIJ et ASCE/ATC du tableau s’explique selon
nous en partie par le fait que celui-ci concerne des voiles clairement identifiés comme
« courts ».

36
Si le 0,75 du ASCE/ SEI correspond à un point situé significativement du côté de la sécurité
pour l’état limite A, on peut par contre s’interroger sur la validité du 0,4 pour un état limite D,
la question portant évidemment sur la notion de réversibilité. En effet, les essais SAFE
montrent une quasi-fermeture des fissures en fin de séisme, mais une fréquence tout à fait
diminuée, notion il est vrai difficile à intégrer, puisque les fréquences chutent dés 0,1 % de
déplacement entre les étages.
L’absence de différenciation des états limite C et D dans le document, quel que soit le type de
voile reste donc surprenant. [Etienne Gallitre, 2008].

2.3.2. Codification des ouvertures des fissures

La plupart des codes de dimensionnement dans le monde expriment les ouvertures des
fissures en termes de contraintes (oudéformation) des aciers aux voisinages de la zone
fissurée, et ils ont négligé la contribution de béton dans la résistance aux ouvertures des
fissures.

2.3.2.1. Eurocode 2

L’Eurocode 2 estime également les ouvertures des fissures comme étant la différence
entre la déformation des aciers tendue et la déformation de béton multiplier par une longueur
appelée la longueur de glissement qui définit comme une distance entre les fissures.

𝐖𝒌 = 𝑺𝒓,𝒎𝒂𝒙 * (𝜺𝒔𝒎 –𝜺𝒄𝒎 ) (2.3)

Avec : 𝑆𝑟,𝑚𝑎𝑥 : Espacement maximal des fissures,


𝜀𝑠𝑚 : Allongement moyen des armatures, sous la combinaison d’actions considérée, tenant
compte de la contribution du béton tendu,
𝜀𝑐𝑚 : Allongement moyen du béton entre les fissures.
𝒇𝒄𝒕,𝒆𝒇𝒇 (𝟐. 𝟒)
𝝈𝒔 − 𝒌𝒕 𝝆 (𝟏 + 𝜶𝒆 𝝆𝒑,𝒆𝒇𝒇 )
𝒑,𝒆𝒇𝒇 𝝈𝒔
𝜺𝒔𝒎 − 𝜺𝒄𝒎 = ≥ 𝟎, 𝟔
𝑬𝒔 𝑬𝑺

Avec :
𝜎𝑠 : Contrainte de l’acier dépende de l’espacement maximal ou de diamètre maximal 𝛷𝑠 des
barres (voir tableaux 2.6 et 2.7)
𝑬 (𝟐. 𝟓)
𝜶𝒆 = 𝒔𝑬𝒄𝒎

Et𝜌𝑝,𝑒𝑓𝑓 ratio d’armature calculée par rapport à une section dite «Effective » du béton, 𝐾𝑡
valant 0,6 pour les chargements de courte durée et 0,4 pour les longues durées.
𝑨
𝝆𝒑,𝒆𝒇𝒇 = 𝒔 (2.6)
𝑨𝒄,𝒆𝒇𝒇

𝐴𝑐,𝑒𝑓𝑓 = b ℎ𝑐,𝑒𝑓𝑓 (2.7) est l’aire de béton entourant l’armature tendue sur une hauteur
(ℎ𝑐,𝑒𝑓𝑓 ).
(𝒉−𝒙) 𝒉 (𝟐. 𝟕)
𝒉 = Min[𝟐. 𝟓(𝒉 − 𝒅);
𝒄,𝒆𝒇𝒇 ; ]𝟑 𝟐

37
Cntainte
de l’acier
Diamètre maximal des barres
(MPa)
(mm)

𝑊𝑘 = 0.4 𝑊𝑘 = 0.3 𝑊𝑘 = 0.2


mm mm mm

160 40 32 25

200 32 25 16

240 20 16 12

280 16 12 8

320 12 10 6

360 10 8 5

400 8 6 4

450 6 5 -

Tableau2. 6:Diamètre maximal Φ_s des barres pour la maitrise de la fissuration


[S.Multon-2012].

contrainte
Espacement maximal des barres
de l’acier
(mm)
(MPa)

𝑊𝑘 = 0.4 𝑊𝑘 = 0.3 𝑊𝑘 = 0.2


mm mm mm

160 300 300 200

200 300 250 150

240 250 100 100

280 200 150 50

320 150 100 -

360 100 50 -

Tableau2. 7: Diamètre maximal Φ_s des barres pour la maitrise de la fissuration


[S.Multon-2012].

38
Figure2. 4: élément sollicité en traction.
-La longueur de glissement :
La distance entre fissures 𝑆𝑟,𝑚𝑎𝑥 est donnée par l’expression suivant :
𝑲𝟏 𝑲𝟐 𝑲𝟒 𝜱 ∗ (𝟐. 𝟖)
𝑺𝒓,𝒎𝒂𝒙 = 𝑲𝟑𝒄 +
𝝆𝒑,𝒆𝒇𝒇

Avec :
𝛷𝑠 : Diamètre des barres, s’il y a plusieurs barres de diamètres différents dans une même
section
𝟐.𝟗 𝟏𝟎(𝒉−𝒅) (𝟐. 𝟗)
𝜱∗ = 𝜱𝒔 * *
𝒇𝒄𝒕𝒆𝒇𝒇 𝒉

,𝛷𝑠 doit être remplacé par un diamètre équivalent :


𝒏𝟏 (𝜱𝟏 )²+ 𝒏𝟐 (𝜱𝟐) ² (𝟐. 𝟏𝟎)
𝜱𝒆𝒒 = 𝒏𝟏 𝜱𝟏 +𝒏𝟐 𝜱𝟐

c : enrobage des armatures longitudinales,


𝐾3 = 3,4 Pour des enrobages c supérieurs à 25 mm, le coefficient 𝐾3 est à remplacer par :
𝟐𝟓 𝟐 (2.11)
𝑲 =3,4
𝟑 .
𝑪 𝟑

𝑺𝒓,𝒎𝒂𝒙 = 1,3 (h – x) (2.12)

𝐾1 = 0,8 pour les aciers HA, coefficient fonction des propriétés d’adhérence des barres,
𝐾2 = 1pour la traction pure et 0,5 en flexion,
𝐾4 = 0,425.

2.3.2.2. Le code ACI (American Concrete Institute)

Les codes américains ne proposaient pas des formules pour le calcul ouverture des
fissures pour des éléments sollicité par cisaillement, mais propose des formules pour les
éléments soumis en traction et en flexion.
Flexion:
w=0.091𝒇𝒔 𝟑√𝒕𝒃 ∗ 𝑨𝜷 (𝒇𝒔 -5)*𝟏𝟎−𝟑 (𝟐. 𝟏𝟑)

39
Traction:
w=0.1𝒇𝒔 𝟑√𝒅𝒄 ∗ 𝑨*𝟏𝟎−𝟑 (𝟐. 𝟏𝟒)
.
Avec : W : ouvertures en inch.
𝑓𝑠 : Contrainte dans les aciers supposés reprendre tout l’effort de traction au droit de la fissure
en ksi.
𝑡𝑏 et 𝑑𝑐 : enrobage mesuré au centre des aciers en inch.
𝛽 : Rapport entre la distance axe neutre – face extérieure tendue et la distance axe neutre –
centre des aciers.
A : quantité définie comme l’aire de béton divisée par le nombre de barres d’aciers la
traversant en 𝑖𝑛𝑐ℎ2
Cette dernière quantité peut être remplacée par (π𝛷 2 /4)/ρ). Ce qui montre que l’expression
proposée n’est pas proportionnelle à Φ comme les règlements européens.

2.3.2.3. CEB FIP 1990 [CEB93] :


a. Ouvertures des fissures (𝑾𝒌 ) :
Il propose une formulation en fonction des contraintes dans les aciers au droit des
fissures après apparition de ces dernières.
La corrosion de l’acier n’est pas accrue par la formation des fissures, Donc : on doit avoir une
ouverture des fissures 𝑊𝑘 <𝑊𝑙𝑖𝑚
Avec : 𝑊𝑘 = 0.3 pour les éléments en béton armé sous la combinaison d’actions quasi-
permanente et pour toutes les classes d’exposition.
𝑾𝒌 = 𝒍𝒎𝒂𝒙 ∗ (𝜺𝒔𝒎 − 𝜺𝒄𝒎 + 𝜺𝒄𝒔 ) (𝟐. 𝟏𝟓)

𝑙𝑚𝑎𝑥 : « Longueur de glissement » égale à la distance sur laquelle il existe un glissement entre
l’acier et le béton dû à la mise en œuvre l’adhérence ; 𝑙𝑚𝑎𝑥 peut donc être considérée comme
étant le double de la longueur d’ancrage. Les déformations 𝜀𝑠𝑚 et 𝜀𝑐𝑚 sont les moyennes sur
cette longueur,
𝜀𝑐𝑠 : Déformation du retrait libre (en valeur absolue).
b. Contrainte de traction limite 𝝈𝒔,𝒍𝒊𝒎 dans les armatures :
𝟏 𝑬𝒔 (𝟐. 𝟏𝟔)
𝝈𝒔,𝒍𝒊𝒎 = ( ′ + ) 𝒇𝒕𝒎
𝝆 𝒇 𝑬𝒄𝒊

Pour 𝜌′ 𝑓 faible et un béton de haute résistance(𝐸𝑐𝑖 grand), on peut écrire


𝒇𝒕𝒎 (𝟐. 𝟏𝟕)
𝝈𝒔,𝒍𝒊𝒎 ≈ ′
𝝆𝒇

Si : 𝜎𝑠 ≤ 𝜎𝑙𝑖𝑚 : la formation d’une fissure isolée doit être envisagé,


Si : 𝜎𝑠 > 𝜎𝑙𝑖𝑚 : on admet que la fissuration est stabilisée.

𝑨 (𝟐. 𝟏𝟖)
𝝆′𝒇 =𝑨′ 𝒔
𝒄𝒇

𝜌′𝑓 : Ratio d’armature de la section d’enrobage.


Avec 𝐴′𝑐𝑓 : Représente l’aire efficace définie ci-dessus.

40
L’air efficace 𝐴′𝑐𝑓 tient compte de la distribution non uniforme des contraintes normales
(dans le béton) dues aux forces d’adhérence, au voisinage des extrémités de la longueur de
glissement.
𝐴′𝑐𝑓 : est définie par sa hauteur efficace ℎ′𝑓 (ou largeur 𝑏′𝑓 )
- Pour les dalles par exemples où les armatures de même diamètre sont disposées en un seul
lit :
𝜱 (𝒉−𝒙) (𝟐. 𝟏𝟗)
𝒉′𝒇 = min{𝟐. 𝟓 (𝒄 + ) ; }
𝟐 𝟑

Où h est la hauteur totale et x la hauteur comprimée.


- Pour un élément en traction (voiles, armatures de peau d’une âme de poutre,…) :
𝜱 𝒃𝟎 (𝟐. 𝟐𝟎)
𝒃′𝒇 = min{𝟐. 𝟓 (𝒄 + ) ; }
𝟐 𝟐

Où 𝑏0 est l’épaisseur totale de l’élément considéré.

c. Espacement des fissures – longueur de glissement :


c.1. dans le cas de la fissuration stabilisée : (systématique) :
𝜱 (𝟐. 𝟐𝟏)
𝒍𝒎𝒂𝒙 =𝟑.𝟔𝝆′
𝒇

c.2. dans le cas de la fissuration unique : (accidentelle) : la longueur de glissement


(double de la longueur d’encrage) :
𝝈 𝜱 (𝟐. 𝟐𝟐)
𝒍𝒎𝒂𝒙 = 𝒔 𝑬
𝟐𝝉𝒃𝒌 𝒔 𝝆′
𝑬𝒄𝒊 𝒇

Où 𝜏𝑏𝑘 = Φ. 𝑓𝑡𝑚 : fractile inférieur de la contrainte d’adhérence moyenne,


Φ = 1.8 au premier chargement,
Φ = 1.35 pour un chargement permanent ou répété.

d. Allongement moyen :
-cette expression est valable pour une fissure unique ou pour la fissuration stabilisée :
𝝈
𝜺 - 𝜺 = 𝒔 (𝟏 − 𝜷′ 𝒔,𝒍𝒊𝒎 )
𝝈 (𝟐. 𝟐𝟑)
𝒔𝒎 𝒄𝒎 𝑬𝒔 𝟐 𝝈𝒔

Avec :
𝛽′2 = 0.60 pour le cas général,
𝛽′2 = 0.38 pour la fissuration stabilisée, sous l’effet des charges permanentes ou répétées.
Avec :
𝒇𝒕𝒎 (𝟐. 𝟐𝟒).
𝜺𝒄𝒎 ≈ 𝜷′𝟐
𝑬𝒄𝒊

2.3.2.4. Norme SIA 262

La norme SIA 262 ils ont exprimé l’ouverture des fissures en fonction des contraintes
des aciers en voisinage de la fissure et ils ont négligé la contribution des contraintes de béton.
𝜱 𝝈𝒔 ² (2.25)
ω= 𝟖 𝒇𝒄𝒕 𝑬𝒔

41
Avec :
𝐸𝑠 𝛷Diamètre d’armature.
: Module élasticité d’acier.
𝜎𝑠 : Contrainte d’acier au voisinage de la fissuration.
𝑓𝑐𝑡 : La résistance caractéristique de béton.
ω : Ouverture des fissures.

2.3.2.5. Fib model concrete (MC2010)

Ce document propose une formule pour estimation des ouvertures des fissures en
fonction des contraintes principale et la contrainte de cisaillement au voisinage de la
fissuration.

𝝎 = 2 𝒍𝒔𝒎𝒂𝒙 *(𝜺𝒔 − 𝜺𝒄 ) (2.26)

Avec :

𝟏𝒇 𝜱𝒔
𝒍𝒔𝒎𝒂𝒙 = K*c+𝟒 𝝉 𝒄𝒕𝒎 𝝆 (2.27)
𝒃𝒎𝒔 𝒔,𝒆𝒇𝒇

(𝜺𝒔 − 𝜺𝒄 ) = 0.339*𝟏𝟎−𝟑 *(cos (Ɵ)+sin (Ɵ))-0.264*𝟏𝟎−𝟒 (2.28)

(2.29)

𝟏 𝑵𝒙𝒚
𝜽 = 𝟐arctg (2abs(𝑵 ))
𝒙𝒙 −𝑵𝒚𝒚

𝑙𝑠𝑚𝑎𝑥 : Longueur de glissement.


K : coefficient égale 1.

42
C : enrobage.
𝑓𝑐𝑡𝑚
: Rapport entre la contrainte caractéristique de béton et la contrainte de cisaillement. des
𝜏𝑏𝑚𝑠
armatures égale 1.8.
𝛷𝑠 : Armature acier.

𝜌𝑠,𝑒𝑓𝑓 : Taux d’armature vertical.

𝜀𝑠 : Déformation des aciers.


𝜀𝑐 : Déformation de béton.
𝑁𝑥𝑦 , 𝑁𝑥𝑥 , 𝑁𝑦𝑦 : Les contraintes principales en MPA.
𝜽: angle entre les armatures vertical et la direction des contraintes principales de traction.

2.3.3. Les approches expérimentales pour l’estimation des ouvertures des fissures

De nombreux articles ont été publiés sur ce sujet avec comme objectif d’améliorer les
codes de calcul parmi ces approches :

2.3.3.1. Approche Chi et Kirstein

Il a proposé une formule qui dépend de la contribution de la contrainte d’acier au


droit de la fissuration et le diamètre armature.

𝟓𝝉𝜱 𝟒𝟑𝟖 (2.30)


𝑾𝒔 = (𝒇𝒔 - 𝝉𝜱 )
𝑬𝒔

Avec :

𝑊𝑠 : Ouverture des fissures estimées.


𝐸𝑠 : Module élasticité de acier.
Φ : Diamètre armature d’acier.
𝑓𝑠 : Contrainte d’acier à la droite de fissuration.
𝜏 : Coefficient variant de 0.786 et 1.

2.3.3.2. Approche Etienne Gallitre

D’après les essais de programme SAFE (Structure Armée Faiblement Elancé), Mr


Etienne G (2008) a réussi à proposer une formule d’ouverture des fissures en fonction de la
distorsion mesurée de la forme suivante :

43
𝜱 (2.31)
W = K (𝒇𝟐𝟖 )ɣ
𝝆

Avec : Φ diamètre moyen, k (f28) coefficient fonction de f28 et ρ ratio d’armatures. Cette loi
s’est vérifiée sur une plage très étendue, allant du début de la fissuration jusqu’à la ruine.

2.3.4. La codification des efforts de cisaillement


2.3.4.1. Code BAEL

Le BAEL est basé sur la théorie des états-limites ultimes et services, mais la résistance
à l’effort tranchant n’est exprimée qu’au travers d’une condition d’état limite ultime.
- Le BAEL 91 ne propose pas de règle spécifique pour les voiles en béton armé, mais ils sont
en général considérés comme des poutres soumises à un effort tranchant en satisfaisant les
deux conditions suivantes :

Béton :
𝒇𝒄𝒋 (2.32)
𝝉𝒖 = 𝑴𝒊𝒏 (𝟎. 𝟐 , 𝟓𝑴𝒑𝒂)
𝜸𝒃

Avec :
𝑽𝒖 (2.33)
𝝉𝒖 =
𝒃𝟎 ∗ 𝒅

𝑉𝑢 : Valeur de calcul de l’effort tranchant.


𝑏0 : Largeur de l’âme.
𝑑 : Hauteur utile de la poutre.
𝑓𝑐𝑗 : Résistance caractéristique du béton en compression.
𝛾𝑏 = 1.15 coefficient partiel de sécurité du béton (accidentel).

Armatures horizontales :
𝝉𝒖 = 0.9𝝆𝒉 𝒇𝒆 + 0.3K𝒇𝒕𝒋 (2.34)

Où 0.9𝜌ℎ 𝑓𝑒 représente la contribution de l’acier et 0.3K𝑓𝑡𝑗 celle du béton.

Avec :𝑓𝑡𝑗 Résistance caractéristique du béton en traction,


𝑓𝑒 Limite élastique de l’acier,
𝝈
K = 1 + 3 𝒇𝒄𝒎 (2.35)
𝒄𝒋

𝜎𝑐𝑚 Contrainte moyenne de compression,


𝑨𝒕 (2.36)
𝝆𝒉 = 𝒃
𝟎 𝑺𝒕

44
𝐴𝑡 Section d’un cours d’armatures horizontales,
𝑆𝑡 Espacement des cours successifs d’armatures horizontales,
ɣ𝑠 = 1 coefficient partiel de sécurité de l’acier (accidentel)

Les armatures verticales : sont justifiées par la règle dite des ‘coutures’ qui s’exprime de la
façon suivante :

𝝉𝒖 = 𝝈𝒖 + 𝝆𝒉 𝒇𝒆 (2.37)

Où 𝜎𝑢 est la contrainte normale, sa valeur est considérée ici positive en compression.

2.3.4.2. PS 92

C’est un code parasismique français qui s’applique aux constructions nouvelles


relevant du risque dit « normal », donc pas a priori aux constructions nucléaires relevant du
risque « spécial ».

- Conventions:

𝑓𝑡 : Résistance en traction du béton,


𝑓𝑐 : Résistance en compression du béton, (Pour le béton, le coefficient de sécurité
ɣ𝑏 Est pris égal à 1.15 et pas 1.5
𝑓𝑦 : Limite élastique des aciers,
𝜎𝑛 : Contrainte normal,
b: longueur du mur,
a: largeur du mur (panneau central),
𝑆𝑤𝑒𝑏 : Aire du panneau central,
𝐴𝑤𝑒𝑏 : Section totale d'acier dans le panneau central,
𝑆𝑓𝑙 : Aire d'une aile,
𝐴𝑓𝑙 : Section totale d'acier dans chaque aile.

𝑨𝒘𝒆𝒃 (𝟐. 𝟑𝟖)


𝝆𝒘𝒆𝒃 =
𝑺𝒘𝒆𝒃

𝑨𝒇𝒍 (𝟐. 𝟑𝟗)


̅ = 𝟏𝟎𝟎 ∗ (
𝝎 )
𝑺𝒘𝒆𝒃

Ce dernier ratio est défini pour des sections rectangulaires,𝐴𝑓𝑙 étant la section des aciers de
flexion concentré à chaque extrémité.
M: Moment fléchissant à la base,
T: Effort tranchant à la base (la contrainte de cisaillement est calculée en négligeant la
surface des ailes) :

𝑻 (𝟐. 𝟒𝟎)
𝝉=
𝑺𝒘𝒆𝒃

N : Effort normal,

45
h : Hauteur du mur,
h’ : Hauteur correspondant au point d'application des forces horizontales,
q : Coefficient de comportement.

-Efforts tranchants de fissuration :


Deux contraintes de cisaillement sont calculées, la première 𝜏1 correspondant à un
mécanisme de flexion et la seconde 𝜏2 à un mécanisme d'effort tranchant.

a. Mécanisme de flexion
𝑴𝒍𝒊𝒎 (𝟐. 𝟒𝟏)
𝝉𝟏 = 𝝉∗ ∗
𝑴

-La contrainte de cisaillement 𝜏 ∗ vaut :


𝝉∗ =(((1+q))/2)* T) / (𝑺𝒘𝒆𝒃 ) (2.42)

Et
𝑴 = 𝑻 ∗ 𝒉′ (𝟐. 𝟒𝟑)

1+𝑞
-On retrouve le coefficient pénalisant ( ) pour "annuler" son effet : ce facteur ne doit
2
s'appliquer qu'aux mécanismes de cisaillement et l'effort tranchant 𝜏1 correspond au
mécanisme de flexion.
-Le moment limite de fissuration vaut :
𝟐 ∗ 𝑰𝒏 𝒇𝒕 (𝟐. 𝟒𝟒)
𝑴𝒍𝒊𝒎 = ∗ ( + 𝝈𝒏 )
𝒃 𝒃

Avec :

Pour une section en (H) :

𝒂 ∗ 𝒃𝟑 𝑺𝒇𝒍 ∗ 𝒃𝟐 𝟔 ∗ 𝑺𝒇𝒍 𝒂 ∗ 𝒃𝟑 𝟔 ∗ 𝑺𝒇𝒍 𝑺𝒘𝒆𝒃 ∗ 𝒃𝟐 (𝟐. 𝟒𝟓)


𝑰𝒏 = + = (𝟏 + )∗ = (𝟏 + )∗
𝟏𝟐 𝟐 𝑺𝒘𝒆𝒃 𝟏𝟐 𝑺𝒘𝒆𝒃 𝟏𝟐

Pour une section rectangulaire :

𝒂 ∗ 𝒃𝟑 𝑺𝒘𝒆𝒃 ∗ 𝒃𝟐 (𝟐. 𝟒𝟔)


𝑰𝒏 = =
𝟏𝟐 𝟏𝟐

D’où :
𝟏+𝒒 𝟐𝑰𝒏 𝒇𝒕 (2.47)
𝝉𝟏 = ( ) ( + 𝝈𝒏 )
𝟐 𝒃𝑺𝒘𝒆𝒃𝒉′ ɣ𝒃

Avec :
Pour une section en (H) :
𝟐𝑰𝒏 𝑺𝒇𝒍 𝒃 (𝟐. 𝟒𝟖)
= (𝟏 + 𝟔 )
𝒃𝑺𝒘𝒆𝒃 𝑺𝒘𝒆𝒃 𝟔𝒉′

Pour une section rectangulaire :


𝟐𝑰𝒏 𝒃 (𝟐. 𝟒𝟗)
=
𝒃𝑺𝒘𝒆𝒃 𝟔𝒉′
.

46
b. Mécanisme d'effort tranchant
-La formule suivante est utilisée :
𝟐 (𝟐. 𝟓𝟎)
𝝉𝟐 =0.45√𝒇𝒕 (𝒇𝒕 + 𝟑 𝝈𝒏 )

Cette formule est valable pour un coefficient de sécurité du béton valant 1.5.
Dans le cas accidentel (ɣ𝑏 =1.15), il faut la modifier :
𝟏.𝟓 𝟏.𝟓 𝟐 (𝟐. 𝟓𝟏)
τ2= 0.45√𝟏.𝟏𝟓 𝒇𝒕 (𝟏.𝟏𝟓 𝒇𝒕 + 𝟑 𝝈𝒏 )

c. Résistance limite associée


-Une contrainte de cisaillement limite est associée aux deux mécanismes de fissuration :
𝝉𝟑 = (𝟏 + 𝟑𝝎 ̅̅̅) ∗ 𝒎𝒊𝒏(𝝉𝟏 , 𝝉𝟐 ) + 𝟎. 𝟏𝟓 ∗ 𝝈𝒏 (𝟐. 𝟓𝟐)

-Cette formule ne tient compte que des aciers de flexion concentrés aux extrémités du voile.
Les aciers répartis apportent une résistance supplémentaire valant:
𝜟𝝉 = 𝟎. 𝟗 𝒇𝒚 𝝆𝒘𝒆𝒃 (𝟐. 𝟓𝟑)

2.3.4.3. E.C.2

EC2 c’est un code européen pour le dimensionnement des structures en béton armé,
pour la partie évaluation de la résistance des éléments.
-En ce qui concerne la résistance à l’effort tranchant, l’Eurocode 2 propose des règles pour les
poutres ou consoles courtes qui peuvent être appliquées aux voiles en béton armé.

a. Résistance au cisaillement sans ferraillage réparti qui correspond à la contribution du


béton :
La formule ci-dessous qui prend en compte l’effet du ferraillage de flexion :
𝟐𝒅
𝑽𝑹𝒅,𝒄 = [ 𝑪𝒓𝒅 K(𝟏𝟎𝟎𝝆𝒍 𝒇𝒄𝒌 )𝟏/𝟑 (𝒂 ) + 0.15𝝈𝒄𝒑 ]𝒃𝒘 d ≤ 0.5𝒃𝒘 dν𝒇𝒄𝒅 (2.54)
𝒗

Avec :
𝒄𝒌 𝒇 (2.55)
ν = 0.6 [1-𝟐𝟓𝟎 ]

d : hauteur utile de la section en millimètre.


𝑎𝑣 : Bras de levier de l’effort tranchant, si 𝑎𝑣 < 0.5d, on prend

𝒂𝒗 = 0.5d. (2.56)

𝑏𝑤 : Largueur de section en millimètre.


H : Hauteur de la section en millimètre.
𝐴𝑠𝑙 : Section d’armature de flexion
𝑓𝑐𝑘 : Résistance du béton en compression à 28 jours.

𝟎. 𝟏𝟖 (2.57)
𝑪𝑹𝒅,𝒄 =
ɣ𝒄

𝟐𝟎𝟎 (2.58)
K=1+√ ≤ 2.0
𝒅

47
Taux de ferraillage pour la flexion :

𝑨
𝝆𝒍 = 𝒃 𝒔𝒍𝒅 ≤ 0.02 (2.59)
𝒘

Contrainte normale en MPa :

𝑵
𝝈𝒄𝒑 = 𝒃 𝑬𝒅𝒉 ≤ 0.2𝒇𝒄𝒅 (2.60).
𝒘

𝑁𝑒𝑑 : Effort normal en N, positif en compression.


Résistance de dimensionnement du béton en Compression en MPa.

𝒇𝒄𝒌 (𝟐. 𝟔𝟏)


𝒇𝒄𝒅 =
ɣ𝒄

Pour les états-limites ultimes en situation accidentelle : ɣ𝑠 = 1 et ɣ𝑐 = 1.2

b- La résistance à l’effort tranchant avec des armatures transversales réparties (non


inclinées) est la suivante :

𝑽𝑹𝒅 = 𝑽𝑹𝒅,𝒄 + 𝑨𝒔𝒘 𝒇𝒚𝒘𝒅 ≤ 𝟎. 𝟓𝒃𝒘 𝒅𝝂𝒇𝒄𝒅 (𝟐. 𝟔𝟐)

Avec :
𝑉𝑅𝑑,𝑐 : La contribution du béton calculée plus haut.
𝐴𝑠𝑤 : La somme des sections des armatures réparties sur largeur de 0.75𝑎𝑣 .
𝑓𝑦𝑤𝑑 : Limite élastique des armatures réparties.

0.5*𝑏𝑤 *d*ν*𝑓𝑐𝑑 représente la limite supérieure de la résistance à l’effort tranchant et


correspond probablement à l’écrasement du béton de l’âme.
-Pour la vérification du non-glissement horizontal des voiles, on peut utiliser la formule
concernant la résistance au cisaillement des reprises de bétonnage :
𝑽𝑹𝒅𝒊 = 𝒄𝒇𝒄𝒕𝒅 + 𝝁𝝈𝒄𝒑 + 𝝆𝒇𝒚𝒅 (𝟐. 𝟔𝟑)

Où :
𝑉𝑅𝑑𝑖 : Résistance au glissement en MPa.
𝑓𝑐𝑡𝑑 : Résistance en traction du béton.
ρ : taux de ferraillage traversant l’interface.
𝑓𝑦𝑑 : Limite élastique nominale des armatures.
c et μ : facteurs dépendant de la rugosité de l’interface, par exemple c=0.45 et μ = 0.7 pour
une interface rugueuse.

2.3.4.4. E.C.8

L’Eurocode 8 est applicable sur des constructions à risque normal. C’est une norme
qui se base comme les PS 92 sur le principe de dimensionnement en capacité.
-EC8 pour la partie évaluation des sollicitations, la détermination de la résistance et pour les
dispositions constructives supplémentaires.
-L’Eurocode 8 distingue 3 niveaux de ductilité (basse, moyenne et haute) :

48
a. le niveau le plus bas : la résistance au cisaillement est strictement calculée selon
l’Eurocode 2 (un coefficient de ductilité de 1,5 est autorité sans aucune disposition
constructive particulière dans les zones de faible sismicité).

b. les voiles à ductilité moyenne : on distingue deux classes, celle faiblement armé
fonctionnant en flexion + cisaillement avec élévation du centre de gravité et les autres dits «
ductiles » qui sont plus fortement armés. Dans les deux cas, les longueurs d’ancrage dans les
zones potentielles de glissement sont majorées de 50%.

c. Le niveau de ductilité élevée : (au-delà de la notion de coefficient de ductilité) introduit


des conditions de résistance différentes de celles données par l’Eurocode 2.

2.3.4.5. ACI 349-01

L’ACI-349-01 (American Concrete Institute, version 2001) est un code


spécifiquement dédié aux structures en béton armé des installations nucléaires. Ce texte
précise notamment les combinaisons de chargements accompagnant le séisme dont la valeur
la plus élevée.
Concernant la résistance à l’effort tranchant des voiles en béton armé, il utilise le coefficient
de réduction de la capacité des éléments (Φ) (varie entre 0.85 et 0.9).

2.4.3.6. JEAG40611987

Le règlement japonais qui a servi pour le dimensionnement des centrales nucléaires.


Pour les voiles en béton armé, ce code permet de calculer la résistance d’effort tranchante par
la formule suivante :
𝟎.𝟎𝟔𝟕𝟗𝑷𝟎.𝟐𝟑
𝒗 ∗(𝑭𝒄 +𝟏𝟖𝟎) (2 .64)
𝝉𝒖 = 𝑴
+ 𝟐. 𝟕√𝑷𝒉 𝝈𝒘𝒚 + 𝟎. 𝟏𝝈𝒗 .
√𝑸𝑫+𝟎.𝟏𝟐

Avec :
M : Moment de flexion.
Q : Effort tranchante.
D : Largeur de voile.
𝑃ℎ : Pourcentage d’armature horizontale.
𝑃𝑣 : Pourcentage d’armature verticale.
𝐹𝑐 : Résistance caractéristique de béton en bars.
𝜎𝑤𝑦 : Limite élastique des armatures transversales.
𝜎𝑣 : Contrainte normale de compression en bars.

De nombreux articles ont été publiés expliquant le calcul effort tranchant des voiles en béton
armé. Parmi ces publications, l’article publié par M.J.N. Priestley, F.Seible et G.M. Calv cette
publication prend en compte la contribution de béton et la contribution de l’acier et la force de
frottement sous effort normale :
𝑽𝒏 = 𝑽𝑪 +𝑽𝑺 + 𝑽𝑷 . (2.65)

𝑽𝒄 = K√𝒇𝒄 𝑨𝒆 . (2.66)

𝑨𝒗 𝒇𝒚 𝑫′𝒄𝒐𝒕Ɵ (2.67)
𝑽𝑺 = .
𝑺

49
𝑽𝑷 = p tgα. (2.69)

Avec :
K= 0.29.
𝐴𝑒 : Section d’effort tranchant= 0.8 A.
𝐷′ : Bras de levier.
S : Espacement des armatures transversales.
α : Angle de frottement.
Ɵ : Inclinaison de bielle en général 30°.

2.4. Conclusion

D’après l’analyse des différents codes de calcul des ouvertures des fissures, on a constaté que
la plupart des codes estiment les ouvertures au voisinage des contraintes des aciers en
négligeant la contribution de béton. Les ouvertures des fissures dépendent de plusieurs
paramètres tels que le taux armature et la durée de chargement.
Dans les codes de dimensionnement, l'effort tranchant V repris par une section est divisé en
un terme apporté par le béton 𝑉𝐶 («concrete contribution») et un terme apporté par l'acier
transversal 𝑉𝑆 («steel contribution»).
La majorité des codes de calcul classifient les degrés de dommage en fonction des ouvertures
des fissures et le coefficient de déplacement entre les étages.

50
Chapitre 03

Validation des approches numérique sur des voiles courts soumise aux
cisaillements.

51
3.1. Introduction
Les voiles sont des structures courant dans les constructions ordinaires et spéciales (les
installations nucléaires par exemple) et sont généralement utilisés pour reprendre les efforts
de cisaillement (dans le cas sismique ou statique).
Par ailleurs, les voiles en cisaillement et leur comportement vis-à-vis des diverses
sollicitations deviennent un axe de recherche très important. Parmi les programmes qui ont
été lancés ces dernières années pour prédire le comportement non linéaire des voiles le
programme (SAFE Structure Armée Faiblement Elancé ,1998).
Le but de ce chapitre est d’évaluer le comportement non linéaire des voiles vis-à-vis du
cisaillement : La réponse statique monotone (push over) et cyclique-statique est évaluée en
utilisant différentes lois de comportement pour le béton et l’acier. L’approche semi-global
(multifibre et multi couche) et l’approche macroscopique sont utilisées. Les simulations sont
faites à l’aide du logiciel élément finis CAST3M.
L’objectif est de voir si les deux approches sont capables de représenter la réponse globale et
locale des structures en béton armé, les résultats de simulation seront confrontés aux résultats
des expérimentations qui sont ont été effectuées dans le cadre du programme SAFE.

3.2. Présentation du programme SAFE

Les voiles fortement armés faiblement élancés constituent une partie importante dans
les ouvrages spéciaux tels que les ouvrages nucléaires d’où l’importance de suivi de leur
sécurité l’évalué de leur vulnérabilité. Ceci a conduit les entreprises responsables a réaliser
des programmes soit expérimental ou numérique pour estimer le degré des dommages des
structures soumis à différentes sollicitations tel que les efforts de cisaillement dynamique
(séisme par exemple).
Le programme SAFE (Structure Armé Faiblement Elancé) a été créé par un groupe
d’exploitant nucléaire (COGEMA et EDF) en 1998 avec aide de la communauté européenne
sur une série des voiles en béton armé. L’objectif principal de l’action de recherche autour de
ces essais c’est d’améliorer les codes de calcul sismique des ouvrages nucléaires.
Les essais présentés ci-dessous ont les particularités suivantes (Pegon et al. 1998a) :
1. De tester une petite zone de voile soumise uniquement à du cisaillement, et non un voile
complet en flexion ;
2. De se placer au plus près des conditions de réalisation effectives (dispositions constructives
et taux de ferraillage correspondant à l'évolution en cours dans le nucléaire) :
3. D'utiliser une excitation sismique ;
4. D'étudier la variation de fréquence propre de la structure avec l'endommagement et
l'influence du rapport fréquence propre sur fréquence d'excitation sur les effets dynamiques ;
5. De mieux décrire et quantifier la phase de début d'endommagement correspondant à
l'apparition des premières fissures et au début de la plastification des armatures. [Panagiotis
K. ,2000].

3.2.1. Caractéristiques de la maquette

Les murs, ses raidisseurs ainsi que les longrines supérieure et inférieure assurant le
chargement sont représentés sur la figure 3.1. Les principales caractéristiques des 13
maquettes du programme sont regroupées dans le tableau 3.2

52
figure3. 1: plan de ferraillage (T10).
Maquette a b h s M ρv ρh σv τd K ξ f

m m m m² (tonne) (%) (%) (MPa) (MPa) (MN/m) (%) HZ


T1 3.0 0.16 1.2 0.48 480 0.8 0.8 4 0 5696 7 17.3

T2 3.0 0.16 1.2 0.48 5700 0.8 0.8 4 0 5696 7 5.03

T3 3.0 0.16 1.2 0.48 9018 0.8 0.8 4 0 7120 7 4

T4 3.0 0.16 1.2 0.48 1102 0.8 0.8 4 0 7120 7 12

T5 3.0 0.4 1.2 0.6 2818 0.8 0.8 4 0 7120 7 8

T6 3.0 0.4 1.2 0.6 11252 0.6 0.4 3 1 7120 7 12

T7 3.0 0.4 1.2 0.6 11272 0.6 0.4 3 1 7120 7 4

T8 3.0 0.4 1.2 0.6 11252 0.4 0.4 2 0 7120 7 12

T9 3.0 0.4 1.2 0.6 11272 0.4 0.4 2 0 7120 7 4

T10 3.0 0.4 1.2 0.6 11272 0.6 0.6 3 0 7120 7 4

T11 3.0 0.4 1.2 0.6 11272 0.4 0.4 2 0 7120 7 4

T12 3.0 0.4 1.2 0.6 11272 0.11 0.11 1.44 1 7120 7 4

Tableau3. 1: description de la maquette SAFE.

53
𝒇𝒄𝟐𝟖 = 𝟒𝟎𝑴𝑷𝒂 , 𝒇𝒚 = 𝟓𝟎𝟎𝑴𝑷𝒂.
𝒇𝒄𝟐𝟖 , 𝒇𝒚 : La résistance caractéristique de béton, nuance d’acier.
a, b, h : La géométrie des voiles.
𝝆𝒗 et 𝝆𝒉 : taux d’armature suivant les deux directions.
𝝈𝒗 ,𝝉𝒅 et K contrainte de compression, contrainte de cisaillement rigidité conventionnelle.
ξ, f : amortissement, fréquence excitation.

figure3. 2:plan de ferraillage (T10).

figure3. 3: plan de ferraillage (T12).

3.2.2. Condition aux limites

Tous les murs sont encastrés à la base. Des masses additionnelles posées en tête de
voile sont à l'origine de la contrainte verticale moyenne à la base. Un dispositif de chargement
a été conçu et construit de telle façon que la rotation de la longrine supérieure soit empêchée.
Ce dispositif est réutilisé pour chacun des tests (Figure 3.4). On remarque le contrefort
additionnel du côté Est, qui permet de solliciter le dispositif de chargement des deux côtés à la
fois. On remarque également que le dispositif de blocage à la base, qui empêche le spécimen
de glisser. Des contrevents latéraux ont été introduits, de façon à empêcher les éventuels
mouvements hors plan.

54
figure3. 4: illustration de principe du dispositif de chargement.

3.2.3. Essai de pseudo dynamique

La méthode pseudo dynamique (PSD) est une méthode hybride


numérique/expérimentale qui combine le calcul du déplacement de la structure et la mesure
de la force employée pour imposer ce déplacement. Pour des applications dans le domaine du
génie parasismique, un accélérograme naturel ou généré artificiellement est introduit dans
l'ordinateur équipé d'un algorithme PSD. Pour des petits pas de temps on calcul les
déplacements correspondants aux degrés de liberté des points où la masse de la structure est
supposée concentrée. Ces déplacements sont ensuite imposés à la structure par l'intermédiaire
de vérins hydrauliques fixés sur le mur de réaction. Les réactions sont mesurées et introduites
dans l'ordinateur pour le pas suivant. Les forces d'inertie (et éventuellement d'amortissement
visqueux) sont calculées numériquement, ce qui permet de réaliser des tests de structures de
taille réelle avec une échelle de temps dilatée et des besoins hydrauliques modérés (Donea et
al. 1995). Il devient ainsi très simple de varier la fréquence propre du système en changeant la
valeur de la masse fictive en translation M2 (Figure 3.2). La méthode PSD constitue une
approche complémentaire aux essais dynamiques sur table vibrante où seulement des
maquettes à échelles réduites de structures peuvent être testées.

3.3. Modélisation numérique

Les simulations numériques sont effectuées sur deux voiles proposés dans le cadre de
programme SAFE (T10 et T12) qui sont sollicités par cisaillement statique et dynamique par
deux approches de modélisation (multicouche et multifibre) et macroscopique afin de trouver
la réponse globale et locale des voiles et aussi faire une comparaison entre les deux approches
en terme ouverture des fissures et la réglementation.
Les deux voiles sont en béton armé et sollicité en cisaillement statique et un cisaillement
statique-dynamique, les deux spécimens ont une portée de 3m, une hauteur de 1.20m et une
largeur de 20 cm avec un taux de ferraillage différent.
Les caractéristiques géométrique et mécanique sont représentées dans le tableau (tableau 3.2)

Maquette T10 T12


Type de test Essai PSD Essai PSD
Condition aux limites Rotation en tête bloquée Rotation en tête

55
bloquée
Élancement 0.46 0.46
Section du voile m² 0.52 0.52
Section de nervure m² 0 .128 0.128
Renforc.Horiz. % 0.6 0.11
Renforc .Vert. % 0.6 0 .11
Espacement mm 125 250
Résistance en compression MPa 40 40
Résistance en traction MPa 3. 3.
Module de Young béton MPa 30000 30000
Limite élasticité MPa 570 570
Module de Young acier MPa 200000 200000
Cont.normale en base MPa 0.34 1
Mass (long .sup+extra Kg 25000 25000
mass)

Tableau3. 2:caractéristique géométrique et mécanique des maquettes SAFE.

La modélisation est effectuée sur des éléments en multifibres avec une poutre de type
Timochenko à deux nœuds, et des éléments en multicouches avec coque mince. Quant à la
modélisation macroscopique, tous les éléments sont présentés (Voir Figure 3.5)
Les voiles en multi fibre sont divisés sur la hauteur en 8 sections, chaque section comporte 20
fibres en béton comportant les aciers longitudinaux pour chaque section suivant toute la
longueur de voile. On suppose qu’il y’ a une adhérence parfaite entre l’acier et le béton.
Les voiles en multicouche sont divisé par trois couches de béton et 2 nappes de ferraillage
longitudinal avec des armatures transversales.
Dans la modélisation macroscopique, des éléments finis tridimensionnels CUB8 sont utilisés
pour le béton, et des éléments poutres pour les aciers.
Les lois de comportement sont basées sur le mécanique endommagement de béton et de la
plasticité parfait avec écrouissage pour les aciers.

(a)
(b)

56
(c)

figure3. 5: différentes approches utilisées dans la modélisation (a)multifibre (b)


multicouche (c) approche macroscopique.

3.4. Les modèles utilisés pour la simulation


3.4.1. Le modèle de Laborderie

Pour la modélisation multifibre, afin de modéliser le comportement non linéaire du


béton on a choisi, le modèle unilatéral de la Borderie qui établit par la borderie 1991 qui
prend en considération les effets unilatéraux et est caractérisé par :
- Deux variables qui contrôlent l’endommagement (une variable pour la traction et une
variable pour la compression).
- Déformations anélastiques couplées avec l’endommagement.
- Phénomènes unilatéraux. Ouverture et refermeture des fissures gérées par une fonction
de rigidité.

-
- figure3. 6:réponse uni axiale du modèle [la borderie, 1991]. [Panagiotis.K, 2006].

57
Les paramètres utilisés pour la loi unilatérale sont :

Les paramètres de modèle Valeur


YS1 (Seuil d’endommagement en traction) 3*10²
YS2 (Seuil d’endommagement en compression) 5*104
A1 (Paramètre d’endommagement en traction) 5*10−3
A2 (Paramètre d’endommagement en compression) 6*10−6
B1 (Paramètre pour la traction) 1.4
B2 (Paramètre pour la compression) 1.5
BETA1 (Activation de la déformation permanente en traction) 1*106
BETA2 (Activation de la déformation permanente en -40*106
compression)
SIGF (Contrainte de refermeture de fissure) 3.5*106
Youn (module de Young) 30000
NU (coefficient de poisson) 0.2

Tableau3. 3: Tableau des paramètres du modèle Unilatéral utilisés dans la simulation.


[D. COMBESCURE, 2001].
3.4.2. Le modèle multicouche

Paramètres de modèle valeur

YOUN² (module de Young) 3*10² (MPa)

NU (coefficient de poisson) 5*104

FT (Résistance en traction) 3*106 (MPa)

FC (Résistance en compression) 30*106 (MPa)

𝐺𝑓𝑡 (énergie de fissuration en traction) 200 J.m

𝐺𝑓𝑐 (énergie de fissuration en compression) 2000J.m

Tableau3. 4:les paramètres de modèle multicouche.

58
3.4.3. Le modèle fichant (micro-iso)

Pour la modélisation macroscopique on a choisi le modèle de fichant et al (1996)


basée sur la théorie de plasticité et de l’endommagement avec prise en compte la refermeture
des fissures, ce modèle permet de présenter correctement l’effet de cisaillement pour les
structures en béton armé. Le tableau (3.5) donne les paramètres de modèle fichant :

Paramètre de modèle Valeur


'EPSDO' Seuil en déformation pour la traction. 1.E-4
'FC' résistance en compression simple. 42Mpa
‘𝐺𝑓 ′ énergie du rupture. 100J/m²
‘ALFA’ facteur influence. 40
‘MP1’ paramètre de plasticité. 10*𝑌𝑔

Tableau3. 5: les paramètres de fichant utilisé dans la modélisation

3.4.4. Loi de comportement de l’acier


Dans notre modélisation on a choisi de modéliser les aciers par le modèle élasto-
plastique avec écrouissage linéaire, la figure donne des illustrations sur le modèle parfait-uni.

figure3. 7:loi de comportement de l’acier élasto-plastique. [D. COMBESCURE,


2001].

3.3.5 Les paramètres utilisés pour la loi élasto-plastique sont :

Paramètre de modèle Valeur


Young (Module de Young) 2.1E+11
SIGY (Limite d’élasticité) 570E+6
H (Module d’écrouissage linéaire) 10E+8

Tableau3. 6:Tableau des paramètres du modèle élasto-plastique utilisés dans la


simulation.

59
3.5. Histoire de chargement

Dans notre modélisation on a appliqué deux types de chargement sur les deux voiles
T10, et T12.Le premier chargement c’est un chargement statique monotone en déplacement
imposé (push-over) appliqué en tête des voiles. Le deuxième chargement c’est un chargement
statique-cyclique proposé par le programme expérimental SAFE et le point d’excitation c’est
identique comme le déplacement imposé. Les 2 figures présentent le chargement imposé
appliqué sur les voiles T10 et T12.

2,00E+01

1,50E+01

1,00E+01

5,00E+00
depl (mm)

0,00E+00
0,00E+00 1,00E+03 2,00E+03 3,00E+03 4,00E+03 5,00E+03 6,00E+03
-5,00E+00

-1,00E+01

-1,50E+01

-2,00E+01
temps (s)

figure3. 8:: déplacement –temps (T10).

T12-
5
T12-4

T12-1 T12-2 T12-3

figure3. 9: : déplacement –temps (T12).

60
3.6. Résultats et discussions

Dans ce qui suit, on présente la réponse globale (effort-déplacement) obtenus avec les
résultats expérimentaux. Dans un premier temps, les résultats des modélisations simplifiés
(Multicouche/Multifibre) sont comparés avec l’expérimentation.

3.6.1. Résultat push-over (Statique monotone)


3.6.1.1. Voile T10

5,00E+02
Effort(KN)

4,00E+02
3,00E+02
2,00E+02
1,00E+02
0,00E+00
-2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-1,00E+02
-2,00E+02
-3,00E+02
-4,00E+02 expérimental
PUSH OVER t10
-5,00E+02
Déplacement (mm) a

8,00E+02

6,00E+02
Effort(KN)

4,00E+02

2,00E+02

0,00E+00
-2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01 2,50E+01
-2,00E+02

-4,00E+02
expérimental
pushover T10
-6,00E+02
Déplacement(mm) b

Figure3.10 la réponse global (expérimental avec numérique (statique)) pour voile


T10 (a) Approche Multifibre (b) Approche Multicouche.

61
3.6.1.2. Voile T12

5,00E+02
Effort (KN)

4,00E+02

3,00E+02

2,00E+02

1,00E+02

0,00E+00
-2,50E+01 -2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-1,00E+02

-2,00E+02

-3,00E+02 expérimental
PUS OVER T12
-4,00E+02

Déplacement(mm) a
Effort(KN)

8,00E+02

6,00E+02

4,00E+02

2,00E+02

0,00E+00
-2,50E+01 -2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-2,00E+02
expérimental
push-over T12
-4,00E+02
Déplacement (mm) b

Figure3.11 la réponse global (expérimental vs numérique (statique)) pour voile


T12 (a) Approche Multifibre (b) Approche Multicouche.

3.6.2. Résultat cyclique


Les mêmes comparaisons sont réalisés pour le cas de chargement statique-Cyclique :

3.6.2.1. Approches multi couches

62
6,00E+02

4,00E+02

2,00E+02
Effort(KN)

0,00E+00
-2,00E+01-1,50E+01-1,00E+01-5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-2,00E+02

-4,00E+02
expérimental
-6,00E+02 numérique

-8,00E+02 a
Déplacement (mm)

6,00E+02

4,00E+02
Effort
(KN)

2,00E+02

0,00E+00
-2,50E+01-2,00E+01-1,50E+01-1,00E+01-5,00E+000,00E+005,00E+001,00E+011,50E+012,00E+01
-2,00E+02

-4,00E+02
cyc exp T12
-6,00E+02 cyc num T12

-8,00E+02
Déplacement b
(mm)

figure3. 10:la réponse global (expérimental vs numérique (cyclique)) pour voile (a)
T10 et (b) T12-multi couches.
3.6.2.2 Approche multi fibres :

63
5,00E+02
4,00E+02
3,00E+02
Effort (KN)

2,00E+02
1,00E+02
0,00E+00
-2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-1,00E+02
-2,00E+02
-3,00E+02 exp-cyc
-4,00E+02 cyc-num

-5,00E+02

Déplacement(mm) a
4,00E+02

3,00E+02
Effort (KN)

2,00E+02

1,00E+02

0,00E+00
-2,50E+01 -2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-1,00E+02

-2,00E+02

exp-cyc
-3,00E+02
num-cyc
-4,00E+02

Déplacement (mm) b

figure3. 11:la réponse global (expérimental avec numérique (cyclique)) pour voile T10
(a) et T12 (b)-multi fibres.

64
3.6.3 Comparaison approche macroscopique Vs. Expérimentation
3.6.3.1 Courbe effort-déplacement (statique-monotone)
8,00E+03
Effort (KN)

6,00E+03

4,00E+03

2,00E+03

0,00E+00
-2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-2,00E+03

-4,00E+03
exp
-6,00E+03 macro
semiglobal
-8,00E+03

Déplacement (mm) a
5,00E+03
Effort(KN)

4,00E+03

3,00E+03

2,00E+03

1,00E+03

0,00E+00
-2,50E+01 -2,00E+01 -1,50E+01 -1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01 1,50E+01 2,00E+01
-1,00E+03

-2,00E+03

-3,00E+03
exp
-4,00E+03 macro
semi-global
-5,00E+03
Déplacement (mm) b
figure3. 12:courbe de capacité macro avec semi-global (multicouche).
Un chargement statique constitue une sorte de courbe enveloppe à la réponse non linéaire
statique cyclique (figure 3.14). L’intérêt principal de ce calcul est d’analyser le comportement
de la structure pour des niveaux de non linéarité (endommagement) proches de ceux
rencontrés dans le cadre du calcul dynamique non linéaire.
La confrontation de la réponse globale expérimentale avec les résultats numériques montre
que l’effort ultime estimé par le modèle passe par la réponse globale expérimentale. On peut
dire aussi qu’après la plastification des armatures longitudinale la forme des cycles évolue
suivant leurs comportements (comportement des aciers).

3.6.3.2. Courbe effort déplacement (cyclique (T10))

65
2,00E+03

Effort(KN)
1,50E+03

1,00E+03

5,00E+02

0,00E+00
-1,00E+00 -5,00E-01 0,00E+00 5,00E-01 1,00E+00
-5,00E+02

-1,00E+03

-1,50E+03
cyc exp 1
-2,00E+03
Déplacement(mm)
2000
Effort(KN)

1500

1000

500

0
-10 -5 0 5 10
-500

-1000

-1500
numcyc1
-2000
Déplacement(mm)
5,00E+03
Effort(KN)

4,00E+03
3,00E+03
2,00E+03
1,00E+03
0,00E+00
-1,00E+01 -5,00E+00 0,00E+00 5,00E+00 1,00E+01
-1,00E+03
-2,00E+03
-3,00E+03
cyc exp 2
-4,00E+03
-5,00E+03
Déplacement(mm)

66
1500

Effort(KN)
1000

500

0
-10 -5 0 5 10

-500

-1000
numcyc2

-1500
Déplacement(mm)
8,00E+03
Effort(KN)

6,00E+03

4,00E+03

2,00E+03

0,00E+00
-1,50E+01-1,00E+01-5,00E+000,00E+005,00E+001,00E+011,50E+01
-2,00E+03

-4,00E+03 cyc exp 3

-6,00E+03
Déplacement(mm)
-8,00E+03

67
3,00E+03

Effort (KN)
2,00E+03

1,00E+03

0,00E+00
-1,50E+01-1,00E+01-5,00E+000,00E+005,00E+001,00E+011,50E+01

-1,00E+03

-2,00E+03
cyc num 3
-3,00E+03

Déplacement (mm)

figure3. 13:comparaison différent cycles numériques avec expérimentales SAFE (T10)


par approche macroscopique.
3.5.3.3 Courbe effort déplacement (cyclique (T12)) :

68
1,50E+03

Effort (KN)
1,00E+03

5,00E+02

0,00E+00
-3,00E-01-2,00E-01-1,00E-010,00E+00 1,00E-01 2,00E-01 3,00E-01
-5,00E+02
cyc exp 1
-1,00E+03

-1,50E+03
Déplacement (mm)
1500
Effort (KN)

1000

500

0
-0,2 -0,1 0 0,1 0,2 0,3

-500
cycnum1

-1000

Déplacement (mm)

69
2,00E+03

Effort (KN)
1,50E+03

1,00E+03

5,00E+02

0,00E+00
-1,00E+00 -5,00E-01 0,00E+00 5,00E-01 1,00E+00
-5,00E+02

-1,00E+03

-1,50E+03
cyc exp 2
-2,00E+03
Déplacement (mm)

1500
Effort (KN)

1000

500

0
-1,00E+00 -5,00E-01 0,00E+00 5,00E-01 1,00E+00

-500

-1000

numcyc2
-1500
Déplacemen(mm)

70
3,00E+03

Effort (KN)
2,00E+03

1,00E+03

0,00E+00
-4,00E+00 -2,00E+00 0,00E+00 2,00E+00 4,00E+00
-1,00E+03

-2,00E+03
cyc exp 3
-3,00E+03
Déplacement (mm)

600
Effort (KN)

400

200

0
-3,00E+00 -2,00E+00 -1,00E+00 0,00E+00 1,00E+00 2,00E+00 3,00E+00

-200

-400
numcyc3
-600
Déplacement (mm)

figure3. 14: comparaison différent cycles numériques avec expérimentales SAFE


(T12) par approche macroscopique.

3.6.3.4. Discussions des résultats globaux :

La superposition des résultats expérimentaux et calculées montre que les efforts


ultimes ont été bien estimés, mais les résultats numériques montrent aussi que les cycles 3 et
4 dissipent légèrement, une fois il y a une plastification des aciers on remarque qu’il y a une
dissipation d’énergie.
On peut noter que la diminution de la raideur globale augmente avec l’amplitude des cycles,
signe d’une macro-fissuration ainsi que l’augmentation de l’énergie dissipée au sein de
chacun des cycles simultanément à l’accroissement de la fissuration.

71
3.7. Résultats en termes comportement local (endommagement-fissuration)
3.7.1. Cas d’un chargement push-over

Avec le modèle d’endommagement du béton, la dégradation mécanique peut être


représentée par l’évolution de l’endommagement dans la structure. La figure (3,17) et figure
(3.18) donnent des illustrations sur le processus de l’endommagement de traction dans le
béton pour différents déplacements pour les voiles T10 et T12.

a. Voile T10 :

72
figure3. 15:évolution endommagement par traction (cisaillement) voile T10.

b. Voile T12

73
figure3. 16: évolution endommagement par traction voile T12.

3.7.2. Cas cyclique


Après la validation des différentes approches de modélisation avec l’expérimental en
terme global (effort –déplacement), on va présenter l’évolution de l’endommagement pour les
différents cycles et faire une comparaison avec l’expérimental.
a. Voile T10

74
75
76
figure3. 17:évolution d’endommagement pour différent cycle T10 et comparaison
avec l’expérimental.

77
b. Voile T12

78
figure3. 18: évolution d’endommagement pour différent cycle T12 et comparaison
avec l’expérimental.
3.7.3. Discussion des résultats locaux (cartes d’endommagement)

Les figures (3.19) et (3.20) montrent l’évolution de l’endommagement comparée avec


l’ordre de répartition des ouvertures des fissures de l’expérience. Les réponses sont
qualitativement similaires entre la simulation et l’expérience. On peut dire que le modèle est
capable à prédire les zones d’endommagement critiques.

79
3.8. Conclusion

Les approches utilisées dans ce chapitre pour traiter le comportement non linéaire des
voiles courtes simulent correctement le comportement global des voiles jusqu'à la ruine, mais
pour les approches semi-globale son n’a pas pu reproduire l’endommagement des voiles de
maniéré fiable.
On a utilisé l’approche macroscopique pour définir la réponse globale et locale en terme
endommagement, afin d’améliorer les résultats expérimentaux, les essais sur les voiles ont été
effectués par des instruments de mesure de champs très développés, afin de déterminer les
ouvertures de fissures le long du chargement. L’objectif de chapitre suivant est d’évaluer les
ouvertures de fissures et de les comparer avec les résultats mesurés et celles prédites par la
réglementation.

80
Chapitre 04

Evaluation des ouvertures des fissures et confrontation avec réglementation


internationale.

81
4.1. Introduction

Les codes de calcul et de dimensionnement des structures en béton armé proposent


des formules empiriques pour l’évaluation des ouvertures des fissures dans les états limite de
service (état limitation des ouvertures des fissures) qui prennent en considération plusieurs
paramètres tel que la durée de chargement les classe exposition environnementale ………ex,
ces formules empiriques ne peuvent pas être généralisées pour toutes les structures en béton
armée.
L’objectif de ce chapitre est d’évaluer les ouvertures des fissures et leur propagation par trois
approches, la modélisation semi-globale et la modélisation macroscopique en utilisant la
procédure OUVFISS développée par Matallah et al (2010) et on compare les résultats
numériques avec l’expérimentation (programme SAFE). L’effort de cisaillement est calculé et
les résultats sont comparés avec les différents règlements internationaux.

4.2. L’évolution des ouvertures de fissures OUVFISS

Une approche de modélisation continue basée sur la mécanique d’endommagement est


employée pour calculer les ouvertures de fissure dans une structure en béton. M. Matallah et
al (2010) ont proposés une méthode de post traitement basée sur la régularisation énergétique
pour extraire les ouvertures de fissure à partir d'un calcul d’endommagement continu en
éléments finis. Cette méthode peut être appliquée à tous les modèles continus
d’endommagement ou/et plasticité [Medjahed Amina, 2012].
Dans les modèles basés sur les approches continues, on suppose que la fissure se localise dans
une bande de largeur h sur laquelle, on considère que les microfissures sont uniformément
distribuées.
La dissipation de l’énergie est constante, et donnée par la formule suivante :

(4.1)
𝑮𝒇 = ∫ 𝝈 𝒅𝜹
𝟎

Avec d𝛿 : saut déplacement.


Le saut de déplacement au niveau de la fissure est calculé comme étant le produit de la
déformation de rupture 𝜀 𝑓 et la largeur de la bande h.
Dans la zone de rupture, la déformation totale est définie comme la somme de la déformation
élastique et une déformation inélastique.
On considère que dans chaque élément il y a une fissure, la déformation de rupture totale est
sous la forme suivante :
𝒐𝒖𝒇 (4.2)
𝜺𝒊𝒋 = 𝜺𝒆𝒊𝒋 + 𝜺𝒊𝒋

Lors d’un calcul par élément fini, on peut estimer les déformations des ouvertures de fissures
à partir des contraintes anélastiques.
𝝈𝒊𝒏
𝒊𝒋 = 𝝈𝒊𝒋 − 𝝈
̃𝒊𝒋 (4.3)

Les contraintes effectives sont calculées à partir des déformations élastiques :


𝝈̃𝒊𝒋 = 𝑪𝒊𝒋𝒌𝒍 𝜺𝒌𝒍. (4.4)

La déformation d’ouverture de fissure est donnée par :


𝒐𝒖𝒇
𝜺𝒊𝒋 = 𝑪−𝟏 𝒊𝒏
𝒊𝒋𝒌𝒍 𝝈𝒊𝒋 .
(4.5)

L’ouverture normale à la direction de la fissure 𝜎𝑛 est donnée par :


𝝈𝒏 = 𝒏𝒊 𝝈𝒊𝒋 𝒏𝒋 (4.6)

82
L’ouverture des fissures est calculée en considérant que chaque élément est traversé par une
seule fissure, la valeur moyenne du déplacement est :
𝒐𝒖𝒇 𝒐𝒖𝒇 (4.7)
𝜹𝒏 = ∫ 𝜺𝒏 = ∫ 𝒏𝒊 𝜺𝒊𝒋 𝒏𝒋 𝒅𝒏.

Les paramètres du modèle sont donnés dans le tableau suivant :

Paramètre de modèle Valeur


3e+6 Pa
Kconst (Seuil d’endommagement en traction)
BETA (Seuil d’endommagement en compression) -3e+6 Pa
SSI (Constante du matériau) 4e+4
𝐺𝑓 (Energie de rupture) 100N/m
EPSR (La déformation à la rupture) 0.267
SIGF (Contrainte de la refermeture de fissure) 1.5e+6Pa

Tableau4. 1: identification paramètre de modèle (OUVFISS) utilisée dans la


simulation.

4.3. Les résultats du comportement local (Ouverture des fissures)


4.3.1. Cas statique

Les quatre figures donnent des illustrations sur l’évolution des ouvertures des fissures
par l’approche macroscopique, cette approche permet de représenter correctement la position
des ouvertures des fissures contrairement autres approches (multi fibre et multicouche).

83
a. T10

84
Figure4. 1:Illustration sur évolution des ouvertures des fissures pour voile T10.

85
b. T12 :

86
Figure4. 2:Illustration sur évolution des ouvertures des fissures pour voile T12.

4.3.2. Cas cyclique


a. T10

Figure4. 3:La carte de l’ouverture de fissures pour un déplacement de 0.78mm.

87
Figure4. 4:La carte de l’ouverture de fissures pour un déplacement de 5.43mm.

88
Figure4. 5: La carte de l’ouverture de fissures pour un déplacement de 11.57mm.

Figure4. 6:6 La carte de l’ouverture de fissures pour un déplacement de 13.50mm.

b. T12 :

Figure4. 7:La carte de l’ouverture de fissures pour un déplacement de 0.25 mm.

89
Figure4. 8:La carte de l’ouverture de fissures pour un déplacement de 0.85mm.

Figure4. 9: La carte de l’ouverture de fissures pour un déplacement de 2.60mm.

90
Figure4. 10:La carte de l’ouverture de fissures pour un déplacement de 7.8mm.

4.4. Comparaison des ouvertures des fissures avec la réglementation

Dans cette partie on expose les différentes formules de calcul des ouvertures de
fissures données par différents règlements. Les résultats des ouvertures obtenus
numériquement sont comparés à ceux évalués par la réglementation pour les deux voiles T10
et T12 (voir tableau 4.2).

 Norme SIA262

Les données :

Φ=10mm,𝐸𝑆 =210000MPa,𝑓𝑐𝑡 = 0.6% ∗ 350 = 2.1𝑁/mmet𝜎𝑠=350𝑀𝑃𝑎 .


10 3502
𝑤= ∗ 2.1∗210000 =0.347mm.
8

𝒘 = 𝟎. 𝟑𝟒𝟕𝒎𝒎.

 Fib model concrete:

𝜔 = 2 𝑙𝑠𝑚𝑎𝑥 *(𝜀𝑠 − 𝜀𝑐 )
𝑁 ′ 𝑥𝑥 𝑀𝑥𝑥
𝑁𝑥𝑥 = + 1*6* = -1074199.95 N/m.
𝑒 𝑒
𝑁 ′ 𝑦𝑦 𝑀𝑦𝑦
𝑁𝑦𝑦 = + 1*6* = 637512.5 N/m.
𝑒 𝑒
𝑁′ 𝑥𝑦 𝑀𝑥𝑦
𝑁𝑥𝑦 = 𝑒 + 1*6* 𝑒 = -164256.00 N/m.
Nous avons extraire les valeurs de : 𝑁 ′ 𝑥𝑥 , 𝑁 ′ 𝑦𝑦 , 𝑁 ′ 𝑥𝑦 , 𝑀𝑥𝑥 , 𝑀𝑦𝑦 et 𝑀𝑥𝑦 à partir de logiciel
cast3m et e = 0.20 m
1𝑓 𝛷
𝑙𝑠𝑚𝑎𝑥 = K*c+4 𝜏𝑐𝑡𝑚 𝜌 𝑠
𝑏𝑚𝑠 𝑠,𝑒𝑓𝑓
ℎ 𝛷𝑠
ℎ𝑐,𝑒𝑓𝑓 = Min [2 ; 2.5 ∗ (𝑐 + )] = 0.075 mm
2

91
Avec : h = 200 mm et c = 25 mm
𝑚
𝐴𝑠 (1 )
𝜌𝑠,𝑒𝑓𝑓 = 𝑙
= 0.019
ℎ𝑐,𝑒𝑓𝑓
1𝑓 𝛷𝑠
𝑙𝑠𝑚𝑎𝑥 = K*c+4 𝜏𝑐𝑡𝑚 𝜌 = 98.10 mm
𝑏𝑚𝑠 𝑠,𝑒𝑓𝑓
𝑓 1
Avec : K = 1 et 𝜏𝑐𝑡𝑚 = 18
𝑏𝑚𝑠
1 𝑁𝑥𝑦
Ɵ = 2arctg (2abs(𝑁 )) = 5.432 °
𝑥𝑥 −𝑁𝑦𝑦
(𝜀𝑠 − 𝜀𝑐 ) = 0.339*10−3*(cos (Ɵ)+sin (Ɵ))-0.264*10−4 = 3.450*10−4

𝝎 = 0.0677 mm

 EC2

𝑊𝑘 = 𝑆𝑟,𝑚𝑎𝑥 * (𝜀𝑠𝑚 –𝜀𝑐𝑚 ) <𝑊𝑙𝑖𝑚 = 0.3 mm

Avec :
𝑓𝑐𝑡,𝑒𝑓𝑓
𝜎𝑠 − 𝑘𝑡 𝜌 (1 + 𝛼𝑒 𝜌𝑝,𝑒𝑓𝑓 )
𝑝,𝑒𝑓𝑓 𝜎𝑠
𝜀𝑠𝑚 − 𝜀𝑐𝑚 = ≥ 0,6
𝐸𝑠 𝐸𝑆

𝜎𝑠 (T10) = 320 (MPa) (à partir de tableau 2.6 (chapitre 02))


25 2
𝐾𝑡 = 0.4, c = 25, 𝐾1 = 0.8, 𝐾2 = 0.75, 𝐾4 = 0.425 et 𝐾4 = 3,4 𝐶 3 = 2.267
𝐴𝑠 𝐴𝑠
𝜌𝑝,𝑒𝑓𝑓 = 𝑏 ℎ = (ℎ−𝑥) ℎ = 0.020
𝑐,𝑒𝑓𝑓 𝑏 𝑀𝑖𝑛[2.5(ℎ−𝑑); ; ]
3 2

𝑓𝑐𝑡,𝑒𝑓𝑓 = 3 (MPa)
𝐸 202000
𝛼𝑒 =𝐸 𝑠 = = 6.3125
𝑐𝑚 32000
𝜀𝑠𝑚 − 𝜀𝑐𝑚 = 1.250*10−3
𝐾 𝐾 𝐾 𝛷
𝑆𝑟,𝑚𝑎𝑥 = 𝐾3 𝑐 + 1𝜌 2 4 = 184.175 mm
𝑝,𝑒𝑓𝑓

𝑾𝒌 = 0.230 mm

 Formule propose par Etienne Gallitre (2008)

À l’aide de logiciel cast3m : ɣ = 4.27345*10−3 .


𝟒𝟎.𝟖𝟏 (𝟐) 𝜱
W = K*( 𝒇 ) 𝟑 *ɣ* 𝝆 = 0.52 mm.
𝒄𝟐𝟖
Avec : K = 1 et 𝑓𝑐28 = 40 MPa.

Maquette T10 Maquette T12


(Φ=10mm, 𝜌=0.6%) (Φ=12mm, 𝜌=0.11%)
Réglementation Ouverture des fissures Ouverture des fissures
(T10) (T12)
EC2
0.230 0.273

92
MC2010 0.0677 0.10
Norme SIA 0.347 1.136
SAFE 0.63 1.02
Numérique multicouche 0.903 0.85
Numérique
0.0584 0.1034
macroscopique

Tableau4. 2: comparaison entre différent réglementations avec les approches de


modélisation.
Le tableau (4.2) donne une comparaison pour les différentes formules proposées dans les
codes de calcul. On peut constater que les valeurs données par l’approche macroscopique sont
similaires avec ceux du modèle fibre (Mc2010). La norme SIA donne une augmentation
d’ordre de 33.71% par apport à l’EC 2 pour le voile T10 et une augmentation de cinq fois plus
par rapport à l’EC2 pour le voile T12. La formule proposée par SAFE donne des valeurs 3
fois plus grandes par apport de l’EC2. La modélisation multicouche surestime les valeurs de
fissuration.

4.5. Comparaison les ouvertures des fissures avec expérimental

Maquette T10 Maquette T12


(Φ=10mm, 𝜌=0.6%) (Φ=12mm, 𝜌=0.11%)
Cycles Dpl-max W W W Cycles Dpl-max W W W
mm (mm). (mm). (mm). (mm) (mm). (mm). (mm).
exp num- num- exp num- num-
mul- macro mul- macro
cou cou
T10-1 0.779 0.2 0.144 0.967 T12-1 0.246 0 0.0260 0.0027
T10-2 5.426 0.7 0.950 2.940 T12-2 0.85 0.3 0.156 2.624
T10-3 11.568 1.4 5.100 4.343 T12-3 2.6 0.7 0.430 5.73
T10-4 15.229 3.5 7.203 4.701 T12-4 7.8 1.7 4.49780 8.789

Tableau4. 3: comparaison les ouvertures des fissures (MC2010, macro, multicouche


et expérimentale) pour T10 et T12.

4.5.1. Discussion

D’après le tableau 4.3 on a constaté que les valeurs des ouvertures des fissures sont
très différentes. Le principal résultat quant peut extraire dans ce tableau c’est que l’ouverture
des fissures augmente lorsque le déplacement appliqué est important (une relation de
proportionnalité), les approches utilisées pour estimation des ouvertures des fissures
surestiment les valeurs de fissuration.

93
4.6. Comparaison effort-tranchant avec la réglementation internationale

Dans cette partie on va expliquer comment calculer les contraintes de cisaillement


selon les codes de conception. Pour le voile T10, pour calculer la contrainte de cisaillement
on a choisi trois réglementations : un règlement japonais (JEAG), américain (FEMA 365) et
un règlement européen (EC2) et après on a présenté deux tableaux récapitulatifs sur les
contraintes de cisaillement pour différents codes de calcul.

Pour les données des efforts internes (moment, effort normal) on les a extraits à l’aide du
logiciel CAST3M.
Les données nécessaires pour le calcul des contraintes de cisaillement sont
Effort normal : -35171N, Moment = 9424 N.m, Effort tranchant = 2.3936*105 .
La résistance caractéristique de béton a 28j :40MPa.
La limite élasticité d’acier : 500MPa.
ρh = ρv = 0.6%.
 JEAG :
σwy = 5000bars.
Fc = 400 bars.

−35171
σv = = −0.586 bars.
3000 ∗ 200
0.0679 ∗ 0.6%0.23 ∗ (400 + 180)
τu = + 2.7√0.6% ∗ 5000 + 0.1 ∗ (−0.586)
9424
√ + 0.12
2.3963 ∗ 105
= 32.5794bars.

𝛕𝐮 = 𝟑. 𝟐𝟓𝟕𝐌𝐏𝐚.

 FEMA 365 :

Vc = 0.25√30 ∗ 145.038 ∗ (200 ∗ 3000 − 3600)=11356605.52 Psi


Vs = 0.6% ∗ 500 ∗ 145.038 ∗ (200 ∗ 3000 − 3600)=259501989.6 Psi.
V=11356605.52+259501989.6=270858595.1 Psi.
270858595.1
V= 145.038 =18677500.897 N.
17677500.897
τu = =3.1125 MPa.
200∗3000

𝛕𝐮 = 𝟑. 𝟏𝟏𝟐𝟓 𝐌𝐏𝐚.

 EC2 :

0.18
CRd,c = 1.15= 0.157.
200
K = 1 + √(3000−30) = 1.25 ≤2.0 …CV.
40
fcd = 1.2=33.33 MPa.

94
25.104
σcp = 3000∗200 = 0.416 ≤ 6.666 …CV.
40
ν = 0.6 [1-250] =0.504.
Effort tranchant par apport les armatures de flexion :
0.5bw dνfcd =0.5*200*(3000-30)*0.504*33.34= 4990597.92.
2
VRd,c = [0.157 ∗1.25(100 ∗ 0.6% ∗ 40)1/3(0.5) + 0.150 ∗ 416]∗ 200 ∗ 3000=13 96039.095N
...CV
1396039.095
τuflex = 200∗3000 =2.327MPa.
Effort tranchant par apport les armatures de cisaillement :
VRd =1396039.095+1080*500=1936039.095 N… CV.
1936039.095
τucisai = 200∗3000 =3.226 MPa.
τuc = max(τucisai , τuflex) .

𝛕𝐮𝐜 = 𝟑. 𝟐𝟐𝟔𝐌𝐏𝐚.

4.6.1. Cas statique

Maquette T10 Maquette T12


(Φ=10mm, 𝜌=0.6%) (Φ=10mm, 𝜌=0.11%)
Règlement 𝜏𝑢 (Mpa). 𝜏𝑢 𝜏𝑢 (Mpa).
3.257 2.62
JAPAN
PS92 4.618 2.1043
EC2/EC8 3.226 1.72
BAEL 4.30 2.38
ACI 318-05 2.33 2.1165
UBT 3.93 2.42
FEMA365/ATC 3.1125 2.822

Tableau4. 4: comparaison les efforts de cisaillement statique pour différent code de


calcul.

4.6.2. Cas cyclique

Maquette T10 Maquette T12


(Φ=10mm, 𝜌=0.6%) (Φ=10mm, 𝜌=0.11%)
Règlement 𝜏𝑢 (Mpa). 𝜏𝑢 𝜏𝑢 (Mpa).
3.47 2.62
JAPAN
PS92 4.618 2.1043
EC2/EC8 3.89 1.72
BAEL 4.30 2.38
ACI 318-05 2.334 2.1168

95
UBT 3.93 2.42
FEMA365/ATC 3.1125 2.8224

Tableau4. 5: comparaison les efforts de cisaillement dynamique pour différent code


de calcul.
4.6.3. Commentaire

D’après le tableau 4.3 Et 4.4 On constate que les contraintes de cisaillement sont
identiques d’un code à l’autre et on a remarqué que les valeurs de cisaillement statique sont
les mêmes que les valeurs de cisaillement dynamique, car les réglementations donnent des
formules standards pour la vérification des efforts de cisaillement dans les voiles quel que soit
la sollicitation (dynamique ou statique) , et on a remarqué dans les calculs de l’effort de
cisaillement que les aciers transversaux participe à la contribution aux effort tranchant, et on
observe que le taux d’armature des aciers influe considérablement sur la résistance aux efforts
de cisaillement.

4.7. Conclusion

La procédure OUVFISS implémentée dans cast3M nous a permis de calculer les


ouvertures des fissures à différente échelle de modélisation. Le calcul numérique nous a
démontré que les zones critiques sont bien définies sur les cartes de fissuration pour le calcul
statique et dynamique en utilisant une approche macroscopique, par contre l’approche semi-
global ne donne pas les informations suffisant sur les ouvertures des fissures.
Pour les efforts de cisaillement, on a montré que la réglementation donne des formules
pratiquement identiques pour le dimensionnement des structures en béton armé sollicité aux
efforts de cisaillement.
Pour les ouvertures des fissures, on a trouvé une divergence des résultats entre l’expérimental,
le numérique et la réglementation.

96
Conclusion et perspective

Le travail de ce mémoire s’intéresse à l’analyse de comportement non-linéaire des


murs voiles sous chargement sévère (cas sismique). L’évaluation de la dégradation en termes
d’ouverture de fissures est un défi pour les constructeurs et les concepteurs.
Notre objectif était de pouvoir estimer l’endommagement et les ouvertures de fissures afin de
les comparer aux résultats expérimentaux et celles des codes de calcul sismique. L’objectif
final est de pouvoir proposer des améliorations à ces codes ou proposer des formules qui
n’existent pas pour le moment.
Le programme expérimental qui a servie à l’étude est le programme SAFE. L’avantage est
que nous disposons de tous les enregistrements et mesures nécessaire à une analyse
dynamique et statique.
Trois modélisations ont été réalisées Multifibres-Multicouches et Macroscopique.
Les approches simplifiées (Multifibres/Multicouches) sont capables de reproduire
correctement les réponses globales. Cependant, les réponses locales sont mal évaluées. La
modélisation macroscopique est celle qui donne les meilleurs résultats à l’échelle globale et
locale.
Les différents codes de conception donnent une estimation de l’effort de cisaillement ultime.
Une petite différence est observée entre ces codes. L’inconvénient c’est que la réglementation
donne des formules standard pour le calcul d’effort tranchant des murs voile soumis au
cisaillement. Quand à l’estimation des ouvertures de fissures, la plupart des codes ne
proposent pas de formules pour les ouvertures à l’exception de certains codes. La
modélisation numérique reste un outil nécessaire pour ce type de prédiction.

Les perspectives dans ce domaine sont nombreuses, on peut recommander les suivantes :
Contribution de l’effet échelles sur le comportement des murs voile soumise à des
chargements sévères.
Influence de l’inclinaison des armatures verticales suivant les contraintes principales sur la
vulnérabilité des murs voiles.
L’Etude du comportement des murs voile élancé dans le cas sismique.
Prise en compte de l’effet de l’adhérence acier –béton sur la réponse des murs voile dans
évènement sismique.

Références bibliographiques
97
 ACI 224R-01 Control of Cracking in Concrete Structures reported by
Committee 224 American Concrete Institute 2001.
 ACI318-05, ACI318R-05 Building code requirements for structural concrete
and commentary American Concrete Institute, 2005.
 ASCE-ACI Joint Task Committee. “Shear Strength of Reinforced Concrete
Members.” J. Struct. Engr. 1973. 1091–1187.
 ATC 40 Seismic evaluation and retrofit of concrete buildings Applied
Technology Council, November 1996.
 BAEL, 1999, Règles techniques de conception et de calcul des ouvrages et
constructions en béton armé suivant la méthode des états limites - BAEL 91
ré1-visé 99.
 Cédric Giry «Modélisation objective de la localisation des déformation et de la
fissuration dans les structures en béton armé.» Thèse de doctorat-Institue
national polytechnique de Grenoble en France année 2006.
 Chetra Mang «Modélisation de la liaison acier béton dans le calcul des
structures en béton armée.» Thèse de doctorat de l’université paris ouest
année 2015.
 D-Combescure rapport DM2S « Modélisation des structures de génie civil
sous chargement sismique à l’aide de CASTEM 2000 », 122 pages, année
06/03/2001.
 Etienne Gallitre «fissuration des voiles courte en béton armé soumis a des
sollicitations sismique-évaluation des ouvertures des fissures : contribution à
l’amélioration des connaissances, des données expérimentales aux pratiques
d’ingénierie.» Thèse de doctorat-Institue national des sciences appliquées de
Lyon en France année 2008.
 Eurocode 2 : Design of concrete structures Norme Européenne pr EN 1992-1-
1 (November 2002).
 Eurocode 8 : Design of structures for earthquake resistance Norme
Européenne pr EN 1998.
 FEMA 356 : Prestandard and commentary for the seismic rehabilitation of
buildings Federal Emergency Management Agency, November 2000.
 Ghezali sabah «modélisation numérique de la fissuration des structures en
béton armé : application à la simulation des enceintes de confinement.»
Magister-université technologie Tlemcen en Algérie année 2012.
 Jacky Mazars & al. Cracking in walls .In : Control of Cracking in Reinforced
Concrete Structures ,2015 ,p.93- 122, ISBN 978-1-78630-052-2.
 JEAG 4601 1987 : Technical Guidelines for aseismic Design of Nuclear Power
Plan (traduction NRC ref NUREG/CR -6241).
 Khuong Le Nguyen «contribution a la compréhension du comportement des
structures renforcées par FRP sous séisme.» Thèse de doctorat-Institue
national des sciences appliquées de Lyon en France année 2015.
 Labbé P., Noé H. (1992), “Safe, An Experimental Evidence Of The Frequency
Dependence Of Shear Wall Seismic Design Margins”, Proceedings of the 11th
World Conference on Earthquake Engineering, Canadian Association for
Earthquake Engineering.
 M. Matallah. " Modélisation numérique de l’endommagement anisotrope et
unilatéral dans les structures en béton». Thèse de Doctorat. Université de Pau
et des Pays de l’Adour,Octobre 2006.

98
 Medjahed Amina «approche semi global 3D pour le calcul des structures en
béton armé.» Magister-université technologie Tlemcen en Algérie année 2012.
 Mohamed Ibrahim «Modélisation et calcul des murs porteurs en béton armé
soumise aux effets sismique par la méthode des éléments finis : plaque
multicouche multicellulaire.» Thèse de doctorat-université Mohammed V
_agdal Rabat année 2013.
 Nicolae Ile «contribution a la compréhension du fonctionnement des voiles en
béton armé sous sollicitation sismique : apport de l’expérimentation et de la
modélisation a la conception.» Thèse de doctorat-Institue national des
sciences appliquées de Lyon en France année 2000.
 P. Pegon, G. Magonette, F. J. Molina, G. Verzeletti, T. Dyngeland, P. Negro,
D. Tirelli et P.Tognoli «Programme SAFE : rapport du test T10» Note
Technique du CCR no. I.98.71, 1998.
 P. Pegon, G. Magonette, FJ. Molina, G.Verzeletti, T.Duyneland, P.Negro,
D.Tirelli et P.Tognoli (1998). Programme SAFE : Rapports des tests T10 à
T12 European Commission, Joint Research Centre, ISPRA, technical notes.
 Panagiotis Kotronis «Cisaillement dynamique des murs en béton armé :
modèle simplifié 2D et 3D.» Thèse de doctorat-école normal supérieure de
Cachan en France année 2006.
 RAHMANI Abdallah Yacine «Contribution à la modélisation de la fissuration
dans les structures en béton par la méthode XFEM» Magister-université
technologie Tlemcen en Algérie année 2012.
 Saad Ammar «Modélisation du comportement, jusqu’à la rupture, en
cisaillement des voiles en béton armé.» Magister-université mouloud
Mammeri de Tizi-Ouzou en Algérie année 2013.
 Soufflet, M. &al. «Modélisation de structures voiles en béton armé soumise à
du cisaillement. »,2013.
 Uniform Building Code, 1997, Volume 2. Structural engineering design
provisions.
 Xuan Huy Nguyen «vulnérabilité des structures en béton armé à voiles
porteurs : expérimentation et modélisation.» Thèse de doctorat-Institue
national polytechnique de Grenoble en France année 2006.
 Xuan Huy Nguyen «vulnérabilité des structures en béton armé à voiles
porteurs : expérimentation et modélisation.» Thèse de doctorat-Institut
national polytechnique de Grenoble en France année 2006.

99
100

Vous aimerez peut-être aussi