These Iddir Abdelkader PDF

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
Ministère de l’Enseignement Supérieur et de la Recherche Scientifique
Université Mouloud Mammeri de TIZI-OUZOU
Faculté du Génie et de la Construction

Département de Génie Civil
THESE DE DOCTORAT
Spécialité : Génie Civil
Option : Structures et Matériaux
Présentée par
IDDIR ABDELKADER
Sujet
Modélisation des éléments de structure de section

Circulaire en béton armé confiné (Analyse de la
fissuration)
Devant le jury composé de :
Abdelkader TAHAKOURT Professeur, U. Béjaia Président

Mohand Said KACHI Professeur, UMMTO Directeur
Youcef BOUAFIA Professeur, UMMTO Co-Directeur
Mohamed KADRI Professeur, U. Boumerdès Examinateur
Abdellatif ZERIZER Professeur, U. Boumerdès Examinateur
Soutenue le 06 / 06 / 2016
REMERCIEMENTS
Une thèse de doctorat est par définition le fruit d’un travail personnel, mais qui pour ma part
n’aurait jamais abouti sans le concours d’un certain nombre de personnes que je souhaite
remercier ici.
Je remercie en premier lieu mon directeur de thèse Pr KACHI Mohand Said pour son
encadrement, ces précieux conseils, et de m’avoir fait profiter de son expérience et de son savoir-
faire ainsi que d’avoir mis à ma disposition tous les moyens nécessaires pour l’aboutissement de
cette thèse. J’ai beaucoup apprécié sa disponibilité et été touchée par son honnêteté intellectuelle.
Je souhaite exprimer ma reconnaissance à mon co-directeur de thèse le Pr BOUAFIA Youcef

qui a suivi de près mon travail tout au long de ces cinq années. Je le remercie particulièrement
pour son soutien et son aide. Je ne saurais oublier les orientations qu’il m’a indiquées pour faire
avancer ce travail et ses nombreux conseils qui ponctuent cette thèse. Je tiens à le remercier tant
pour ses compétences scientifiques que pour ses qualités humaines.
Je remercie le Pr TAHAKOURT Abdelkader qui me fait l’honneur de présider le jury, je

remercie également le Pr KADRI Mohamed et le Pr ZERIZER Abdellatif d’avoir accepté de juger
ce travail.
Je remercie amicalement tous les doctorants que j’ai rencontrés au cours de ma thèse pour leur
ouverture d’esprit et leurs amitiés.
Je remercie mes proches pour leurs soutiens et leurs encouragements.
Que tous ceux qui ont contribué de près ou de loin à l’aboutissement de cette thèse et trouvent ici
l’expression de ma profonde gratitude.
i
DEDICACES
A la mémoire de mon grand-père et de ma mère

À mon fils Aghiles
ii
Table des matières
INTRODUCTION GENERALE..........................................................................................1
CHAPITRE I ........................................................................................................................4
SYNTHESE BIBLIOGRAPHIQUE ....................................................................................4
I-1.Première partie : Le confinement ...................................................................................4
I.1.1. Introduction .................................................................................................................4
I.1.2. Principe du confinement .............................................................................................4
I.1.3. Types de confinement .................................................................................................5
I.1.3.1. Le confinement interne ........................................................................................5
I.1.3.2. Le confinement externe........................................................................................6
I.1.4. Poteaux mixtes ............................................................................................................7
I.1.5. Quelques modèles de comportement du béton confiné en élasticité non linéaire ....10
I.1.5.1. Modèle de Richart..............................................................................................10
I.1.5.2. Modèle d’Hognestad..........................................................................................10
I.1.5.3. Modèle Bousalem et al .[19]..............................................................................12
I.1.5.4. Modèle d’Ahmad et Shah ..................................................................................12
I.1.5.5. Modèle de Toutanji ............................................................................................13
I.1.5.6. Modèle Mander et al. .........................................................................................13
I.1.5.7. Modèle Saatcioglu et al. [25] .............................................................................14
I.1.5.8. Modèle Hoshikuma et al. [26] ...........................................................................15
I.1.5.9. Modèle de Fardis et Khalili ..............................................................................15
I.1.5.10. Modèle de Karbhari .........................................................................................16
I.1.5.11. Modèle de Mirmiran ........................................................................................16
I.1.5.12. Modèle de Xiao et Wu ....................................................................................17
I.1.5.13. Modèle de Harmon ..........................................................................................17
I.1.5.14. Modèle Fujikake ..............................................................................................17
I.1.5.15. Modèle de Beque ............................................................................................19
I.1.6. Discussion .................................................................................................................20
I.1.7. Modèle proposé par Bouafia.Y, Iddir. A et Kachi. M.S ...........................................20
I-2. Deuxième partie : Méthodes de calcul non linéaire ....................................................22
I.2.1. Introduction ...............................................................................................................22
I.2.2. Proposition de VIRLOGEUX et M’RAD.................................................................23
iii
I.2.3. Proposition de MULLER J [48]...............................................................................25
I.2.4. Proposition de Martins ..............................................................................................27
I.2.5. Proposition de Kachi.................................................................................................30
I.2.5.1 Domaine d’application et hypothèses générales.................................................31
I.2.5.2. Etude de l’équilibre d’une poutre en élasticité non linéaire ..............................32
I.3.Conclusion ....................................................................................................................33
CHAPITRE II.....................................................................................................................34
MODELISATION AVEC PRISE EN COMPTE DU CONFINEMENT.........................34
II.1 Introduction..................................................................................................................34
II.2 Domaine d’application et hypothèses générales..........................................................35
II.3. Comportement des matériaux .....................................................................................36
II.3.1. Béton comprimé ......................................................................................................36
II.3.2. Comportement du béton en traction ........................................................................37
II.3.3. Les aciers .................................................................................................................39
II.3.3.1 Comportement des aciers passifs.......................................................................39
II.3.3.2 Comportement des aciers actifs.........................................................................40
II.3.3. 3. Comportement des tubes en aciers ..................................................................41
II.3.4. Comportement du béton confiné .............................................................................42
II.3.4.1.Cas d’une section circulaire...............................................................................42
II.3.4.2. Cas d’une section rectangulaire........................................................................44
II.3.4.3. Cas d’une section circulaire confiné par tube..................................................46
II.3.4.4. Modèle ..............................................................................................................49
II.4. Modélisation des sections.......................................................................................51
II.4.1. Le béton ...................................................................................................................51
II.4.2. Les aciers .................................................................................................................51
II.5. Formulation du principe de calcul ..............................................................................52
II.5.1. Séquence de calcul...................................................................................................54
II.6. Élément poutre............................................................................................................57
II.6.1. Étude d’une section droite .......................................................................................58
II.6.2. Calcul des sollicitations internes .............................................................................60
II.6.3 Déformation de la section à l’équilibre ....................................................................65
II.6.4. Méthode itérative de résolution et critères de convergence ....................................66
II.6.5. Organigramme de la recherche de l’équilibre d’une section...................................67
II.6.6.Calcul de la rigidité de cisaillement de la section ....................................................68
II.6.7. Domaine d’application ...........................................................................................70
II.6.8.Contrainte transversale dans le béton .......................................................................71
II.7. Principes de la modélisation et équations générales...................................................71
iv
II.7.1. Milieux continus équivalents...................................................................................71
II.7.1. 1. En déformations...............................................................................................71
II.7.1.2. En contraintes ...................................................................................................72
II.7.2. Compatibilité partielle des déformations du béton avec les déformations moyennes
............................................................................................................................................72
II.7.3. Equations d’équilibre...............................................................................................73
II.7.3. 1. Equilibre local des couches .............................................................................73
II.7.3. 2. Equilibre global des sections ...........................................................................73
II.7.4. Lois de comportement des matériaux .....................................................................73
II.7.4.1. Comportement du béton en compression dans la direction D2.......................73
II.7.4.2. Comportement du béton en traction dans la direction D1...............................76
II.7.4.3. Comportement de l’acier .................................................................................77
II.7.5. Calcul de la distribution des contraintes tangentielles (y) ....................................77
II.7.6. Calcul de la distorsion moyenne.............................................................................79
II.7.7. Résolution du problème...........................................................................................79
II.7.7.1. Transformation des équations générales .........................................................79
II.7.7.2. Algorithme et organigramme de calcul ............................................................80
II.8.Expression de la matrice de rigidité de l’élément dans son repère intrinsèque.......83
II.9. Expression de la matrice de rigidité de l’élément dans le repère local ..................86
II.10. Expression de la matrice de rigidité de l’élément dans le repère global lié à la
structure ..........................................................................................................................92
II.11. Expression des accroissements des efforts en fonction des accroissements des
déplacements dans le repère global................................................................................94
II.12. Conclusion ............................................................................................................94
CHAPITRE III ...............................................................................................................95
APPLICATION DU MODELE ET CONFRONTATION DES RESULTATS ............95
III.1. Applications et confrontations dans le cas des sections en béton armé ................95
III.1.1. cas des sections circulaires .....................................................................................95
III.1.2. cas d’une section rectangulaire ..............................................................................99
III.1.3. cas d’une section béton enrobées de tube en aciers (CFST : Concrete filled steel
tubes) ................................................................................................................................101
III.1.4.Influence des paramètres de confinement sur le calcul des sections ....................103
III.1.4.1. Influence des armatures transversales ...........................................................103
III.1.4.2. Influence du diamètre des armatures transversales......................................106
III.1.4.3. Influence de l’espacement des armatures transversales ...............................107
Plus l’espacement des armatures transversales est réduit, on à un gain important de la
capacité portante...........................................................................................................109
III.1.4.4. Influence du diamètre des armatures longitudinales.....................................110
v
III.1.4.5. Influence de la résistance à la compression du béton ..................................111
III.1.4.6. Influence du diamètre de la section du béton...............................................113
III.2. Calcul global des poutres armées avec prise en compte du confinement ...........114
III.2.1 Essais du CEBTP – Poutres hyperstatiques HZ4..............................................114
III.2.2.Etude paramétrique ...........................................................................................118
III.3.Calcul en cisaillement d’une section de poutre....................................................120
III.3.1. Essais de l’Université de Toronto –Comportement de section .......................120
III.4.Etude de la fissuration ..........................................................................................125
III.5.Pieu en béton armé soumis à une flexion composée...........................................129
III.6. Conclusion...........................................................................................................131
CONCLUSION GENERALE ......................................................................................133
REFERENCES BIBLIOGRAPHIQUES .....................................................................135
vi
LISTE DES FIGURES
Figure I.1: béton armé confiné par cadres étriers transversaux [23-24]. .................................5
Figure I.2: béton armé confiné par cerce ou spirale en acier [23],[24]....................................6
Figure I.3: béton armé confiné par FRP et tube en acier. ........................................................7
Figure I.5 : Courbe contrainte-déformation du modèle Hognestad [18]................................11
Figure I.6: Courbe contrainte-déformation du modèle Bousalem et al [19]. .........................12
Figure I.7: Courbe contrainte-déformation du modèle Mander et al. ....................................14
Figure I.8: Loi de comportement réel du béton en compression proposée par ......................21
Figure I. 9 : courbe Charge – flèche obtenue par DEFLECT sur l’essai du CEBTP [49] .....26
Figure I.10: L’équilibre d’un voussoir. ..................................................................................29
Figure II.1 : Comportement du béton en compression loi de Sargin [38] .............................37
Figure II.2: Comportement du béton avant fissuration. .........................................................38
Figure II.3 : Contribution du béton tendu après fissuration A.Grelat [61]. ..........................38
Figure II.4: Comportement des aciers naturels. .....................................................................39
Figure II.5: Comportement des aciers écrouis. ......................................................................40
Figure II.6: Courbe contrainte-déformation des tubes en aciers ............................................41
Figure II.7 : section circulaire confinée par spirale en acier..................................................42
Figure II.8: béton effectivement confiné (cas d’une section circulaire confinée par spirale en
acier). ......................................................................................................................................43
Figure II.9 : béton effectivement confiné (cas d’une section rectangulaire confinée par
armatures transversales). .......................................................................................................45
Figure II.10 : béton effectivement confiné (cas d’une section circulaire confinée par tube en
acier). ......................................................................................................................................46
Figure II.11 :L’équilibre des forces dans une section circulaire confinée par tube en acier).
................................................................................................................................................48
Figure II.12 : Loi de comportement réel du béton en compression .......................................49
Figure II.13 : Section de béton découpée en tables trapézoïdales..........................................51
Figure II.14 : Représentation de la méthode de substitution utilisant les matrices cordes. ...54
Figure II.15: Description d’une section typique d’une poutre en béton armé et précontraint
................................................................................................................................................58
Figure II.16 : Organigramme simplifié de la recherche de l’équilibre de l’état d’équilibre
d’une section...........................................................................................................................68
Figure II. 17: Comportement du béton en traction.................................................................76
Figure II.18:Forces agissant sur une couche de béton entre deux sections [55]. ...................78
Figure II. 19:Organigramme général de l’étude d’une section. .............................................81
Figure II.20:Equilibre d’une couche à x et  connus. ...........................................................82
Figure II.21:Organigramme de la méthode simplifiée de calcul d’une section .....................82
Figure II.22: Repères d’un élément poutre.............................................................................83
Figure II. 23: Elément poutre déformé...................................................................................84
Figure III.1:Confrontation de la courbe théorique avec la courbe ........................................96
Figure III. 2:Confrontation de la courbe théorique avec la courbe .......................................96
Figure III. 3 : Confrontation de la courbe théorique avec la courbe .....................................96
Figure III.4: Confrontation de la courbe théorique avec la courbe .......................................97
Figure III.5: Confrontation de la courbe théorique avec la courbe .......................................97
vii
Figure III.6 : Superpositions des courbes théoriques  = f(e)...............................................98
Figure III.7 : Confrontation de la courbe théorique avec la courbe ......................................99
Figure III.8: Confrontation de la courbe théorique avec la courbe .....................................100
Figure III.9: Confrontation de la courbe théorique avec la courbe .....................................100
Figure III.10: Confrontation de la courbe théorique avec la courbe ...................................100
Figure III.11 : Confrontation des courbes théoriques avec les courbes ..............................101
Figure III.12: Confrontation des courbes théoriques avec les courbes ...............................102
Figure III.13: Confrontation des courbes théoriques avec les courbes ...............................102
Figure III.14: Confrontation des courbes Moment-Courbure calculées avec et sans
confinement ..........................................................................................................................104
Figure III.15: Confrontation des courbes Effort Normal -Déformation calculée sans et avec
confinement ..........................................................................................................................105
Figure III.16 : Confrontation des courbes frontières efforts normal -moment calculée ......105
Figure III.17: Confrontation des courbes Moment-Courbure confinée en T8 et T10..........106
Figure III.18 : Confrontation des courbes Effort Normal –Déformation confinée en T8 et
T10........................................................................................................................................107
Figure III.19: Confrontation des Courbes Moment-Courbure confinée avec espacement
variable .................................................................................................................................108
Figure III.20: Confrontation des Courbes Effort Normal –Déformation.............................109
Figure III.21:Confrontation des courbes d’interactions effort normal -moment confinée
avec espacement variable .....................................................................................................109
Figure III. 22 : Confrontation des courbes Effort Normal -Déformation calculée avec les
différents diamètres ..............................................................................................................110
Figure III. 23: Confrontation des courbes d’interaction avec les différents diamètres........111
Figure III.24: Confrontation des courbes Moment -Courbure calculées avec les différentes
classes du béton ....................................................................................................................112
Figure III.25: Confrontation des courbes Effort Normal -Déformation calculée avec les
différentes classes du béton..................................................................................................112
Figure III.26: Confrontation des courbes d’interactions calculées avec les différentes classes
du béton ................................................................................................................................113
Figure III.27: Confrontation des courbes d’interaction avec les différents diamètres de la
section...................................................................................................................................114
Figure III.28: Dimensions et détails de ferraillage des poutres HZ [70]. ...........................115
Figure III.29 : Disposition des armatures longitudinales pour les poutres HZ [70]. .........115
Figure III.30 : Evolution de la flèche en fonction de la charge pour les poutres HZ...........117
Figure III.31: Evolution de la flèche en fonction de la charge pour les poutres HZ4 pour des
espacements 10, 15,20 et 30cm ............................................................................................119
Figure III.32 : Comparaison des résultats du calcul aux résultats expérimentaux et à ceux124
Figure III. 33 : caractéristiques géométriques et détails de ferraillage de la poutre HZ4....125
Figure III.34: Distribution des déformations longitudinales moyennes sur la hauteur de la
section avant fissuration de la poutre HZ4...........................................................................126
section au début de la fissuration de la poutre HZ4 .............................................................126
section à la rupture la poutre HZ4 ........................................................................................127
Figure III.37: hauteur décomprimée de la section en fonction de l’espacement a la rupture
de la poutre HZ4...................................................................................................................128
Figure III.38 : Evolution de la flèche en fonction de la charge pour le pieu de Zhan. ........130
viii
LISTE DES TABLEAUX
Tableau III.1:Caractéristiques des données expérimentales des sections circulaires en béton.

................................................................................................................................................95
Tableau III.2: Caractéristiques des données expérimentales des sections circulaires en béton.
................................................................................................................................................98
Tableau III.3: Caractéristiques des données expérimentales des sections rectangulaires en
béton. ......................................................................................................................................99
Tableau III.4: Caractéristiques des données expérimentales des sections circulaires en béton
..............................................................................................................................................101
Tableau III.5: Caractéristiques du béton section 40cm armé avec 8T12 et des armatures ..104
Tableau III. 6: Caractéristiques des aciers. ..........................................................................104
Tableau III. 7: Caractéristiques du béton section 40cm armé avec 8T12 et des armatures .106
Tableau III.8 : Caractéristiques des aciers. ..........................................................................106
Tableau III. 9 : Caractéristiques du béton section 40cm armé avec 8T12 et des armatures
transversales T8 avec un espacement variable.....................................................................108
Tableau III.10: Caractéristiques des aciers. .........................................................................108
Tableau III.11 : Caractéristiques du béton (section 40cm) armé avec 8T12 ,8T14 ,8T16 et
des armatures transversales T8.............................................................................................110
Tableau III.12: Caractéristiques de la section en béton (section 40cm) avec les différentes
classes ...................................................................................................................................111
Tableau III.13: Caractéristiques du béton avec les différents diamètres de la section ........113
Tableau III.14 : Caractéristiques du béton et des armatures passives [70]. .........................116
Tableau III.15 : Comparaison des valeurs ultimes expérimentales et calculées de la charge Q
Pour la poutre HZ4. ..............................................................................................................117
Tableau III.16. Caractéristiques du béton et des armatures passives poutre HZ4 [70]........118
Tableau III. 17: Comparaison des valeurs ultimes expérimentales et calculées de la charge Q
..............................................................................................................................................118
Tableau III.18 : Propriétés des différentes poutres [71],[63],[66]. .....................................121
Tableau III. 19: Comparaison des valeurs ultimes expérimentales et calculées de l’effort
tranchant. ..............................................................................................................................124
Tableau III.20 : Résultats calculés des efforts au début des premières fissures de la Poutre
HZ4.......................................................................................................................................128
Tableau III.21 : Résultats calculés des efforts à la rupture de la Poutre HZ4.....................129
Tableau III.22 : Caractéristiques du béton et des armatures passives [72]. .........................129
Au début de la première fissure et à la rupture. ...................................................................130
ix
LISTE DES SYMBOLES
Aaj : section d’une barre d’ordre j

Ac : section interne comprise entre les spirales.
Ae : section effectivement confinée.
Asl : section d’armature longitudinale.
Ast : section d’armatures transversale.
Asp : section des aciers transversaux.
bj : base inférieure du trapèze d’ordre i
bj+1 : base supérieure du trapèze d’ordre i+1
b(y) : base de la section à une côte y
ds : diamètre de la section pris entre centre de gravité d’une spirale
Eb0 : module d’élasticité longitudinale.
Esc : Paramètre de la loi de comportement qui régit la partie descendante poste pic.
Es : module élastique de l’acier
fcc : la contrainte du béton confiné.
fc0 : la contrainte du béton non confiné à l’âge j
fl : pression latérale du fluide.
f’l : pression de confinement latérale effective.
fyh : contrainte limite des aciers.
Fcj : résistance à la compression du béton.
Ftj : résistance à la traction du béton
fe : Limite élastique de l'acier
IG : moment d’inertie de la section
IGH : moment d’inertie de la section homogénéisée
kc : Paramètre de la loi de Sargin avec les paramètres de confinement
ke : coefficient de confinement effectif.
k : paramètre qui ajuste la branche ascendante de loi de la loi de sargin
k’ : paramètre qui ajuste la branche descendante de loi de sargin
k’c : RC: Paramètre de la loi de SARGIN avec les paramètres de confinement
k1, k2 : coefficients, fonctions du mélange du béton et de la pression latérale.
Kb : paramètre de la loi de Sargin
K’b : paramètre de la loi de Sargin
Miext : moment initial extérieur
x
Mint : moment intérieur
Niext : effort normal initiale extérieur
Nint : effort normal intérieur
Nbt : nombre de barres d’aciers
Ne : nombre de trapèzes
s : distance entre les extrémités de deux spirales adjacentes.
s’ : distance entre les centres de gravité de deux spirales adjacentes.
S : aire de la section
Sh : aire de la section homogénéisée
yaj : ordonnée de la barre d’ordre j
YG : ordonnée du centre de gravité de la section
yj : ordonnée de la fibre inférieure du trapèze d’ordre i
yj+1 : ordonnée de la fibre inférieure du trapèze d’ordre i+1
YGH : ordonnée du centre de gravité de la section homogénéisée
w’i : la distance entre deux armatures longitudinale adjacente
N : accroissement de l’effort normal
M : accroissement du moment fléchissant
V : accroissement de l’effort tranchant
u : déformation axiale d’une section
w : rotation différentielle d’une section
 : Déformation du béton non confiné.
εc : déformation du béton confiné
εcc0 : la déformation correspondante a la contrainte fcc du béton confiné
ε0 : déformation correspondante à fc0.
rt : déformation de rupture du béton en traction
e : déformation élastique limite de l’acier
a : déformation de l’acier
cu : déformation de rupture du béton en compression
ccu : déformation de rupture du béton confiné en compression
au : déformation de rupture de l’acier passif
i : déformation de la fibre inférieure du trapèze d’ordre i
i+1 : déformation de la fibre supérieure du trapèze d’ordre i +1
ft : déformation de fissuration du béton
65 : déformation du béton confiné correspondant à 0.65fcc
xi
η : coefficient d’équivalence acier-béton
 : Opérateur permettant le passage des déformations aux efforts internes
bt : restriction de l’opérateur  au béton
at : restriction de l’opérateur  aux aciers
() : loi de comportement d’un martiaux
b() : loi de comportement du béton
a() : loi de comportement de l’acier passif
b, s : coefficients de sécurité
e : limite élastique de l’acier passif
σc : contrainte du béton en compression.
σcc : contrainte du béton confiné en compression
ρcc : rapport entre la section des armatures longitudinales et la section du béton.
ρs : rapport entre le volume des armatures de confinement transversale et le volume du béton
confiné.

Fs : Incrément de la force extérieure
[Ks] : Matrice de rigidité de la section
θ : Iinclinaison de la contrainte principale de compression
∆u : Incrément de déformation axiale
∆w :Iincrément rotation différentielle
∆moy : Incrément de la distorsion moyenne
σc1 : Contrainte principale de traction.
σc2 : Contrainte principale de compression.
1 : Déformation principale de traction.
2 : Déformation principale de compression.
 : Contrainte de cisaillement.
∆i : Incrément de la contrainte de cisaillement.
Ps : Action de précontrainte initiale
∆ : Incrément de la déformation de la section.
xii
Résumé
Dans le domaine de génie civil le calcul des sections circulaires n’est pas aisé compte tenu de la
disposition des armatures qui induit plusieurs inconnues dans les équations d’équilibre.
L’expérience acquise, les avantages et les inconvénients ainsi que les problèmes à résoudre ont été
largement discutées par de nombreux auteurs. Le calcul de telles sections présente des incertitudes
sur certains aspects de leur comportement. L’effet de l’influence du béton tendue même après
fissuration est d’une grande importance. L’effet confinement du béton en général lui confère une
plus grande résistance et une plus grande ductilité en fonction du niveau du confinement. Dans le
cadre de cette étude nous proposons un modèle de calcul basé sur les équations d’équilibres et de
compatibilités qui sont formulées en termes de contrainte et déformation moyenne selon des
hypothèses. Ce modèle est capable d’analyser des sections quelconques avec tous les détails de
ferraillage sous l’ensemble des sollicitations (effort normal, moment fléchissant et effort tranchant).
Une loi de comportement contrainte- déformation du béton confiné a été proposée. Le confinement
du béton sera introduit par le biais d’un coefficient de confinement effectif, et d’une contrainte de
confinement latérale. La méthode de calcul est injectée dans un logiciel de calcul non linéaire des
éléments de structures afin d’étudier la contribution et l’influence du confinement du béton sur le
comportement global des éléments de structures de section circulaire et sur la rigidité au
cisaillement dans le cas des poutres en béton armé et précontrainte. Enfin, le modèle de prise en
compte de l’effet de l’effort tranchant dans le domaine non linéaire a été également proposé, ce
dernier permet de prédire le comportement des poutres en cisaillement combiné à l’existence des
contraintes normales, ce qui permet alors d’estimer la rigidité réelle de cisaillement des poutres,
pour un niveau de charge donné.
Mots-clés : Élasticité non linéaire, Confinement, section circulaire, effort tranchant, rigidité au
cisaillement, armatures discrètes, équations d’équilibres, équations de compatibilités, rupture.
xiii
Abstract
In the civil engineering, the calculation of circular sections is not easy because there are many
unknowns in the equations of equilibrium. The required experience, advantages and disadvantages
and problems to solve have been widely discussed by many authors. The calculation of such
sections is difficult to predict exactly there behavior. The effect of the tensioned concrete, even after
cracking, has a great importance. The effect of the confinement of concrete generally gives this one
greater resistance and ductility according to the level of confinement. In this study we propose a
model which is based on the equilibrium, compatibility equations and whish are formulated by
average stresses and average strains. According to the assumption, this model is able to analyze
sections having unusual forms or reinforcing details, loaded in combined bending, axial load and
shear. In the case of the confined concrete, a stress-strain relation of this one is proposed. In this
case, for the purpose of investigating confinement effect, an effectively confined ratio was
introduced and the effects according to concrete compressive strength. Predictions of the model are
compared with several calculated unconfined reinforced or prestressed beams and are shown to
estimate the effect of the confinement on the shear stiffness of the sections in the case of the elastic
linear domain, after concrete cracking and after reinforcement yielding. The results of the
comparison give a good satisfaction.
Keywords: nonlinear elasticity, confinement, circular section, shear force, shear stiffness, ties
discrete, equilibrium, compatibility, cracks.
xiv
INTRODUCTION GENERALE
Dans le domaine de génie civil le calcul des sections circulaires n’est pas aisé compte tenu
de la disposition des armatures qui induit plusieurs inconnues dans les équations d’équilibre.
L’expérience acquise, les avantages et les inconvénients ainsi que les problèmes à résoudre ont été
largement discutées par de nombreux auteurs.
Le calcul de telles sections présente des incertitudes sur certains aspects de leur
comportement. Un effort de recherche tant théorique qu’expérimental est donc nécessaire pour
permettre de les comprendre et de les résoudre.
L’effet de l’influence du béton tendue même âpres fissuration est d’une grande importance.
De plus, l’influence de la prise en compte du confinement dans la modélisation et le calcul

des structures en béton armé et en béton précontraint permet de mieux approcher le comportement
réel. Et, la prise en compte de l’effet de l’effort tranchant sur le comportement des structures est
aussi d’une grande importance. Dans le cas particulier des poutres, cet effet est pris en compte par
l’adoption du modèle de calcul en cisaillement des poutres développé par Kachi et al [52], [55-60] .
L’effet confinement du béton en général lui confère une plus grande résistance et une plus
grande ductilité en fonction du niveau du confinement. Plusieurs études ont été effectuées sur le
confinement du béton [1], [18-37]. L’ensemble des modèles proposés dans ces études, pour prédire
le comportement du béton confiné, présentent des coefficients de confinement difficilement
maitrisables.
Dans le cas particulier des sections circulaires, les abaques de Davidovicci, basées sur les
principes du BAEL91 et de l’EUROCODE2 et d’une répartition uniforme de l’acier sur toute la
section, permettent de déterminer la section d’armatures longitudinales en flexion composée dans le
cas d’un effort normal de compression. Une méthode de calcul par les armatures discrètes permet le
calcul des armatures longitudinales et ce qu’elle que soit la nature de l’effort normal. Le
confinement interne du béton y est introduit et des confrontations ont été faites sur des sections et
sur des poutres en béton armé et en béton précontraint.
Dans le cadre de cette étude nous proposons une méthode de calcul des éléments de section
circulaire basée sur la méthode de calcul par les armatures discrètes développée par Kachi et al [52],
1
[55-60] en y introduisant différentes formes de confinement internes et externe. L’étude sera menée
en élasticité non linéaire, il sera tenu compte du comportement non linéaire des différents matériaux
constituant la section transversale. Le confinement du béton sera introduit par le biais d’une loi de
comportement également non linéaire, du béton confiné et du béton armé confiné tenant compte du
niveau de confinement, du type de confinement et de l’ensemble des sollicitations (effort normal,
moment fléchissant et effort tranchant). Enfin, la méthode de calcul sera injectée dans un logiciel de
calcul non linéaire des éléments de structures afin d’étudier la contribution et l’influence du
confinement du béton sur le comportement global des éléments de structures de section
quelconques.
En ce qui concerne les non linéarités, le calcul proposé prend en compte deux types de non
linéarités :
- la non linéarité mécanique du comportement réel des matériaux
(Fissuration, plasticité, etc),
- les non linéarités géométriques, liées au comportement d’ensemble de la
structure (effets de second ordre).
Le modèle de prise en compte de l’effet de l’effort tranchant dans le domaine non linéaire
développé par Kachi et al [52], [55-60] a été adopté, ce dernier permet de prédire le comportement
des poutres en cisaillement combiné à l’existence des contraintes normales, ce qui permet alors
d’estimer la rigidité réelle de cisaillement des poutres, pour un niveau de charge donné.
La présente étude est structurée comme suit :

- au chapitre 1 est présentée d’une part une étude bibliographique portant sur le
confinement du béton. Cette étude porte principalement sur les types de confinement internes et
externes, une comparaison entre les différentes méthodes de calcul à ce sujet et quelques modèles
de confinement également à la fin de cette première partie. D’autre part, des méthodes de calcul non
linéaire et des études, portant sur le comportement des éléments en béton armé et précontraint, sont
aussi présentées dans ce chapitre.
- le chapitre 2 est consacré au développement du modèle de calcul non linéaire du

comportement des poutres en béton armé et précontraint en cisaillement ainsi que le comportement
des différents matériaux est aussi développé dans ce chapitre.
2
- le chapitre 3 porte sur les résultats de calcul mené sur les sections et sur des éléments de
structure en béton confiné. Une suite de résultats est présentée aussi bien sur les sections sous forme
de loi de comportement reliant la contrainte à la déformation proposée comparée à celle mesurée
expérimentalement. Nous avons présenté également les résultats obtenus dans l’étude du
comportement globale des poutres en béton armé et en béton précontraint, une comparaison aux
résultats mesurés expérimentalement a été également effectuée.
- Enfin on termine par une conclusion de la présente étude.
3
Chapitre I Synthèse Bibliographique
CHAPITRE I
SYNTHESE BIBLIOGRAPHIQUE
I-1.Première partie : Le confinement
I.1.1. Introduction
Le béton est un matériau dont les performances en traction sont faibles par rapport à celles
de la compression. Pour l’utilisation efficace d’un tel matériau, il a été nécessaire de le
renforcer par des aciers longitudinaux qui reprennent les efforts de traction .Ce renforcement
a permis le développement du béton armé dans le domaine de la construction. Aussi, en plus
de ces armatures longitudinales, le béton a été renforcé avec d’autres types d’armatures,
appelées ; armatures transversales ou armatures de confinement (c’est le béton confiné).
I.1.2. Principe du confinement

C’est une autre appellation du béton armé renforcé par des armatures transversales ou par
jaquette externe qui exercent des pressions latérales sur le béton. Les études récentes ont
prouvé que le comportement du béton en compression dépend de sa richesse en confinement ;
généré par une disposition adéquate des armatures transversales ou une épaisseur
conséquente de la jaquette externe. Il en résulte par voie de conséquence un accroissement
important de la résistance en compression ainsi que la déformation correspondante.
4
I.1.3. Types de confinement

Actuellement le confinement des sections est une méthode efficace pour augmenter la
capacité portante des sections en béton armé ou en béton de fibres ; on distingue deux types
de confinement à savoir :
I.1.3.1. Le confinement interne

Le confinement interne se fait par des cerces en aciers ou spires dans le cas des sections
circulaires et des cadres ou étriers transversales dans le cas des sections rectangulaires. Le
confinement interne du béton peut être aussi réalisé avec des fibres de carbone ou de verre-
époxydes. Le Béton confiné ne sera pas considéré comme une structure, mais en tant que
matériau multiple (composite), qui sera défini par les propriétés globales.
Figure I.1: béton armé confiné par cadres étriers transversaux [23-24].
5
Figure I.2: béton armé confiné par cerce ou spirale en acier [23],[24].
I.1.3.2. Le confinement externe

Les fibres composites polymères de renforcement FRP (Fibre renforced polymer) sont
maintenant très utilisées pour une variété d’applications. Elles sont particulièrement
attrayantes pour le renforcement et la réadaptation des structures en béton. Il est prouvé
qu’elles sont souvent utilisées pour le confinement des colonnes en béton. En effet, de
nombreuses études ont clairement démontré qu’envelopper avec des feuilles de FRP des
colonnes circulaires est une technique très efficace pour augmenter leur résistance et ductilité.
Le confinement latéral est une méthode bien connue pour augmenter la capacité des
colonnes renforcées pour porter les charges. Dès 1929, Richart et al. [1] ont prouvés qu'une
contrainte de compression latérale peut mener à un perfectionnement significatif de la
résistance à la compression du béton.
Le confinement externe peut être fourni par divers moyens, tels que ; les vestes, les
cercles, l'emballage en acier de FRP (Fibre renforced polymer) autour de la colonne ou par
6
tubes en acier (CFST : Concrete filled steel tubes); de nombreux modèles analytiques sont
proposés pour prédire un tel confinement.
Dans le cas du confinement externe par tube en acier, aux premières phases du
chargement, le coefficient de Poisson du béton est inférieur à celui de l'acier, et le tube en
acier n'a aucun effet sur le noyau béton. Par conséquent, l'expansion latérale du béton devient
graduellement plus grande que celle de l'acier. Une pression radiale se développe à l'interface
acier le retenant ; de ce fait le noyau béton en développant une tension circulaire dans le tube.
À ce stade, le noyau béton est soumis aux contraintes triaxiales et le tube en acier aux
contraintes bi axiales, de sorte qu'il y ait un transfert de charge à partir du tube au noyau,
comme le tube ne peut pas soutenir l'effort de flexion longitudinalement en présence d'une
tension circulaire. La charge correspondant à ce mode de rupture peut être considérablement
plus grande que la somme de celle de l'acier et du béton, mais la rupture de cisaillement peut
intervenir avant que le transfert de charge soit complet [2].
Figure I.3: béton armé confiné par FRP et tube en acier.
I.1.4. Poteaux mixtes

I.1.4.1. Poteaux mixtes remplis de béton [2]
Les poteaux mixtes remplis de béton ont été l’objet de plusieurs études .Ce genre de
structure à été largement discuté par de nombreux auteurs .Y.C.Wang [3] à proposé une
méthode basée sur le code Britannique BS5950 en faisant l’étude sur plusieurs type de
7
poteaux remplis de béton. Vu la souplesse de sa méthode Y.C.Wang [3] à conclu qu’elle

prédit d’une façon claire le comportement des sections mixte et elle est semblable à celle
proposée dans l'EC4 et meilleure que celle du BS 5400.
Kefeng.T, John M.N, et Xincheng P [4] ont étudiés des colonnes tubulaires de section
circulaire remplis de béton. Les résultats expérimentaux ont montrés que le confinement du
béton par tube en acier peut de manière significative améliorer les propriétés mécaniques de
la colonne mixte comparée à celle en béton armé. Cette amélioration était directement
proportionnelle au coefficient de confinement.
Les travaux de M. Mouli, et H. Khelafi [5],[8] basé sur l’expérimental concerne les
poteaux rectangulaire en acier remplis de béton ordinaire. Les propriétés du béton de
remplissage utilisé étaient le facteur important de l’étude. La confrontation des résultats
obtenus été satisfaisante. Il on est de même pour les travaux de Dalin Liu [6], [9], Ehab. E,
Ben Young. D, et Dennis. L [7] qui ont pris pour leurs travaux des bétons de haute résistance.
Les poteaux mixtes de section rectangulaire et carrée avec raidisseurs ont pris également
une grande partie dans les projets de recherches avancés, dans le but de mieux prédire le
comportement de ces éléments sous n’importe quel type de chargement ou sollicitation [9]
Dans la pratique, des colonnes sont soumises à la compression axiale et aux moments de
flexion. Dans l’étude de telle structure non seulement les propriétés matérielles non linéaires
qui doivent être prises en considération mais également les effets géométriques non linéaires.
I.1.4.2. Effets de chargement, et de sollicitations extérieurs [2]
Le comportement d’un poteau mixte soumis à des sollicitations externes est conditionné
par plusieurs facteurs citant : le type du poteau, son élancement, sa section droite, et le rapport
entre la longueur du poteau et les dimensions de sa section droite.
Dans la conception, le comportement interactif pour des poteaux de longueur
intermédiaire est pris en considération par une constante empirique. La prévision de tels
constants en fonction de la géométrie, les propriétés de matériaux, et les imperfections du
poteau était l'objet de recherche d’Ever .J. Barber [11]. Qui a trouvé que les sections mixtes
restent ainsi linéaires pour de grandes contraintes en tenant compte de l'interaction entre les
différents modes de voilement [11-12] et de l’élancement des poteaux.
8
Le comportement des poteaux est défini par l’approche d’EULER, si les imperfections
sont négligés, quand en augmente la longueur du poteau on risque d’avoir un problème
d’instabilité a cause de l’augmentation de la valeur du moment fléchissant secondaire à mi-
hauteur, ce qui amène à la rupture du poteau par moment fléchissant, et non pas par
compression.
L’influence des charges axiales de compression sur les sections droites des poteaux courts
est moins importante comparée aux poteaux élancés [13].
Actuellement, les poteaux mixtes de section circulaire sont de plus en plus employés dans
la construction des bâtiments modernes. L’information disponible sur les modèles de
comportement de ce type de structure sous les différentes conditions de chargement est
limitée. Les poteaux mixtes de section circulaire sous des charges axiales, et des charges de
flexion cycliquement croissantes ont été l’objet de certains travaux de recherche ou plusieurs
paramètres ont été inclus dans l’étude.
Dans l’EC4, Un model mécanique a été développé sur des poteaux [14], et ce pour le
calcul du coefficient de ductilité des poteaux mixtes sous chargement cyclique.
I.1.4.3. Tube en acier [2]

La totalité des sections classiques des poteaux tubulaires remplis de béton sont produites
de l’acier de différentes nuances, et limites élastiques, des aciers de haute résistance pour
améliorer le comportement des poteaux mixtes, mais dernièrement on a pu développer
d’autres matériaux de construction dont l’objectif est d’améliorer les propriétés physiques,
chimiques, ainsi que comportement mécanique des tubes produits à base de ces matériaux.
Lorsque le béton tend à prendre de l'expansion dans la direction transversale au
chargement dû à l'effet de Poisson, il engendre des efforts dans l'acier qui s'oppose à cette
expansion. Le confinement qui est bénéfique pour le béton peut donc se révéler néfaste pour
la section d'acier. Dans le cadre de la plupart des projets de recherches, le niveau de
confinement que fournit la section d'acier au béton est complexe ; car il n'est pas uniforme
variant selon les directions. Il dépend également d'un ensemble de paramètres pouvant varier
selon le niveau de chargement. Il n'existe donc pas de solution simple pour évaluer ce niveau
de confinement et son influence sur le béton.
Les sections creuses circulaires fournissent une quantité significative de confinement tandis
que cet effet est négligeable dans le cas des sections rectangulaires. La force supplémentaire
se produit en raison de l'augmentation de la résistance à la compression du noyau béton qui
est retenu latéralement par le tube en acier dans le cas des sections circulaires. L’effet
9
d'emprisonnement n'est pas présent dans les sections creuses rectangulaires remplies de béton,
excepté dans les régions de coin, ou une tension de cercle se développe sur les murs latéraux
[15-17].
I.1.5. Quelques modèles de comportement du béton confiné en élasticité non linéaire
I.1.5.1. Modèle de Richart

Professeur Rechart, F.E. de l'Université de l'Illinois à Urbana-Champaign est un pionnier dans
la recherche du béton confiné [1]. Dans ses premiers travaux sur l'effet de confinement
transversal sur le comportement du béton en compression [1928,1929], Richart à découvert
que la résistance du béton à confinement actif par pression latéral hydrostatique était
approximativement la même que pour le béton avec une pression de confinement passive par
spirales en acier circulaires rapprochées provoquant une pression latérale équivalente. La
force et la déformation longitudinale correspondante à la résistance du béton confinée par une
pression de fluide hydrostatique active peuvent être représentées par les relations simples:
'  f ' k f
f cc (1.1)
co 1 l
 f 

 cc   co 1  k 2 l  (1.2)
 f' 
 co 
Où : f'cc et cc sont la résistance du béton confiné et la déformation correspondante,
respectivement. f'co et co sont la résistance du béton non confiné et la déformation
correspondante. fl est la pression latérale de confinement. k1 et k2 sont des coefficients. Basé
sur les résultats des essais, Richart et al. [1], ont trouvé que les valeurs moyennes des
coefficients k1 et k2 sont égales à k1 = 4.1 et k2 = 5k1.
I.1.5.2. Modèle d’Hognestad

Le modèle d’Hognestad, (1951), [18] permet de reproduire l’adoucissement après fissuration,
le comportement unilatéral du béton (fermeture des fissures) ainsi que l’adoucissement après
avoir atteint la résistance en compression et la prise en compte du confinement sur cet
adoucissement et la résistance en compression. Deux fonctionnements d’ouverture et
fermeture de fissures sont possibles, (Figure I.5).
10
Figure I.4 : Courbe contrainte-déformation du modèle Hognestad [18].
Sous chargement monotone, la loi est du type parabole puis linéaire avec adoucissement. Un
plateau avec une contrainte résiduelle non nulle peut-être défini après l’adoucissement. En
utilisant les conventions du signe (déformation positive en compression) :
    
Pour 0<e<ec0   2  
c 0 c 0 c 0

Pour ec0 <e  1  z .    c 0 
c 0
Avec:
- c0  2 . c 0
et de plus il faut prendre la contrainte résiduelle ;
E
-  pt = 20% de la contrainte maximale ;
- (ec0, c0) Sont les coordonnées du pic de la courbe de compression ;
- Z : pente après le pic ;
-  p t : contrainte résiduelle.
11
I.1.5.3. Modèle Bousalem et al .[19]

f cc xn
 cc  (Partie ascendante)
n 1  xn
 cc  f cc  E s  c   cc 0   0.3 f cc (Partie descendante)
4  f c20
Es 
ke   s  f yh
Figure I.5: Courbe contrainte-déformation du modèle Bousalem et al [19].
I.1.5.4. Modèle d’Ahmad et Shah

S.H. Ahmad et S.P. Shah [20-21] ont proposé leur modèle pour prédire la courbe contrainte-
déformation du béton confiné basé sur les propriétés des spirales de renforcement et la
relation constitutive du béton ordinaire.
Pour les essais de pression de confinement actifs, des échantillons de béton sont soumis à une
augmentation de la contrainte axiale tandis que la pression de confinement latéral est
maintenue constante. Dans le modèle de Shah, pour un béton soumis à un état général de
compression triaxiale, la relation contrainte- déformation est :
y
A x  D 1 x2
i i i

i

i

1 A  2 x  D x2
i i i

Ou : y=fl /fcc et x1=i/ci,
f1 = la plus grande contrainte principale de compression,
fp = le plus grande résistance principale à la compression,
i = la déformation dans la ième direction principale (i = 1, 2 ou3),
ip = la déformation au pic dans le Ième direction,
Ai = Ei / Eip, Ei est la pente initiale de la courbe contrainte-déformation;
12
Eip = fp/ip. Di est le paramètre qui régit la partie descendante de la courbe contrainte-
déformation.
La résistance à la compression fp a été déterminée par un critère de résistance basé sur une
théorie octaédrique.
 oct 
 0 .2261  0 .7360 oct
f f
o o
L’équation équivalente est:
'  f '  4 .255 f

f cc co l
I.1.5.5. Modèle de Toutanji

Un modèle d’Ahmad et Shah modifié a été proposé par Toutanji (1998) [22].
A
f 
1  B  C 2
 0 . 84 
  2 tf  
'  f ' 1  2 .2  com  
f cc 
c
 df '  
  c  
 
Où : fcom est la résistance à la traction du composite confiné.

Sur la base des résultats des tests, la relation modifiée de Rishart a été utilisé dans Modèle de
Toutanji.
  f' 

  cc
 cc   co 1  537  l  2.6   
 1 
  f' 
  c 
I.1.5.6. Modèle Mander et al.

Mander et al. [23-24] a mis au point un modèle général pour le béton confiné par divers type
de renforcement transversaux. Ce modèle a été largement utilisé dans l'analyse du
renforcement des colonnes en béton (Xiao, Priestley et Seible, 1993,1994). Dans ce modèle,
l'application de la charge peut être soit statique ou, dynamique et appliquée de façon
monotone ou par des cycles de charge. Le renforcement transversal peut être également de
configuration multiples, par exemple, circulaire ou en spirale, et des cerceaux rectangulaires
13
avec ou sans traverses. La résistance à la compression du béton confiné peut être estimée par
l'équation suivante:
 7,94 f ' l 2 f ' l 

f cc  f c 0   1,254  2,254 1   
 f f c0 
 c0 
La relation unique contrainte déformation proposée par Mander et al. est :
f cc xn
 cc 
n 1 xn
Avec :
c
x 
 cc0 r 
Ec f cc
Ec  Esec , Ec  5,000 f c 0 et E sec 
 f   cc
 cc0   c0  1  5 cc 1
  f c0 
Figure I.6: Courbe contrainte-déformation du modèle Mander et al.
I.1.5.7. Modèle Saatcioglu et al. [25]
 cc  f cc 2 x  ( x) 2 (1 2 K )
1
(Partie ascendante)
 cc  f cc  Z  c   cc 0   0.2 f cc (Partie descendante)
14
0.15  f cc
Z
( 85   cc 0 )
I.1.5.8. Modèle Hoshikuma et al. [26]

 1 
 cc  E c  c 1  x n 1  (Partie ascendante)
 n 
 cc  f cc  E des  c   cc   0.5 f cc (Partie descendante)
11 .2  f c20
E des 
 sh  f yh
Ou :
c
x
 cc0
Avec :
εc : déformation du béton confiné
εcc0 : la déformation correspondante a la contrainte fcc du béton confiné
I.1.5.9. Modèle de Fardis et Khalili

Fardis et Khalili (1981) [27] ont étudié le comportement du béton enrobé de fibre de verre. La
résistance à la rupture et la contrainte ultime de béton confiné par fibres de verre sont
augmenté de manière significative avec le nombre de couches. Un modèle de béton confiné a
été proposé (1982)
  
' '   f com t  
f cc  f c 1  4 .1 
  df '  
  c 
 
 E com t 
 cc  0.002  0.01 
 df ' 
 c 
Où Ecom est le module d'élasticité du confinement dans la direction transversale.

Ce modèle de contrainte-déformation a été proposé sur la base du modèle de Richart et al.
(1928,1929).
15
I.1.5.10. Modèle de Karbhari

Karbhari et al. [28] ont suggéré l'équation suivante pour prédire la contrainte ultime et la
déformation ultime des cylindres composite en béton enveloppés sous une compression
uniforme.
 0 . 87 
   
' '   2 tf com  
f cc  f c 1  2 . 1 
 
  df '  
 c 
 
 
 2tf com 
 cc  0 .002  0 .01 
 df ' 
 c 
Cette loi contrainte-déformation a également été proposée sur la base du modèle proposé par
Richart et al. (1928,1929). Ce modèle ne comprend pas le comportement contrainte-
déformation bilinéaire du béton confiné par FRP.
I.1.5.11. Modèle de Mirmiran

Mirmiran et al. [29] ont développé un modèle pour le béton confiné avec des fibres
composites. Dans ce modèle, la relation à quatre paramètres de Richard et Abbott (1975) a été
utilisée:
 
E1  E 2   E 2 .
n 1 / n


1 
E1  E 2  
 0 
 
Cependant, dans ce modèle E1 et E2 ont été déterminés directement par des résultats
expérimentaux. Ils n’ont pas de signification physique claire.
E t
j j
E 2  245 .61 f c'0 .2  1 .3456 MPa 
D
  0.872 f  0.371 f  6.258 MPa 

o c' r
16
I.1.5.12. Modèle de Xiao et Wu

Xiao et Wu [30] ont étudié le comportement du béton confiné avec matériaux composites
avancés. Ils ont dérivé les expressions de E1 et E2 puis ils ont établi leur sens physique dans
le modèle proposé.
 1 .4 
  f '2  

'  f '  kf  f '   4 .1  0 .45 c  f '  '
f cc c 1 c     r   ro cz
C
 j   l
   
 
0 .9
 f' 
'   E  kC
  10  co  E  kC  ' 2 j
Cj  1 j
 
I.1.5.13. Modèle de Harmon

Harmon et al. (1998) [31] ont développé un modèle pour le béton confiné avec FRP sur la
base du concept de propagation de la fissure dans le béton.
Pour prédire la déformation radiale, Harmon a proposé le concept de frottement interne. Deux
modèles ont été proposés par Harmon et al (1998) basés sur le frottement interne dans le
béton confiné. Ce modèle suppose que la déformation du béton est composée de trois parties.
   elastic   crack   void
Où elastic, crack et void sont la déformation élastique, la déformation due aux fissures, et la
déformation du au vide respectivement.
I.1.5.14. Modèle Fujikake

Sur la base du modèle de comportement orthotrope de Darwin, Fujikake et al (2004) [32] ont
proposé un modèle analytique pour prédire le comportement du béton confiné avec FRP sous
charge de compression axiale.
Un modèle constitutif pour le béton brut a été formulé à base des résultats expérimentaux de
compression triaxiale du béton.
La jaquette FRP a été supposée comme un matériau élastique linéaire. Le modèle de
comportement orthotrope est décrit dans la direction des contraintes principales, lorsque les
axes orthotropes coïncident avec les axes principaux:
17
 d 

 E 1  
1 1 1 23 32
  
E   
1 12 13 32
 
E   
1 13 12 23
 d1 
 d    E 
2    2 21
  
23 31
 
E 1  
2 13 31
 
E   
2 23 21 13
 d 2 

 d 
 3
E 
 3 31
  
32 21
 
E   
3 32 31 12
 
E 1  
3 12 21
  d 
  3 
Ou :
  1   13 31   23 32   12 21   21 13 31   12 23 31
Dans laquelle :
di (i = 1,2,3) = incrément de contrainte principale dans la direction i;
dei (i = 1,2,3)= Incrément de déformation principale dans la direction i;
Ei (i = 1,2,3) = module sécant total de élasticité dans la direction i d'orthotropie ;
νij (i, j = 1,2,3) = coefficient d'anisotropie transversale de déformation dans la direction i
causée par la contrainte dans la direction j. le modèle peut être écrit sous une forme simple tel
que:
d i  D d i 
Où [D] = matrice de rigidité.
La relation équivalente de contrainte-déformation uni-axiale est donnée comme suit :
AX  ( B  1) X 2
i  i i  i  1, 2, 3
p
1  ( A  2) X  BX 2
i i
Où
- A, B = constantes; Xi = déformations uni-axiales équivalentes normalisées;
- p = contrainte ultime correspondante à l’état actuel.
- Les constantes A et B et Xi sont définis par les équations suivantes:
E 0 p
A
p
1 ( A  1) 2   A 2 (1   )   (  1)  
B  
B  1-A
 2  (1   ) 
18
Avec :

  E 0 / E c ,   0 .45 f c' /  p , X i  ui
p
p
B  p
 p  R
ui
m
 
 
X i   ui 
p 
 
Ou :
2
 f' 
 
m  1 .04  2 .0 c 
 100 
 
I.1.5.15. Modèle de Beque

Beque et al. (2003) [33] ont proposé un modèle basé sur le modèle Octaédral de Gerstle
(1981) avec quelques modifications. Les modules de cisaillement et élastiques sont pris
comme fonctions scalaires des invariants. Ainsi, l'état des contraintes peut être défini par les
deux contraintes Octaédrales, la contrainte normale Octaédrale et 0 et la contrainte de
cisaillement Octaédriale 0:
 m  K s  kk
S ij  2G s eij
Avec :
m =  kk / 3 = 0 est la contrainte normale moyenne et Ks et Gs sont appelés le module
sécant module de cisaillement sécant, respectivement.les expressions de Es et νs peuvent
être obtenus par Ks et Gs
 oct E
s
Ks  
3 3(1  2 )
oct s
 E
G s  oct  s
 oct 2(1   s )
19
Ce modèle dépend essentiellement du module Ks et du module de cisaillement Gs.

Modifications faite par rapport aux modèles originaux de Gerstle.
I.1.6. Discussion
La plupart des modèles proposés sont composés de deux phases (partie ascendante et partie
descendante) excepté certains modèles [20-24], [29] qui proposent une seule relation.
La majorité de ces modèles sont développées sur la base des données expérimentales et
contiennent des paramètres donnés par des essais expérimentaux.
Les différentes confrontations effectuées par ces auteurs montrent clairement :
- Deux relations (ascendante puis descendante) ne suffisent pas à elles seules de représenter
correctement le comportement d’une section confinée.
- La relation unique proposée par Mander et al. [23-24] est aussi limitative, l’expression
mathématique n’est pas assez souple pour tenir compte de l’exactitude des paramètres.
- Le confinement est introduit par un coefficient de confinement effectif, le pourcentage

d’armatures transversales et la limite élastique de l’acier dans le cas de confinement par
armature transversales et des propriétés de la jaquette dans le cas de confinement externe.
- Le domaine d’application de ces modèles se limite à des pourcentages de confinement allant

de 20 à 25 ‰, on observe en effet un écart très important (dans la branche poste pic) entre
l’expérimental et le calcul pour des pourcentages inférieurs à 20 ‰.
- Pour tous ces modèles, la partie ascendante est plus au moins représentative, par contre, la
partie descendante après le pic de déformation, s’éloigne des résultats expérimentaux.
I.1.7. Modèle proposé par Bouafia.Y, Iddir. A et Kachi. M.S

Le modèle de confinement proposé par Bouafia,Y., Iddir, A. et Kachi, M.S. [34-37] est
développé suite à une analyse rigoureuse des modèles analytiques existants et une exploitation
des données expérimentales. La courbe définissant la relation contrainte –déformation du
béton confiné est définit par les trois relations suivantes :
La première partie est une courbe ascendante non linéaire, la loi de Sargin [38] a été modifiée
pour tenir compte du confinement :
20
σcc
f cc
f cc
0.65 f cc
Es
Ebc0
 cc0
 65  ccu

Figure I.7: Loi de comportement réel du béton en compression proposéec par
Bouafia, Y. Iddir,A. et Kachi, M.S. [34]
La première partie est une courbe ascendante non linéaire est régie par la loi de Sargin
modifiée :
_2
 
_
kc   c  k  1   c
'
 cc  f cc  c
Pour : 0 ≤ εc ≤ εcc0
_ _2
1  k c  2    c  k   c'
c
Avec :
 c
c 
 cc 0
  f 
 cc 0   c 0  1  5 cc  1
  f c0 
Eb 0   cc 0
kc 
f cc
E bc 0  11000 3 f cc
k c  kc  1
'
21
Suivie d’une partie descendante linéaire après le pic (poste pic), la relation est décrite par :
 cc  f cc  E s  c   cc 0  Pour : εcc0 < εc ≤ ε65
Avec :
6  f c20
Es 
0.35  f cc
 65    cc
Es
Et la troisième partie est un palier linéaire, la relation contrainte –déformation est donnée par :
 cc  0,65 f cc Pour : Pour : ε65 ≤ εc ≤ εccu
fl f
 ccu  0.4 x   cu  0.0035  0.4 x l
f c0 f c0
L’avantage de ce modèle est la proposition de trois relations qui dépendent de plusieurs

paramètres et permet ainsi de mieux approcher le comportement post- pic et de tenir compte
d’une plus large bande de pourcentage d’aciers.
I-2. Deuxième partie : Méthodes de calcul non linéaire
I.2.1. Introduction
Plusieurs logiciels de calcul ont été développés pour étudier le comportement des
structures dans le domaine non linéaire. Certains utilisent la formulation par éléments finis,
d’autres utilisent des approches plus simples. Cependant, dans tous les cas, le domaine
d’étude reste la flexion plane.
22
I.2.2. Proposition de VIRLOGEUX et M’RAD

Plusieurs publications de M. VIRLOGEUX [39-46] font état d’analyses théoriques sur
le comportement des poutres à précontrainte extérieure. Certaines de celles-ci ont servi de
base à l’étude numérique, développée par M’RAD [47], sur le comportement des poutres à
précontrainte partiellement ou totalement extérieure au béton, de sections droites
symétriques par rapport au plan moyen de la structure et chargées dans leur plan en
admettant les hypothèses suivantes :
- conservation de la planéité des sections,
- principe de fissuration : répartie pour les poutres monolithes, et concentrée pour

les poutres à voussoirs,
- absence de glissement entre les différents matériaux composant la section.
Le principe du calcul utilisé est celui pas à pas, du fait du caractère non linéaire du
problème. Les valeurs initiales sont celle obtenues à l’état à vide, les étapes du calcul
itératif sont les suivantes :
- calcul de l’allongement initial des différents câbles de précontrainte,
- calcul des déformations (g, z), sous l’effet des charges extérieures, de la section,
sans considérer la surtension des câbles,
- première évaluation de la surtension des câbles,
- recalcule des déformations (g,z) de la section jusqu’à l’équilibre de la section.
Dans le cas des poutres à voussoirs, le processus est modifié de façon à tenir compte
de la concentration des déformations dans les joints.
La formulation des efforts intérieurs est donnée par :

n
Nint ( g  z ) =   ( b )dydz +   ( sj ) A sj
A j=1
n
Mint ( g  z ) =   ( b ) ydydz +   ( sj ) y jA sj
A j=1
23
Avec : b = g + z y
La résolution du problème consiste à déterminer l’accroissement des déformations en

fonction des efforts extérieurs (Next, Mext), qui sont les sollicitations dues aux charges
extérieures appliquées.
En linéarisant les fonctions des efforts intérieurs au voisinage du point (g , z), dans
l’espace vectoriel des déformations et en supposant que les efforts au point (g ,z) sont
connus. Les efforts au point (g i+1 z i+1) doivent être égaux aux efforts extérieurs (Next,
Mext), l’estimation d’ordre i+1 des déformations est alors donnée par :
-1
  Nint  Nint 
 g i +1   g i     z   Nint - Next 
  =   +  g

 z i +1   z i    Mint  Mint  Mint - Mext 
  g  z 

Le calcul des tensions des câbles se fait par itérations successives à partir des
estimations successives des déformations des sections.
Le modèle de calcul est appliqué aux poutres isostatiques à précontrainte extérieure.

Pour celle-ci les cas suivants sont étudiés :
- câbles extérieur maintenus transversalement sans frottement,
- câble extérieur libre,
- câble dévié.
Ce modèle tient compte uniquement du cas du glissement isolé des câbles sur
déviateur, il utilise alors un processus itératif correcteur pour déterminer la surtension des
câbles extérieurs, ce qui peut poser certains problèmes d’ordre numérique (problème de
convergence). Pour palier à cette difficulté M’RAD [47] a proposé de ralentir le calcul
itératif, de telle façon à n’ajouter qu’une partie de la surtension calculée à la tension du
câble obtenue dans l’itération précédente.
En ce qui concerne les poutres à voussoir, le modèle de fissuration concentré est

appliqué, car la poutre se déforme par l’ouverture des joints et non par une fissuration
régulièrement répartie.
24
Pour trouver une relation qui lie l’ouverture de la fissure à la déformation de la

structure dans la zone du joint ouvert, le modèle introduit des sections de calcul
intermédiaires de part et d’autre du joint fissuré, définies par l’angle de diffusion de la
compression.
La variation de la déformation est alors linéaire entres les sections intermédiaires et

celle du joint.
La surtension est corrigée par un terme qui tient compte de la valeur moyenne des
tensions dans les sections intermédiaires et celle du joint.
I.2.3. Proposition de MULLER J [48]

Le logiciel DFLECT a été développé par le bureau Jean Muller International Inc [48],
il permet d’analyser le comportement des poutres isostatiques et hyperstatiques à
précontrainte extérieure. Dans le cas des ponts à voussoirs, l’accroissement des charges au-
delà de celle de service, ne produit de fissures en dehors de l’ouverture des joints : cela est
vrai dans le cas de joints secs, mais il ne l’est pas pour les joints à résine chimique.
Le programme tient compte du changement de rigidité d’un segment à un autre en

utilisant la relation moment courbure, construite par la méthode des éléments finis pour
chaque segment. L’ensemble de l’analyse est basé sur les hypothèses suivantes :
- la distribution de la contrainte est élastique linéaire dans chaque segment,
- conservation de la planéité de la section du joint,
- le béton en traction entre deux joints successif reste intègre (pas de fissuration),
- le comportement des aciers de précontrainte est élastique non linéaire.
Le programme tient compte également du tracé, des déviateurs et de l’adhérence des

câbles aux déviateurs.
La relation moment-courbure pour les segments est obtenue par la méthode des
éléments finis en supposant que le comportement des matériaux est élastique linéaire. Ce ci
est dû au fait qu’il n y a pas de contrainte limite de fissuration du béton en compression
[49]. La relation moment-courbure est élargie même aux très grandes valeurs. En effet, à la
rupture la courbure de la section est de l’ordre de 450 fois plus grande que celle
correspondant à la même section complètement comprimée. Cela peut s’expliquer par le
25
fait que le programme DEFLECT ne tient pas compte de la plastification du béton à la

rupture.
La validation du logiciel DEFLECT est faite par la comparaison des résultats obtenus
des essais du CEBTP et du CTL. Pour le premier cas, un essai du CEBTP est modélisé
deux fois avec les mêmes caractéristiques géométriques, les mêmes matériaux et le même
tracé des câbles. Dans le premier cas les câbles sont considérés comme extérieurs au béton
et dans le second ils sont considérés comme intérieurs au béton et adhérant sur toute la
longueur de la poutre.
La différence entre les trois cas (CEBTP et les deux modélisations) est insignifiante
pour tous les stades de chargement jusqu’à la charge ultime, figure (I.9)
Figure I. 8 : courbe Charge – flèche obtenue par DEFLECT sur l’essai du CEBTP [49]
Cependant on peut noter que le comportement de la poutre pour les deux cas de
modélisation (précontrainte intérieure et précontrainte extérieure) présente une très petite
différence pour tous les stades de chargement, alors que les résultats expérimentaux
obtenus au CEBTP montrent, cependant, que les rotations sont concentrées dans les joints
où se produit la rupture par écrasement du béton dans le cas des poutres à précontrainte
extérieure, de plus DEFLECT ne tient également pas compte de la plastification du béton à
la rupture, mais il admet un comportement élastique linéaire du béton sans limitation de la
contrainte ou de la déformation.
26
I.2.4. Proposition de Martins

MARTINS [50] a développé un logiciel de calcul permettant d’analyser le
comportement des poutres isostatiques, à voussoirs ou coulées sur place, en béton armé ou
précontraint par des câbles extérieurs au béton ou par un câblage mixte, soumises à un
chargement symétrique. Le tracé du câblage extérieur est limité au cas d’un tracé
symétrique.
Dans son travail Martins utilise l’approche classique de l’élasticité non linéaire pour
résoudre le cas des poutres coulées en place. Pour les poutres à voussoirs, il utilise, pour
tenir compte des voussoirs, l'approche développé par GURIANI [51] pour la distribution
de la déformation du béton dans un élément fissuré.
Ainsi, pour le cas des poutres coulées en place, Martins admet les hypothèses
suivantes :
- les déplacements sont petits et donc la relation moment courbure est donnée par
d 2 v(x) M (x)
=
dx 2
E I(x)
- les sections restent planes après déformation, l’hypothèse de BERNOULLI-
NAVIER reste valable dans tout le domaine de travail,
- l’adhérence acier béton est parfaite jusqu’à rupture,
- les lois de comportement des matériaux sont celle de la sollicitation unidirectionnelle,
- le chargement est quasi statique, monotone et croissant,
- les effets différés sont négligés,
- on tient compte de la contribution du béton tendu entre les fissures.
Le système d’équations à résoudre dans la recherche de l’équilibre d’une section

transversale sollicitée par les efforts (Next,Mext) est le suivant :
-1
  Nint  Nint 
 g i +1   g i     z   Nint - Next 
 =   +   Mint 
g
  Mint  Mint - Mext 
 z i +1   z i  
  g  z 

Nint et Mint étant les efforts résistants correspondant à un état de déformation (g , z)
donné par :
27
n
Nint ( g  z ) =   ( b )dydz +   ( sj ) A sj
A j=1
n
Mint ( g  z ) =   ( b ) ydydz +   ( sj ) y jA sj
A j=1
Avec b = g + z y
sj = g + z yj
b étant la déformation dans le béton au niveau du point d’ordonnée y et sj la déformation

de l’acier d’ordonnée yj.
La section est discrétisée en m fibres de béton et le système d’équations est résolu par
la méthode de NEWTON-RAPHSON.
Dans le cas des poutres a voussoir cette approche est modifiée pour tenir compte des
voussoirs. Une formulation plus générale est alors développée pour étudier les poutres a
voussoirs à précontrainte extérieure en tenant compte de la dégradation de l’adhérence
acier-béton après ouverture des fissures (joints). Cette formulation peut être appliquée aux
poutres monolithes présentant un schéma de fissuration non régulier. Elle est basée sur la
méthode moment-courbure ainsi que sur le modèle de GURIANI [51] pour la distribution
de la déformation du béton dans un élément fissuré.
Cette approche admet qu’une section fissurée (ou de joint) se déforme après ouverture
de la fissure (joint) en tournant autour de l’axe de déplacement nul (uc =0 ) par rapport à la
section médiane entre deux fissures (joints) successives. Elle est basée sur les hypothèses
suivantes :
- Une section fissurée se déforme après ouverture de la fissure en tournant autour de

l’axe de déplacement nul (uc =0 ) par rapport à la section médiane entre deux fissures
consécutives, figure(I.10).
28
Figure I.9: L’équilibre d’un voussoir.

- la section reste plane, au moins dans la partie du béton comprimé après
fissuration,
- le moment fléchissant est constant le long d’un voussoir,
- les effets de l’effort tranchant sont négligés,
- l’aire des diagrammes des déformations de compression dans les sections médiane
et de fissure ou joint sont égales,
- la déformation du béton sur la fibre la plus comprimée évolue selon une loi
parabolique entre la section de fissure ou joint et celle à mi-voussoir,
- la distribution des déformations de compression dans la section du joint est

linéaire ainsi que celle dans la section à mi-voussoir.
Une modification a été introduite dans la formulation classique de l’élasticité non

linéaire, pour tenir compte, du comportement d’un voussoir. Cette modification porte sur la

relation liant les variables classiques (g , z) aux variables  = et f(y) de la
Lv
29
formulation du comportement d’un voussoir, où  est la rotation de la section j du joint par

rapport à la section médiane du voussoir, rapportée à l’unité de longueur du voussoir Lv, et
où f(y) une fonction permettant d’obtenir la déformation des fibres de béton dans la section
du joint (fissure).
La détermination de l’ouverture des fissures (joint) n’est pas simple, les valeurs
calculées par MARTINS sont des valeurs moyennes probables.
En effet, le modèle tient compte de la dégradation de l’adhérence entre l’acier et le

béton, par une loi polygonale d’adhérence - glissement des câbles, proposée par Martins
comme résultat de l’analyse d’essais rapportés dans la littérature technique.
I.2.5. Proposition de Kachi

La méthode de calcul présentée par Kachi [52] a pour objectif l’introduction de la
précontrainte extérieure en tenant compte des différentes dispositions des câbles au droit des
déviateurs (câbles adhérant et n’adhérant pas aux déviateurs), ainsi que l’effet de l’effort
tranchant dans le domaine élastique non linéaire dans un logiciel de calcul des poutres
hyperstatiques développé à partir des travaux de Nait Rabah[53] et de Garcia [54] au CEBTP.
L’effet de l’effort tranchant dans le domaine non linéaire est introduit dans le modèle
de calcul par une estimation du module de cisaillement réel d’une section de poutre par la
méthode présentée et par son introduction dans la matrice de rigidité corde de la section.
La méthode de calcul a pour objectif la simulation du comportement jusqu’à rupture
de poutres hyperstatiques en béton armé, béton précontraint (par des câbles intérieurs ou
extérieurs au béton). Elle tient compte de trois types de non linéarités :
- celles résultant des lois de comportement intrinsèques des matériaux. Ces non-
linéarités conduisent à des relations non-linéaires entres les sollicitations et
déformations ;
- les non-linéarités géométriques dues au comportement d’ensemble de la structure.
Il est indispensable de tenir compte de ce type de non-linéarités dans l’analyse
structurale des systèmes où les effets du second ordre ou bien la stabilité des
éléments ou de l’ensemble du système sont très importants ;
- celles résultant du glissement des câbles extérieurs sur les déviateurs provoquant la
modification dans les surtensions et des déformations dans l’ensemble de la structure.
30
I.2.5.1 Domaine d’application et hypothèses générales
La structure étudiée par Kachi, MS. [52], Kachi, MS., Bouafia, Y., Fouré, B., Muller,
P. [55-60] est modélisée en une ossature constituée par deux types d’éléments : un élément
poutre et un élément câble lié aux éléments poutre par des éléments rigides, qui font dans leur
ensemble l’objet d’une analyse non linéaire.
La procédure de calcul est basée sur la méthode des déplacements pour résoudre
l’équation d’équilibre de la structure qui lie les accroissements des déplacements aux
accroissements de sollicitation de charges extérieures et à l’action éventuelle de la
précontrainte, par des câbles intérieurs à la mise en tension.
La structure étudiée est considérée comme soumise à un système de charges
instantané. Les charges appliquées sont supposées concentrées aux nœuds de la structure.
Le système de chargement considéré comporte :
- des charges Q01 et Q0é fixes, appliquées avant et après précontrainte
- des charges  Q v variables, d’intensité proportionnelle à un paramètre.
L’analyse de la structure est effectuée pas à pas pour des valeurs croissantes de 
La précontrainte intérieure au béton est étudiée par une méthode basée sur la formulation de
Nait Rabah [53].
L’analyse de l’état des déformations d’un élément est effectuée dans un repère
intrinsèque lié à la position déformée de cet élément. Dans cette analyse les hypothèses faites
sont les suivantes :
- l’analyse est effectuée en élasticité non linéaire ;
- les déformations sont petites ;
- les déplacements sont petits ;
- la longueur des éléments est faible et donc la non linéarité géométrique attachée à
la déformation de l’élément est négligeable ;
- les déformations des nœuds situés à la jonction de plusieurs éléments sont
négligeables ;
- le glissement des câbles extérieurs est nul.
Les caractéristiques de déformabilité d’un élément sont obtenues en étudiant le
comportement mécanique d’une ou plusieurs sections droites de cet élément.
31
L’analyse non linéaire concerne aussi bien les sollicitations dues à l’effort normal et au
moment fléchissant que celle due à l’effort tranchant. Les effets différés et à longue durée
ne sont pas pris en compte.
I.2.5.2. Etude de l’équilibre d’une poutre en élasticité non linéaire

En élasticité non linéaire, l’étude de l’équilibre d’une poutre en béton armé et précontraint
passe par la résolution d’un système d’équations de la forme :
[P] = [K]. [U]
Où [P] est le vecteur représentant l’accroissement du chargement appliqué à la poutre, [U]

est le vecteur représentant l’accroissement des déplacements des nœuds de la poutre et [K]
représente la matrice de rigidité globale de la poutre, qui est construite à partir des matrices de
rigidité des sections [Ks].
L’expression de cette matrice adoptée par les auteurs Kachi, MS. [52], Kachi, MS.,
Bouafia, Y., Fouré, B., Muller, P. [55-60] pour tenir compte de l’effort tranchant. En
élasticité non linéaire, pour une poutre en béton armé, est :
 N N 0 
  u  w 
[Ks] =  M M 
  u  w 0 
 V 
 0 0
  moy 
 
Où N est l’accroissement de l’effort normal de la section, M est l’accroissement du

moment fléchissant de la section, u est l’accroissement de la déformation au niveau centre
de gravité de la section et w est l’accroissement de la rotation différentielle (courbure) de
la section.
Ces déformations répondent à l’hypothèse de planéité des sections (Navier-Bernoulli) ; la

déformation longitudinale x est une fonction linéaire de l’abscisse transversale y :
x (y) = u + w. y
32
V est l’accroissement de l’effort tranchant de la section et moy est l’accroissement de la

distorsion moyenne de la section.
L’expression adoptée ci-dessus pour [Ks] est déjà une simplification. En effet, les termes
de couplage entre N et w, et entre M et u, sont incontournables dans l’analyse de
flexion composée, mais les auteurs ont supposé que les termes de couplage entre N et moy
(ou V et u), et entre M et moy (ou V et w), étaient négligeables. Le couplage entre
les déformations dues à V et celles dues à (N , M) évoqué dans l’introduction va se trouver
dans la dépendance de moy avec les contraintes et déformations dues à N et M.
Pour les poutres en béton précontraint, la matrice [Ks] est complétée par des termes dus à
l’effet de la précontrainte, dans l’équilibre de flexion composée et dans l’équilibre d’une
partie de l’effort tranchant par la composante transversale de l’effort des câbles inclinés –
Nait-Rabah (1990) [53].
I.3.Conclusion
Le développement des calculs aux états limites impose l’élaboration de méthodes de calcul
tenant compte du comportement non linéaire des matériaux .Le développement parallèle de
l’outil informatique permet de mettre au point des méthodes numériques assez complexes
pour mieux caractériser le comportement des sections particulièrement en béton confiné.
Cette synthèse sur les modèles de confinement et le calcul non linéaire des structures nous
permet de mieux appréhender son application aux cas des structures de sections circulaires.
Vu la complexité des modèles existants et la difficulté de maîtriser certain de leur paramètres,
nous sommes amenés à proposer un modèle simple dont les paramètres sont facilement
maîtrisable.
Le confinement est introduit par un coefficient de confinement effectif, le pourcentage
d’armatures transversales et la limite élastique de l’acier ainsi que l’épaisseur des tubes en
acier dans le cas des poteaux mixtes.
33
Chapitre II Modélisation avec prise en compte du confinement
CHAPITRE II
MODELISATION AVEC PRISE EN

COMPTE DU CONFINEMENT
II.1 Introduction
La méthode de calcul qu’on présente a pour objectif l’introduction du confinement en
tenant compte des différentes dispositions des armatures transversales ou des tubes, ainsi que
l’effet de l’effort tranchant dans le domaine élastique non linéaire dans un logiciel de calcul
des poutres hyperstatiques développé à l’université de Tizi-Ouzou Kachi [52], Kachi, Fouré,
Bouafia, Muller [52-60]. Le confinement est introduit par une loi de comportement proposée
par Bouafia, Y. et Iddir, A. et Kachi, M.S. [34-37] qui permet de mieux simuler le
comportement des sections en béton armé confiné.
Certains problèmes spécifiques, tels que les effets du second ordre, sont pris en
compte dans la formulation du système d’équations d’équilibre de la structure.
L’effet de l’effort tranchant dans le domaine non linéaire est introduit dans le modèle
de calcul par une estimation du module de cisaillement réel d’une section de poutre par la
méthode développée par Kachi et al. [52], [55-60] et par son introduction dans la matrice de
rigidité corde de la section.
La méthode de calcul a pour objectif la simulation du comportement jusqu’à rupture
de poutres hyperstatiques en béton armé avec prise en compte du confinement, Elle tient
compte de deux types de non linéarités :
- celles résultant des lois de comportement intrinsèques des matériaux. Ces non-
linéarités conduisent à des relations non-linéaires entres les sollicitations et
déformations.
34
- les non-linéarités géométriques dues au comportement d’ensemble de la structure.

Il est indispensable de tenir compte de ce type de non-linéarités dans l’analyse
structurale des systèmes où les effets du second ordre ou bien la stabilité des
éléments ou de l’ensemble du système sont très importants.
II.2 Domaine d’application et hypothèses générales
La structure étudiée est modélisée en une ossature constituée par deux types
d’éléments : un élément poutre et un élément câble lié aux éléments poutre par des éléments
rigides, qui font dans leur ensemble l’objet d’une analyse non linéaire.
La procédure de calcul est basée sur la méthode des déplacements pour résoudre
l’équation d’équilibre de la structure qui lie les accroissements des déplacements aux
accroissements de sollicitation de charges extérieures et à l’action éventuelle de la
précontrainte, par des câbles intérieurs à la mise en tension.
La structure étudiée est considérée comme soumise à un système de charges
instantané. Les charges appliquées sont supposées concentrées aux nœuds de la structure.
Le système de chargement considéré comporte :
- des charges Q01 et Q0é fixes, appliquées avant et après précontrainte
- des charges  Q v variables, d’intensité proportionnelle à un paramètre.
L’analyse de la structure est effectuée pas à pas pour des valeurs croissantes de 
La précontrainte intérieure au béton est étudiée par une méthode basée sur la formulation de
Nait Rabah [53].
L’analyse de l’état des déformations d’un élément est effectuée dans un repère
intrinsèque lié à la position déformée de cet élément. Dans cette analyse les hypothèses faites
sont les suivantes :
- l’analyse est effectuée en élasticité non linéaire ;
- les déformations sont petites ;
- les déplacements sont petits ;
- la longueur des éléments est faible et donc la non linéarité géométrique attachée à
la déformation de l’élément est négligeable ;
- les déformations des nœuds situés à la jonction de plusieurs éléments sont
négligeables ;
35
Les caractéristiques de déformabilité d’un élément sont obtenues en étudiant le

comportement mécanique d’une ou plusieurs sections droites de cet élément.
L’analyse non linéaire concerne aussi bien les sollicitations dues à l’effort normal et au
moment fléchissant que celle due à l’effort tranchant. Les effets différés et à longue durée
ne sont pas pris en compte.
Le calcul des déplacements aux extrémités de l’élément dans le repère intrinsèque est
effectué par intégration des déformations des sections de calcul.
La longueur des éléments et le nombre de sections droites de calcul ont une influence
sur les déplacements calculés. Grelat [61] a montré que, pour les éléments de faible
longueur, la prise en compte de trois sections de calcul est en générale suffisante.
L’approche utilisée pour le calcul des déplacements est basée sur la formulation de Nait
Rabah [53] conduisant à la formation de la matrice de souplesse reliant les accroissements
des déplacements aux accroissements d’efforts aux extrémités de l’élément. Les effets du
second ordre dus aux déplacements des nœuds sont introduits essentiellement par la
transformation non linéaire des déplacements des extrémités de l’élément dans le passage
du repère intrinsèque à un repère fixe.
II.3. Comportement des matériaux

Les lois de comportement uni axiales sont usuellement utilisées dans l’étude du
comportement des matériaux. Cette hypothèse est acceptable dans le cas d’un chargement
instantané croissant proportionnellement à un paramètre.
II.3.1. Béton comprimé

Pour décrire le comportement du béton en compression, on admet un comportement
élastique non linéaire qui sera décrit par la loi de Sargin [38] (figure II.1).
La contrainte est donnée par :
2
k b   ( k b,  1)
  fcj 2 (2.1)
1  ( k b  2)  k b' 

Où  =
0
- Le paramètre kb ajuste la branche ascendante de la loi et est donné par:
E b0 . 0
kb =
f cj
. fcj : résistance à la compression du béton à l’âge j,
36
. 0 : déformation de pic correspondant à fcj,

. Eb0 : module élastique du béton à l’origine.
Le paramètre K’b permet d’ajuster la forme de la branche décroissante de la courbe.
- Avec K’b = 0, on obtient des courbes avec une branche qui descend brutalement, ce qui
correspond à un béton ayant un comportement élastique fragile.
- Avec K’b = 1, on obtient au contraire des courbes assez plates qui correspondent à un béton
ayant un comportement ductile.
Pour un béton normal, on prend généralement k’b = kb – 1
σc
k b  0
kb kb1
f c0 k  1
b
Eb 0
0  cu 
Figure II.10 : Comportement du béton en compression loi de Sargin [38]
II.3.2. Comportement du béton en traction

Généralement, dans la pratique (règlements), la résistance à la traction du béton ainsi
que la contribution du béton tendu entre deux fissures sont négligées.
La prise en compte de la contribution du béton tendu entre les fissures, permet une
représentation convenable de la déformabilité réelle d’un élément fissuré, ainsi qu’une
estimation meilleure de la redistribution des efforts dans les structures hyperstatiques et dans
les cas sensibles aux effets du second ordre. Elle permet ainsi de mieux estimer les
déplacements.
Plusieurs modèles sont proposés dans la littérature. Dans le cadre de la présente étude
nous avons utilisé la loi de Grelat [61]. Elle attribue au béton tendu fissuré un diagramme
fictif triangulaire à partir de l’axe neutre (figure II.2).
37
bc bc
zone comprimé
h
n o o n
zone tendue
b bt bt
Section déformation contrainte
Figure II.11: Comportement du béton avant fissuration.

Les fibres du béton tendu sont affectées d’un module de déformation Ebt qui est défini à partir
de la contrainte et de la déformation instantanée du bord tendu.
bt
Ebt = (2.2)
 bt
Au - delà de la fissuration en traction on tient compte d’une participation du béton
tendu situé entre deux fissures successives. La contrainte ne s’annule pas brusquement, mais,
décroît selon une loi parabolique (figure II.3). On étudie ainsi le comportement moyen d’une
zone de poutre.
ftj
Eb0
ft rt 
Figure II.12 : Contribution du béton tendu après fissuration A.Grelat [61].
ftj : résistance du béton à la traction,

ft : déformation de traction correspondant à ftj,
rt : déformation correspondant à la plastification de l’acier le plus tendu,
Ebo : module d’élasticité longitudinal du béton.
 bt = Ebo . bt si bt   ft (2.3)
 bt  rt et bt   ft fissuration avec participation du béton tendu.
38
(   rt ) 2
bt = - fft (2.4)
( rt   ft ) 2
 bt   rt fissuration sans participation du béton tendu

bt = 0
II.3.3. Les aciers

La loi contrainte-déformation est supposée identique en traction et en compression
sous chargement monotone, elle dépend de la nature de l’acier.
II.3.3.1 Comportement des aciers passifs

On distingue les aciers naturels et les aciers écrouis.
A) Aciers naturels: Ils sont caractérisés par une loi élasto-plastique parfaite (figure II.4)

e
-u - e
e u 
- e
Figure II.13: Comportement des aciers naturels.
 = Ea  pour   e
 = e pour e    u
=0 pour   u (2.5)
Du coté de la traction les relations se déduisent par symétrie.
Les déformations extrêmes sont fixées par le règlement BAEL à 10 o/oo
Ea : Module longitudinal de l’acier,
e : Déformation limite élastique de l’acier,
e : Contrainte limite élastique de l’acier,
u : Déformation ultime de l’acier.
39
b) Aciers écrouis
On adopte la loi préconisée par les règles BAEL 83 qui supposent un comportement
linéaire élastique jusqu’à 70  de la limite élastique (e). Au - delà, une courbe  = f (  )
du cinquième degré (figure II.5).
  e
e
  e
    % u 
Figure II.14: Comportement des aciers écrouis.

Elle se traduit par les relations suivantes :
 = Es.  pour   0.7  e
 
= + 0.823 ( - 0.7 )5 pour 0.7 e    1.1e
Ea e
 = 1.1 e pour 1     u (2.6)
(Avec les relations qui se déduisent en traction par symétrie).
II.3.3.2 Comportement des aciers actifs

Le comportement des aciers de précontrainte est représenté par la loi, du même type
que celle des aciers écrouis du béton armé, recommandée par le règlement BPEL 83.
  a p pour     
p p
p = + 100 (  0.9 ) 5 pour 0.9e  p  1.06e
a fe
p = 1.06 e pour    (2.7)
(Avec les relations qui se déduisent par symétrie du coté de la traction).
40
II.3.3. 3. Comportement des tubes en aciers
La section métallique est habituellement faite à partir de l'acier de construction. Pour

expliquer les propriétés mécaniques de l'acier de construction, l'essai le plus approprié est un
essai de traction sous chargement statique. La relation de contrainte- déformation pour un
spécimen en traction d'un tel acier est montrée sur la figure II.6 [2].
Au commencement, la relation de contrainte- déformation est linéaire avec un module
élastique Ea qui pour la plupart des qualités des aciers est de près de 200 GPa. Si l'effort est
dans la phase élastique, l'acier récupère parfaitement en déchargement. Le comportement
élastique linéaire continue jusqu'à ce que la force limite fy soit atteinte, à une déformation
limite ay.
Autres déformations résultent en phase d'écoulement plastique avec peu ou pas
d'augmentation d'effort jusqu'à ce que l'écrouissage soit atteint à une contrainte oh Alors,
l’effort dans acier augmente jusqu' atteindre la résistance limite en traction fu.
La déformation commence alors à se localiser, ayant pour résultat une grande réduction de
l'aire de la section droite qui cause une diminution de contrainte jusqu'à la rupture.
Cependant, la dernière partie de la relation de contrainte- déformation est légèrement
fallacieuse concernant le comportement du matériau.
Figure II.15: Courbe contrainte-déformation des tubes en aciers
41
II.3.4. Comportement du béton confiné
II.3.4.1.Cas d’une section circulaire
- Résistance à la compression du béton confiné f’cc
Enrobage
Spirale en acier
Béton effectivement confiné

Figure II.16 : section circulaire confinée par spirale en acier.
La résistance à la compression du béton confiné est donnée par [23-24] :
 7,94 f ' l 2 f ' l 

f cc  f c 0   1,254  2, 254 1    (2.8)
 f f c0 
 c 0 
- f c 0 : résistance à la compression du béton non confiné.
- f’l : pression de confinement latérale effective.

f 'l  f l  k e (2.9)
- fl : pression latérale du renforcement transversale supposer uniformément distribuée
sur la surface du béton.
- ke : coefficient de confinement effectif.
Ae
k e (2.10)
Acc
- Ae : section effectivement confinée.
Acc  Ac 1   cc  (2.11)
- ρcc : rapport entre la section des armatures longitudinales et la section du béton.
- Ac : section interne comprise entre les spirales.
42
ds
A A
Enrobage
Section B-B Béton effectivement confiné

Enrobage
s' s
B B
ds - s’/2
Béton effectivement
Confiné
ds
Section A-A
Figure II.17: béton effectivement confiné (cas d’une section circulaire confinée par spirale en
acier).
2
 s' 
Ae   ds  
4 2
2
 2 s' 
Ae  ds 1   (2.12)
4  2ds 
 2
Acc  ds 1   CC  (2.13)
4
Donc :
 s' 
1  
Ae  2ds 
Ke  
ACC 1   CC
 s' 
1  
 2ds 
Ke  (2.14)
1   CC
43
L’équilibre des forces dans une section donnée est :
2 f yh  Asp  f l  s  ds (2.15)
- fyh : contrainte limite des aciers transversaux.

- Asp : section des aciers transversaux.
- fl : contrainte de confinement latérale.
- s : distance entre les extrémités des spirales (voir section A-A).
- s : rapport entre le volume des armatures de confinement transversales et le volume
du béton confiné.
4 Asp 4 Asp
s  ; s  ds 
ds  s s
2 f yh  Asp
fl  (2.16)
s  ds
1
Donc f l   S  f yh (2.17)
2
En remplaçant (2.14) et (2.17) dans (2.9) :
 s' 
1  
1  2ds 
f 'l   S  f yh 
2 1   cc
1
f 'l   S  ke  f yh (2.18)
2
II.3.4.2. Cas d’une section rectangulaire

Pour une section rectangulaire la section effectivement confinée est celle obtenue après
avoir soustraire toutes les paraboles non confiné.
Pour une parabole la section est de (’i) 2/6.
- (’i) : la distance entre deux armatures longitudinales adjacente.
44
bc
w’
x
dc
dc -s’/2
y y

y
Section z-z
z s’ z
s
Enrobage
bc - s’/2 Béton effectivement

Confiné
bc
Section y-y
Figure II.18 : béton effectivement confiné (cas d’une section rectangulaire confinée par
armatures transversales).
Pour une section de n barres longitudinales la section non confiné à soustraire est :
( ' )
2
n
Ai   i (2.19)
i 1 6
 n    
2
Ae   bc  dc   i   1  s '   1  s '  (2.20)

6   2bc   2dc 
 i 1

Acc  bc  dc1   cc  (2.21)

1  n ωi' 
2

  1  s '   1  s ' 
 i 1 6  bc  dc  
  
2bc   2dc 

Ae  
ke  
A cc 1  cc 
45

1  n

 '
i
2
 

  1  s '    1  s ' 
 i 1 6  bc  dc   2bc   2dc 
ke    (2.22)
1  cc 
On peut exprimer :
Asx
x 
s  dc
Asy
y 
s  bc
Asx
f lx  f yn   x  f yh (2.23)
s  dc
Asy
f ly  f yn   y  f yh (2.24)
s  bc
On aura :
f 'lx  ke   x  f yh (2.25)
f 'ly  ke   y  f yh (2.26)
La distribution de la contrainte de confinement latérale se fait selon un angle de 45°, la valeur

de fl est donnée par :
f 'lx  f ' ly
f l 45 
2
II.3.4.3. Cas d’une section circulaire confiné par tube (CFST :Concrete filled steel tubes)
ds
A A
Tube
tj
Section B-B
Figure II.19 : béton effectivement confiné (cas d’une section circulaire confinée par tube en
acier).
- f’l : pression de confinement latérale effective.
46
f 'l  f le  k e (2.27)
- fle : pression latérale du renforcement transversale supposer uniformément distribuée
sur la surface du béton.
- ke : coefficient de confinement effectif.
Ae
k e (2.28)
Acc
- Ae : section effectivement confinée.
2
 s' 
Ae   ds  
4 2
Dans le cas d’un tube l’espacement s’, est égale à zéro (s’=0), donc :
Ae = ds2 /4
- Acc  Ac 1   cc  (2.29)
Ac : section interne comprise entre les spirales
Ac = ds2 /4
- ρcc : rapport entre la section des armatures longitudinales et la section du béton.
ρcc = Asl /Ac
- Asl : section des aciers longitudinaux.
Dans le cas d’un tube la section des aciers longitudinaux, est égale à zéro (Asl =0),
Donc : ρcc=0,
D’où : Acc  Ac ,
Ae
Finalement : k e 1 ,
Ac
Et pression de confinement latérale effective : f ' l  f l e
47
L’équilibre des forces dans une section donnée est :
ds
tj
flu
fjuforces dans une section
Figure II.20 :L’équilibre des fju circulaire confinée par tube en acier).
2 f ju  t j  f le  ds
2 f ju  t j
f le 
ds
- fju : contrainte limite du tube.

- fl : contrainte de confinement latérale.
- tj : épaisseur du tube (voir figure II.11).
- ds : diamètre de la section du béton
Introduisant le pourcentage de confinement j : rapport entre le volume du tube et le volume

du béton confiné.
  ds  t j  l 4 t
j  Donc :   j
  ds 2 j
ds
l
4
1
D’où : f le   j  f ju
2
1 1
Finalement : f 'l   ke  j  f ju   j  f ju (2.30)
2 2
48
II.3.4.4. Modèle
Le modèle de confinement utilisé dans cette étude est celui développé Bouafia, Y. et Iddir, A.
et Kachi, M.S. [34-37].
La courbe définissant la relation contrainte –déformation du béton confiné est définit par trois
relations (figure II.12) ;
σcc
f cc
f cc
0.65 f cc
Es
Ebc0
 cc 0  65  ccu
c
Figure II.21 : Loi de comportement réel du béton en compression
modifiée :
_2
 
_
kc   c  k  1   c
'
 cc  f cc  c
Pour : 0 ≤ εc ≤ εcc0 (2.31)
_ _2
1  k c  2    c  k   c '
c
Avec :
 c
c 
 cc 0
  f 
 cc 0   c 0  1  5 cc  1
  f c0 
Eb 0   cc 0
kc 
f cc
Ebc 0  11000 3 f cc
k c  kc  1
'
49
 cc  f cc  E s  c   cc 0  Pour : εcc0 < εc ≤ ε65 (2.32)
Avec :
6  f c20
Es 
0.35  f cc
 65    cc
Es
 cc  0,65 f cc Pour : Pour : ε65 ≤ εc ≤ εccu (2.33)
Le béton non armé possède une capacité de déformation à la rupture εcu très limitée, de l’ordre
de 3,5‰. Cette valeur est 14 fois plus petite que la déformation εau, de l’armature (50‰). La
déformation à rupture εcu est de l’ordre de 2 fois le raccourcissement élastique maximum εc0
du béton. Soit 20 fois moins que celle de l’armature.
εcu peut être relevé par un facteur de l’ordre de 2 à 4 si le béton est bien confiné par une
armature transversale, circonstance qui augmente aussi sa résistance. L’armature de
confinement empêche l’expansion transversale du béton (effet Poisson) ; elle est soumise à
une pression perpendiculaire à la surface, qui la met en traction (de type « circonférentielle »).
La déformation à la rupture εccu d’un béton confiné d’après: L’Eurocode 2 est donnée :
εccu= εcu + 0,2 fc0 / fl
Où, fl est la contrainte de compression latérale effective due au confinement.
La ductilité du béton est alors de l’ordre de 4 ou plus, mais seulement dans la zone intérieure
au confinement ; la zone extérieure à l’armature transversale se désagrège en cas de
déformation importante. La déformation à rupture εccu, du béton confiné peut être estimée en
fonction de l’armature transversale. La contrainte de confinement est l’ordre 1/10 de fc0.par
conséquence, et afin de faire participer l’acier jusqu'à sa plastification puis rupture autour de
40‰ en terme de déformation, nous proposons la relation suivante :
fl f
 ccu  0.4 x   cu  0.0035  0.4 x l (2.33 Bis)
f c0 f c0
50
II.4. Modélisation des sections
II.4.1. Le béton
Pour mieux approcher le domaine délimité par le contour d’une section transversale de
forme quelconque, on définit la section de béton par une succession de tables trapézoïdales.
Chaque table est définie par les largeurs des fibres extrêmes bj , bj+1 et les ordonnées de ces
fibres yj , yj+1 par rapport à un axe de référence (figure II.13).
bj+1
bj yj+1
yj
z
Figure II.22 : Section de béton découpée en tables trapézoïdales.

Le procédé d’intégration étant numérique, il est nécessaire d’exprimer la largeur du trapèze
d’ordre j à une cote y par une relation compatible avec le calcul numérique.
B(y) = bj +(bj+1 - bj )(y - yj ) / (yj+1 -yj ) (2.34)
II.4.2. Les aciers

La section de chaque armature est concentrée en son centre de gravité. Les armatures
passives sont donc définies comme une succession de lits d’armatures. Chaque lit est défini
par son ordonnée yai et par l’aire totale des armatures situées à ce niveau Aai.
Les armatures actives sont définies par leur excentricité epk et leur aire Apk, .
ns : nombre de lit d’armatures passives,
np : nombre de câbles de précontrainte.
A la mise en tension, les câbles subissent un allongement initial po. L’action de la
précontrainte sur la section est un effort de compression, égal et opposé à l’effort de traction
dans le câble. Cet effort est appliqué suivant la ligne moyenne du câble au point où il traverse
la section. Quand le câble est rendu solidaire de la poutre (ancrage, injection de la gaine) il
51
subit les mêmes variations de déformation que le béton situé au même niveau. Il devient alors
un élément résistant, et possède une pré-déformation par rapport à la section.
II.5. Formulation du principe de calcul
La structure réelle est modélisée en une ossature formée par des éléments poutre et
des éléments câbles liés aux éléments poutres par des éléments rigides, dont les extrémités
constituent l’ensemble des nœuds de la structure.
La structure est supposée soumise à un système de charge instantané. Les charges
exercées sont appliquées aux nœuds de la structure.
Le système de charges considéré comporte des charges d’intensité fixe et des charges
d’intensité variable proportionnellement à un paramètre.
L’analyse de la structure est effectuée pas à pas pour des valeurs croissantes de.
Le problème non linéaire est résolu en utilisant une méthode itérative, basée sur la
méthode des déplacements, dans laquelle on recalcule à chaque étape la matrice de rigidité
tangente reliant les accroissements des déplacements aux accroissements des efforts.
Les ancrages des câbles sont des nœuds fixes (glissement nul) et tous les autres nœuds
des éléments câble sont identifiés et associés aux éléments câbles de part et d’autre de chaque
déviateur où le glissement doit être testé après chaque étape d’équilibre de la structure. La
rotation des nœuds des éléments câbles est libre.
Pour un accroissement de chargement extérieur Q donné, on construit la matrice
rigidité  K de la structure et le vecteur P dû à l'action éventuelle de la précontrainte à la

mise en tension des câbles intérieurs. Il contient aussi la variation de la tension, dans les
éléments câbles extérieurs, générée par un éventuel glissement de ceux-ci sur les déviateurs.
Pour chaque étape de chargement, le problème consiste à chercher l'accroissement de
déplacement des nœuds U solution du système d’équations non linéaire suivant :
Q + P= K . U (2.35)
La matrice de rigidité [K] et le vecteur P sont formés respectivement à partir des
matrices  KX et des vecteurs Px  des différents éléments dans le repère absolu de la
structure, compte tenu de l'équilibre des nœuds de la structure.
52
La matrice de rigidité KX  et le vecteur Px  de chaque élément de la structure sont

calculés en étudiant le comportement mécanique d'un certain nombre de sections transversales
le long de l'élément dans un repère dit " intrinsèque " lié à la position déformée de l'élément.
La non linéarité géométrique liée à la déformation propre de l’élément est supposée
négligeable. Les effets du second ordre dus aux déplacements des nœuds sont introduits dans
le passage du repère intrinsèque au repère local lié à la position initiale de l’élément, sous
l’hypothèse des petits déplacements. Le passage au repère global, lié à la structure, est
effectué au moyen d’une matrice de rotation TG qui ne dépend que de l’inclinaison initiale de
l’élément.
La recherche de l’équilibre de la structure sous un accroissement du chargement
extérieur est effectuée d’abord en admettant que les glissements des câbles extérieurs sur les
déviateurs sont nuls. Une fois l’équilibre de la structure acquis, on teste la condition de non
glissement sur chaque déviateur. Si cette condition est vérifiée, les résultats obtenus en
efforts, déplacements et déformations sont validés, on passe alors à une nouvelle étape de
chargement. Une étape de correction sera alors nécessaire pour tenir compte de la
modification des tensions dans les câbles due à leur glissement et de leur conséquence sur
l’ensemble de la structure. Dans ce cas, le vecteur Q de l’équation (2.35) est nul.
A chaque accroissement du chargement Q donné, le système d’équation (2.35) est

résolu en utilisant la méthode des substitutions successives basée sur les matrices cordes
reliant l’accroissement des déplacements des nœuds de la structure à celui du chargement
extérieur appliqué (figure 2.14).
  
Cette méthode consiste à construire une suite de solutions U 0 U1 ……….. U i ;
 
U r est calculée à partir de Ui -1 en résolvant le système linéaire suivant :
  
Q  Pi -1  K i 1 U

Où K i 1 représente la matrice de rigidité corde de la structure à l’itération r-1 et Pi -1 le
vecteur tenant compte de l’action éventuelle de la précontrainte à la mise en tension des
câbles intérieurs ou de la surtension des câbles extérieurs résultant du glissement éventuel de
ceux-ci sur les déviateurs.
Pour le démarrage de l’algorithme la condition initiale en déplacements est donnée

 
par U 0 = 0 .
53

Q
 
Étape à équilibrer Q + Q

étape stable Q
   
Us U s + U U
Figure II.23 : Représentation de la méthode de substitution utilisant les matrices cordes.

Le critère de convergence des déplacements est analogue à celui utilisé pour la
résolution de l’équilibre d’une section en déformation. Les critères de rupture utilisés sont
ceux classiques en flexion composée et ceux de rupture en cisaillement (Ecrasement des
bielles de béton, rupture des aciers transversaux).
II.5.1. Séquence de calcul

L'analyse de la structure est effectuée  pas à pas  pour des valeurs croissantes d'un
paramètre. On suppose connu l'état complet de la structure (déplacements, sollicitations,
déformations et glissements des câbles extérieurs sur les déviateurs) à l'étape r où le paramètre
 vaut r.
Soit {P} le vecteur des accroissements des charges extérieures, appliquées aux
nœuds, lorsque le paramètre  passe de r à r+1
r+1 = r +  .
L'organigramme de résolution d'une étape de calcul sous l’accroissement des charges

extérieures, sous l’action éventuelle de la précontrainte à la mise en tension et éventuellement
aux glissements des éléments câbles sur les déviateurs est le suivant :
1. On suppose connue la valeur des accroissements de déplacements U o des
nœuds. Cette valeur est prise égale à {o}, pour la première itération, ensuite c'est celle
fournie par l'itération précédente.
54
2. Le vecteur SN des accroissements de déplacements aux extrémités d'un élément
est calculé à partir de U o . On calcule la valeur des accroissements d'efforts intérieurs
FN aux extrémités des éléments, dans leur repère intrinsèque, en utilisant pour chaque
élément, la relation (2.94). La valeur de la matrice  KN  est celle obtenue à l'itération

précédente.
3. La connaissance de FN et des déformations {} des sections, correspondant à
l’étape précédente, nous permet d'étudier l'équilibre des sections. En effet le vecteur des
accroissements de sollicitations FS dans une section de calcul est déterminé à partir de
FN en utilisant la relation (2.80). L'action éventuelle de la précontrainte sur la section est
donnée par les relations (2.59) et (2.60). L’équilibre d'une section est exprimé en résolvant le
système non linéaire (2.62) par le procédé itératif (méthode de substitution utilisant les
matrices cordes. Cependant, si l'équilibre d'une section n'est pas vérifié, ceci provient du fait
d'une rupture de la section, ou du fait que le nombre d'itérations fixé a priori n'est pas
suffisant pour aboutir à la convergence du processus itératif. Dans ce cas le calcul est repris
avec un pas de chargement P plus petit.
4. Du fait de l'équilibre des sections, l'accroissement des déformations {} est connu.
L'intégration selon la méthode numérique donne la matrice de rigidité K N  de chaque
élément ainsi que le vecteur PN.
5. La matrice  KU et le vecteur AU dans le repère local intermédiaire s’en
déduisent par les équations (2.117) et (2.121).
6. La matrice  KL et le vecteur AL dans le repère local sont ensuite obtenus par les
relations (2.119) et (2.120).
7. Le passage à  KX et AX dans le repère absolu, lié à la structure, est effectué par
l'intermédiaire de la matrice de rotation [RT].
8. La matrice de rigidité [K] et le vecteur second membre {A} dû à l'action éventuelle

de la précontrainte à la mise en tension des câbles sont construits en écrivant l'équilibre des
nœuds dans le repère absolu lié à la structure.
9. Les accroissements de déplacement des nœuds de la structure {U} correspondant à
l'itération en cours s'obtiennent en résolvant le système (2.35) mentionné plus haut, en tenant
compte des conditions aux limites de la structure.
55
10. Après avoir obtenu l'accroissement de déplacement des nœuds {U}, on procède
à la vérification de la convergence. Si le test de convergence n'est pas vérifié, le calcul est
répété avec la nouvelle valeur {U}. Toutefois un défaut de convergence peut apparaître du
fait de l'instabilité d'ensemble ou du fait de ce que le nombre d'itérations fixé a priori n'est
pas suffisant pour assurer la convergence. Dans ce cas ; le calcul est repris avec un pas de
chargement plus petit.
11. À la fin de chaque étape on calcul les déformations totales de chaque section, Les
efforts totaux FU et les déplacements totaux SU aux extrémités des éléments dans le
repère local, ainsi que les déplacements totaux des nœuds {U} et les charges extérieures
totales {P} dans le repère absolu lié à la structure, en cumulant les accroissements de cette
étape avec les valeurs de l'étape précédente.
12. Avant de passer au test de glissement des câbles sur les déviateurs, on effectue
une étape de correction des déplacements des nœuds {U} obtenus. Cette opération est
indispensable pour le cas des structures sensibles aux effets du second ordre. En effet, lors de
la recherche de l'accroissement des déplacements {U} des nœuds, les matrices [B] et [D]
exprimant la non-linéarité géométrique, sont gardées constantes pendant l'étape, avec les
valeurs obtenues à la fin de l'étape précédente. Cette procédure a pour effet de sous-estimer
l'ensemble des déplacements dus aux effets du second ordre. Pour pallier à ce défaut, on
effectue donc une étape de correction dans l'étape en cours, en y introduisant l'effet d'un
accroissement de charge fictif {P*} qui est la différence entre les charges réelles {P} et les
charges {P*} qui sont en équilibre avec les efforts {Fx} calculés, compte tenu des
déplacements des nœuds {U} calculés :
P * = P - P *
Le vecteur des charges {P*}, est obtenu en écrivant, pour chaque nœud, la somme des forces
{Fx} exercées par le nœud sur les éléments qui le croisent. Toutefois, les réactions d'appuis
qui résultent de l'équilibre des nœuds ne sont pas prises en compte dans ce vecteur.
L'algorithme est dispensé de l'étape de correction si la norme euclidienne f des
efforts est inférieure ou égale à une certaine tolérance Tf choisie à priori :
56
P* . P
1/ 2
 f =   Tf
  P . P 
Dans le cas où la norme nf est supérieure à Tf, une étape de correction suffit en
général pour obtenir une valeur correcte des déplacements des nœuds.
13. Avant de continuer avec un nouveau pas de chargement extérieur, on vérifie la
convergence sur les efforts des éléments câble. Si cette convergence n’est pas respectée, une
étape de correction est alors effectuée pour équilibrer la modification des surtensions dans les
éléments câbles induite par leur glissement sur les déviateurs et leur conséquence sur
l’ensemble de la structure.
Pour cela on introduit le vecteur de surtension des éléments câbles due au glissement
{P*} défini par :
{P*}= {P}-{P*}
Où {P} est le vecteur des efforts dans les éléments câble dans l’étape en équilibre et {P*} le
vecteur des efforts dans les éléments câble après leur glissement sur les déviateurs.
Le vecteur {P*} est remplacé dans l’équation qui s’écrit :
{P*}= K  {U}
14. Après cette étape de correction, on vérifie si ces modifications de surtension des
éléments câbles ne produisent pas de nouveaux glissements sur les déviateurs. S’il n'y a pas
de nouveaux glissements, on continue l’analyse avec une nouvelle étape de chargement. Si
non, on recommence avec la séquence 12 jusqu’à garantir la convergence des efforts.
15. La charge de ruine est obtenue lorsque l'algorithme de résolution diverge pour un
accroissement de charge inférieur ou égal à une précision fixée à priori.
II.6. Élément poutre
Les éléments d’une structure sont modélisés par des éléments poutre à six degrés de
liberté. La section transversale de ces éléments peut être en béton armé, en béton confiné, en
béton précontraint (câbles intérieurs) ou mixte. La figure (II.15) montre une section
transversale d’un élément poutre en béton précontraint. Le comportement non linéaire du
béton, des armatures passives et des armatures actives, peut être modélisé en discrétisant la
57
section en une succession de fibres et en adoptant les lois de comportement réelles des
matériaux.
La section est étudiée dans l’hypothèse des petites déformations et du comportement
élastique non-linéaire. Elle est définie dans son repère principal.
La section droite d’un élément est sollicitée en flexion. L’analyse non linéaire
concerne aussi bien les sollicitations normales dues à un effort normal N et à un moment de
flexion Mz que les sollicitations tangentes dues à un effort tranchant Ty.
II.6.1. Étude d’une section droite

Nous nous intéressons à l'étude d'une section droite d'un élément de poutre en béton
armé, en béton confiné et/ou précontraint. La section est étudiée dans l'hypothèse des petites
déformations et du comportement élastique non linéaire.
En admettant que les sections restent planes après déformation, que les aciers passifs et
les câbles de précontrainte intérieurs après mise en tension, sont parfaitement solidaires au
béton, que les aciers passifs et actifs ne reprennent que les efforts axiaux.
Figure II.24: Description d’une section typique d’une poutre en béton armé et précontraint
L'effet de la précontrainte est introduit en tenant compte du phasage éventuel sur la

mise en tension des câbles. A la mise en tension, l'action de la précontrainte, sur la section, est
un effort de compression, égal et opposé à l'effort de traction dans le câble. Cet effort est
appliqué suivant la tangente à la ligne moyenne du câble, au point où il traverse la section.
Les câbles de précontrainte deviennent des éléments résistants, immédiatement après
leur mise en tension, et possédant une pré-déformation par rapport à la section.
58
Dans le cadre de cette étude, l'analyse non-linéaire concerne aussi bien les
sollicitations normales dues à un effort normal N et à un moment de flexion M, agissant dans
le plan zy. Que la sollicitation tangente, dans le béton, due à l'effort tranchant T dans le plan
zy. On admet que la distorsion due à l’effet de T influe peu sur les déformations normales de
la section.
Compte tenu de l'hypothèse des sections planes et de l'adhérence parfaite entre les
matériaux constituants la section et on négligeant l’influence de la distorsion sur les
déformations normale de la section, la déformation longitudinale  comptée positivement dans
le sens des raccourcissements en un point courant de la section est donnée par :
 =  g + y m . z (2.36)
Où :
g est l’allongement unitaire au centre de gravité de la section, on désigne par z la rotation
unitaire de flexion autour de l'axe Gy et ym l’ordonnée du centre de gravité de la fibre d’ordre

m.
La déformation tangente notée, y, désigne la déformation de cisaillement de la section
dans le plan xy.

Les déformations normales (g , z) peuvent s'écrire sous forme d'un vecteur  n , défini par

n = g , z T (2.37)
De plus, la déformation tangente peut s'écrire sous forme :

t =  y (2.38)

Nous notons  n le vecteur contenant l'accroissement de déformation de la section

droite sous sollicitations normales et t la déformation de la section droite sous
sollicitations tangentes.

n = g , z T (2.39)
59

 t =   y (2.40)
 
A une déformation  n et t de la section droite correspond, dans une armature
adhérente et inclinée à la section, une déformation si, défini par Nait Rabah [53]. Dans
l’hypothèse des armatures passives et actives qui ne reprennent que les efforts axiaux son
expression est :

 si  cos 2 i 1 , ysi  n (2.41)
où  est l’angle que fait la ligne moyenne du câble avec l'axe Gx normal à la section droite.
Dans le cas d'une armature active, la mise en tension produit sur chaque section de
l’élément un effort de compression opposé à l’effort de traction du câble, appliqué suivant la
tangente à la ligne moyenne du câble au point de passage de celui-ci dans la section. La
relation (2.41) doit être corrigée pour tenir compte de la pré-déformation p* du câble par
rapport à la section droite. La déformation p* dans le câble correspond à l'état de la
décompression du béton. Elle exprime la différence entre la déformation du câble et celle de
la section, dans la direction du câble, au moment de la mise en tension. La déformation dans
l'armature de précontrainte d'ordre i et également définie par Kachi [52] elle s’écrit dans les
hypothèses admises dans cette étude sous forme :

 pi  cos 2 i 1 , y pi  n   pi* (2.42)
II.6.2. Calcul des sollicitations internes

Les sollicitations normales à la section droite dues aux matériaux constituant la

section, peuvent s'écrire sous la forme d'un vecteur Fn , définit dans le repère principal Gzy , lié
à la section droite par :

  =  1y 
Nm
Fn = M m
A
m dydz (2.43)
60
Où :
m est la contrainte normale en un point courant de la section transversale A de
l’élément.
La section droite est composée généralement de béton et des armatures passives et
actives. On peut écrire comme suit en tenant compte de la contribution de chaque matériau
constituant la section:
   
Fn  Fnb  Fns  Fnp (2.44)

Où F nb représente la contribution du béton. Ce vecteur est donné par :

 
Fnb =  y1m  m dydz
A
(2.45)
La contribution des armatures passives, considérée comme concentrée dans leur lit en

F
leur centre de gravité, est donnée par le vecteur ns défini par :
 nsi
 
Fns =  y1si  si A si
i 1
(2.46)
Où si est la contrainte normale au centre de gravité de la section de l’armature passive

d’ordre i, Asi la section droite de l’armature et nsi le nombre d’armatures passives traversant la
section droite.
La contribution des armatures actives, considérées comme concentrées dans leur lit en

leur centre de gravité, est donnée par le vecteur F np défini par :

 
npi
Fnp =  y1pi  pi A pi cos i (2.47)
i 1
Où :
pi étant la contrainte normale au centre de gravité de la section de l’armature active d’ordre i,
Api la section droite de l’armature,
i l’angle que fait la ligne moyenne du câble d’ordre i avec l’axe x,
npi le nombre d’armatures actives traversant la section droite.

La variation des déformations normales notée,  n de la section droite, produit une
variation de sollicitations normales dans les matériaux de la section.

La variation des sollicitations normales Fnb dans le béton est donnée par :
61
 y m  dy dz   g 
Fnb  E m  1   z  (2.48)
A
 y m y 2m   
Où
Em désigne le Module élastique longitudinal en un point courant de la section du béton. Il relie
l'accroissement de la contrainte normale m à l'accroissement de déformation longitudinale
m correspondant au point courant.
 m = Em  m (2.49)
La section A étant définie comme une succession de couches, les intégrales sur la
section sont effectuées dans chaque couche en utilisant la formule de SIMPSON.
La variation de la sollicitation tangente due au béton est calculée en fonction des
déformations correspondantes en supposant que le comportement est élastique non linéaire.
Elles peuvent s'écrire sous forme :
Tym  K 2b   (2.50)
Où :
 t est la déformation de la section sous sollicitations tangentes
K2b = Tb est la rigidité de la section à l'effort tranchant dans le plan xz.

Où le terme Tb constitue le module corde de cisaillement de la section.


En tenant compte des relations (2.44) et (2.46), nous pouvons écrire sous forme
matricielle la relation liant les accroissements de déformations dans une section droite, et les
accroissements de sollicitations internes développés dans le béton :

 Fnb  K1b o   n 
Fb      (2.51)
 Tb   o T K 2b    
1 y m  dydz
Avec K1b = E
A
m
 y m y 2m 
62
o  0
0
K2b = Tb

La variation des déformations normales {n} de la section droite produit une
variation des sollicitations normales {Fns} dans les armatures passives donnée par :
ns
Fns   E si A si  1
1
 y si
 
y si   g
y si2   z (2.52)
Où Esi est le module élastique au niveau de l’armature passive d’ordre i. Il relie

l’accroissement des contraintes a l’accroissement des déformations longitudinales
correspondante au niveau du centre de gravite de l’armature.
 si  E si  si (2.53)
En négligeant la sollicitation tangente due aux armatures passives, la relation liant les
accroissements des déformations aux accroissements des sollicitations internes développées
dans les armatures passives s’écrit :

Fs  Fns    K1sT
 Ts   o 0 
 
o   n (2.54)
où [K1s] est la matrice symétrique (2x2) définie par :
ns  1 ysi 
K1s  1 E si A si  
 ysi ysi2 
La variation des sollicitations normales à la section droite dues aux armatures de

précontrainte adhérentes inclinées et supposées concentrées en leur centre de gravité, elle peut
s'écrire sous la forme d'un vecteur Fnp , en négligeant l’influence de la distorsion sur les
déformations longitudinales au niveau de l’armature de précontrainte, comme suit :
63
 np 1 y pi   g 
Fnp  1 E pi A pi cos i  y pi y2   
 (2.55)
 pi  z 
L’équation qui peut s’écrire sous la forme :
 
Fnp  K1p  n (2.56)
Avec
ns  1 y pi 
K1p = 1 E pi A pi 
 y pi
 cos  i
y 2pi 
La variation des sollicitations tangentes dans l’armature de précontrainte d’ordre i

est donnée par :

Fnp  E
np
1
pi  
 g
A pi sin i cos 2 i 1 y pi  z (2.57)
L’équation qui peut s’écrire sous la forme :
 
Ftp  K12p  n (2.58)
Avec
np
K 2p   E pi A pi sin i cos i2  1 
1  y pi 
Epi : Module élastique longitudinal de l'armature de précontrainte d'ordre i. Il relie
l'accroissement de contrainte pi à l'accroissement de déformation pi au niveau de
l'armature d'ordre i :
pi = Epi pi
La relation matricielle liant les accroissements de déformations de la section aux

accroissements de sollicitations internes développées dans les armatures de précontrainte
s’écrit:
64
 

 F  K 1p o    n 
Fp  Ftp  K T 0 
np
 2p     (2.59)
Avec 
o  00
II.6.3 Déformation de la section à l’équilibre
La section droite est supposée soumise :

- d’une part aux accroissements d'efforts sollicitant dus aux forces extérieures.
 N 
  
Fs  M z  (2.60)
 
 Ty 
- d’autre part à l'action de la précontrainte, lors de la mise en tension des câbles :
 np2  cos i 
Ps  -  0i A pi  y pi cos i  (2.61)
1  sin i 
Où : oi est la contrainte initiale dans le câble d’ordre i après pertes instantanées, np2 le
nombre de câbles de précontraintes traversant la section concernés par la mise en tension

L'équilibre de la section droite se traduit par l'égalité des accroissements d'efforts
sollicitants et les accroissements d'efforts internes.
    
Fs  Ps  Fb  Fa  Fp (2.62)

Le vecteur Ps contenant l’action de la précontrainte à la mise en tension des câbles est
non nul. Ces composantes sont définies par la relation (2.43). Après la mise en tension des
câbles intérieurs ce vecteur devient nul et ne participe à l’équilibre de la section que si de
nouveaux câbles sont mis en tension à une étape de calcul quelconque. Dans le cas ou l’en

tient compte des effets différés, le vecteur Ps est non nul. Ces composantes sont alors
fonction des contraintes et du temps.
65
La condition d’équilibre de la section, tenant compte des relations (2.48), (2.50) et

(2.52), peut être exprimée par la relation :
  
Fs  Ps  K s   (2.63)
[Ks] : Représente la matrice de rigidité de la section, donnée par :
K1b  K1s  K1p o 

Ks  =   (2.64)
 K T2p K 2b 

   n
y
Pour un accroissement donné d'efforts sollicitant, la résolution en déformation

d’équations (2.55) est itérative. La solution peut s’exprimer par :

  F
  K s
-1
s  Ps  (2.65)
Où K s  représente la matrice de souplesse de la section.

1
II.6.4. Méthode itérative de résolution et critères de convergence
En élasticité non-lineaire, il est facile de calculer les efforts développés par une
déformation donnée. Cependant il n'existe pas de méthode directe permettant de calculer les
déformations produites par des efforts donnés. On est donc contraint à utiliser une méthode
itérative pour la recherche de la solution en déformation.
Considérons un état d'équilibre de la section soumise à des sollicitations extérieures,

qui peuvent s'écrire sous forme d'un vecteur Fs , définit par :

Fs =  N , M z ,Ty T (2.66)

On note par  s , le vecteur de déformation de la section correspondant. Il est donnée par :
66

 s = g ,  ,  T (2.67)

Pour une variation des sollicitations extérieures Fs , et sous l’action éventuelle de la

précontrainte à la mise en tension des câbles, l'accroissement de déformation  s de la
section est recherché, en utilisant la méthode des substitutions successives utilisant les
matrices cordes.
  
Cette méthode consiste à construire une suite de solutions,  0 , 1 ,......,  i
 
 i Étant calculé à partir de  i 1 , en résolvant le système linéaire :
  
FS  PS  K S i -1  i (2.68)
où  KS i 1 est la matrice de rigidité de la section, à l'itération (i-1)
En partant de la solution initiale nulle, le schéma de résolution s'écrit :
 
 0  0
 
  
 i  K S i 1 FS  PS
1
(2.69)
De façon générale l’équilibre d’une section est défini par deux conditions :
- La limitation de la variation des déformations entre deux itérations successives à

une tolérance fixée a priori. Cette variation est déterminée avec la norme
euclidienne des déformations :
i -  i-1 .  i -  i 1 

1/ 2
 =   Tu
 s +i . s +  i  
(2.70)
Où s est la déformation à l’étape équilibrée précédente.

Généralement une dizaine d’itérations suffisent pour assurer la convergence.
II.6.5. Organigramme de la recherche de l’équilibre d’une section

Nous définissons ici l’organigramme simplifié de la recherche de l’état d’équilibre

d’une section transversale lors d’un accroissement de sollicitations extérieures  F s et de

l’action due à la précontrainte PS à la mise en tension des câbles.
67
Etape équilibrée

FS

PS

 0
 
Q  Q /2
Formation de la matrice de rigidité de la
oui section
K s 
  t
non non
Test sur le det de la
matrice K s 
positif
Résolution du système linéaire
divergence   
Qs  Ps  K s  
Convergence de la  
 0  
oui norme des déformations non
Rupture
Sortie
Figure II.25 : Organigramme simplifié de la recherche de l’équilibre de l’état d’équilibre

d’une section.
II.6.6.Calcul de la rigidité de cisaillement de la section

En élasticité non linéaire, l’étude de l’équilibre d’une poutre en béton armé et précontraint
passe par la résolution d’un système d’équations de la forme :
[P] = [K] . [U] (2.71)
68
Où [P] est le vecteur représentant l’accroissement du chargement appliqué à la poutre, [U]

est le vecteur représentant l’accroissement des déplacements des nœuds de la poutre et [K]
représente la matrice de rigidité globale de la poutre, qui est construite à partir des matrices de
rigidité des sections [Ks].
En élasticité linéaire (module de Young E), l’expression de cette matrice est :
 E .A 0 0 
[Ks] =  0 E .I 0  (2.72)
 
 0 0 G.A 1 
Dans laquelle interviennent l’aire A et l’inertie de flexion I de la section, le module

transversal G = E / 2(1+) et l’aire réduite A1 pour l’effort tranchant. En élasticité non
linéaire, pour une poutre en béton armé, l’expression qui sera utilisée est :
 N N 0 
  u  w 
[Ks] =  M M 
  u  w 0 
 V 
 0 0
  moy 
 
Où N est l’accroissement de l’effort normal de la section, M est l’accroissement du

moment fléchissant de la section, u est l’accroissement de la déformation au niveau centre
de gravité de la section et w est l’accroissement de la rotation différentielle (courbure) de
la section.
Ces déformations répondent à l’hypothèse de planéité des sections (Navier-Bernoulli) ; la

déformation longitudinale x est une fonction linéaire de l’abscisse transversale y :
x (y) = u + w . y (2.73)
V est l’accroissement de l’effort tranchant de la section et moy est l’accroissement de la

distorsion moyenne de la section.
Dans beaucoup de logiciels de calcul non linéaire du comportement des poutres en béton
armé, le terme dû à l’effort tranchant est soit totalement ignoré - Grelat (1978) - soit pris égal
à la valeur G.A1 de l’élasticité linéaire - Nait-Rabah (1990).
69
Dans cette étude, on adopte le modèle de calcul développé par kachi et al. [52], [55-60] qui
évalue la rigidité à l’effort tranchant de chaque section :
V
GA1 
moy
À chaque étape de calcul, sans faire intervenir explicitement un module transversal G.

L’algorithme développé simule le comportement en cisaillement jusqu’à rupture d’une
section de poutre en élasticité non linéaire. La rigidité de cisaillement globale de la poutre est
obtenue alors à partir de celles des sections par un processus d’intégration numérique. Il
permet ainsi d’évaluer la rigidité à l’effort tranchant de la poutre à chaque étape du
chargement et les flèches correspondantes.
L’expression adoptée ci-dessus pour [Ks] est déjà une simplification. En effet, les termes
de couplage entre N et w, et entre M et u, sont incontournables dans l’analyse de
flexion composée, mais on a supposé que les termes de couplage entre N et moy (ou V et
u), et entre M et moy (ou V et w), étaient négligeables. Le couplage entre les
déformations dues à V et celles dues à (N, M) évoqué dans l’introduction va se trouver dans la
dépendance de moy avec les contraintes et déformations dues à N et M.
Pour les poutres en béton précontraint, la matrice [Ks] est complétée par des termes dus à
l’effet de la précontrainte, dans l’équilibre de flexion composée et dans l’équilibre d’une
partie de l’effort tranchant par la composante transversale de l’effort des câbles inclinés –
Nait-Rabah (1990).
II.6.7. Domaine d’application
L’étude porte sur des poutres assez élancées, soumises à une flexion plane avec effort
tranchant. La section transversale est supposée constante par tronçons le long de la poutre.
Elle présente une symétrie par rapport à l’axe y, ce qui permet sa décomposition en une
succession de couches selon y. Le principe de Saint Venant est supposé applicable en toute
section, la contrainte agissant dans le sens transversal y est alors nulle. L’hypothèse de Navier
est appliquée aussi en toute section ; on ne peut donc pas traiter des poutres courtes.
Les poutres sont suffisamment armées pour être loin de la fragilité, donc présenter des
fissures bien réparties.
70
II.6.8.Contrainte transversale dans le béton

Kachi et al [52], [55-60] ont montré que la contrainte dans le sens transversale dans le
béton est une compression tan disque celle dans les armatures transversales est une traction,
l’expression reliant ces deux contrainte s’écrit alors :
by = - y ay (2.74)
II.7. Principes de la modélisation et équations générales

II.7.1. Milieux continus équivalents
II.7.1. 1. En déformations
Seules les armatures sont réellement continues. On fait l’hypothèse que les déformations
moyennes du milieu continu équivalent à l’association du béton fissuré et des armatures sont
gouvernées par celles des armatures ax et ay . Les relations générales entre les déformations
moyennes sont :
x = 1 sin2  + 2 cos2  (2.75)

y = 1 cos2  + 2 sin2  (2.76)
 = 2 ( 1 - 2 ) sin  cos  (2.77)
La direction principale D1 correspond à la plus grande déformation principale (valeur

positive = allongement), la direction D2 à la plus petite. L’angle de D2 avec l’axe longitudinal
x est :
 
 = ( GD2 , Gx )
La direction D2 est aussi celle des fissures, du moins au début de la fissuration. Lorsque les
sollicitations croissent, cette direction peut dévier de celle des fissures initiales, ce qui est
physiquement admissible en raison des contraintes de cisaillement transmises par
engrènement des lèvres des fissures. Ceci correspond aussi au concept des bielles
d’inclinaison variable de l’Eurocode 2 (1992).
71
Dans les directions particulières x et y, on a :
ax = x (2.78)
ay = y (2.79)
II.7.1.2. En contraintes
Le béton doit équilibrer seul les contraintes de cisaillement  parallèles à x et y

(on néglige les contraintes qui seraient équilibrées par l’ « effet de goujon » des armatures).
On fait l’hypothèse d’un milieu continu équivalent pour le béton seul. Les relations générales
entre les contraintes du béton sont donc :
bx = b1 sin2 b + b2 cos2 b (2.80)
by = b1 cos2 b + b2 sin2 b (2.81)
b=  = (b1- b2) sin b cos b (2.82)
On fait de plus l’hypothèse de la coïncidence des directions principales des contraintes du

béton avec les directions principales des déformations moyennes :
b =  (2.83)
II.7.2. Compatibilité partielle des déformations du béton avec les déformations

moyennes
On fait l’hypothèse que les déformations du béton dans les directions principales b1 et b2,
qui sont liées aux contraintes principales b1 et b2, sont égales aux déformations moyennes
principales :
b1 = 1 (2.84)
b2 = 2 (2.85)
Dans la direction D1, qui est grosso modo perpendiculaire aux fissures, devraient intervenir
les deux termes de déformation propre du béton et d’ouverture des fissures. On contourne la
difficulté en faisant intervenir une loi de comportement conventionnelle qui relie la contrainte
de traction b1 à la déformation 1 moyenne .
72
Dans la direction D2 qui est grosso modo celle des bielles comprimées parallèles aux
fissures, on a effectivement continuité du béton.
II.7.3. Equations d’équilibre

II.7.3. 1. Equilibre local des couches
Dans les directions x et y, les contraintes x et y résultent de la somme des termes dus
aux armatures et au béton, compte tenu des aires respectives :
x = x ax + bx (2.86)

y = y ay + by = 0
II.7.3. 2. Equilibre global des sections
L’équilibre doit être assuré entre les sollicitations extérieures N, M, V et les résultantes des
contraintes internes Nint, Mint, Vint. Pour la flexion composée, on fait la somme des
contributions des contraintes longitudinales dans les armatures et dans le béton (ce qui est la
traduction de l’équation (2.86)) :
Nint = A
i
xi  axi  b h 
i
i i bxi (2.87)
Mint = A i
xi  axi y ai  b
i
i h i  bxi y i (2.88)
Pour l’effort tranchant, seules les contraintes de cisaillement dans le béton interviennent :
Vint = b
i
i h i  bi (2.89)
II.7.4. Lois de comportement des matériaux

II.7.4.1. Comportement du béton en compression dans la direction D2
La contrainte principale b2 est une fonction des deux déformations principales 1 et 2 :
73
b2 = Eb2 (2 , 1 ) . 2 (2.90)
On se base sur la relation contrainte-déformation uniaxiale de Sargin proposée à l’annexe 1

des Règles BPEL (1999) [62], qui fait intervenir la résistance fc et le raccourcissement
correspondant b0. On la transforme pour traduire le comportement du béton des bielles
comprimées selon D2 et tendu dans la direction D1. La résistance est modifiée selon la relation
proposée par Vecchio et Collins (1986) [63] en fonction de la déformation 1 , mais en
bornant le rapport fc2 / fc :
fc2 1 f c2
 Avec 0.7  1
fc 0.8 - 0.34 1 / b0 fc
Pour les faibles valeurs de l’allongement 1 on doit retrouver la loi de comportement uni-
axiale de compression simple, donc fc2 = fc. On fait l’hypothèse que la déformation
correspondant à fc2 reste égale à b0. La loi de comportement est alors :
_2
 
_
k b    k b'  1  
  (2.91)
_ _2
1  k b  2     k  '
b
Avec : ε  ε 2 et σ  σ b2 ; kb et k’b selon Règles BPEL 1999 [62].

ε b0 f c2
Par analogie pour le béton confiné la déformation correspondant à fc2 =fcc reste égale à cc0.
 7,94 f ' l 2 f ' l    f 

Avec f cc  f c 0   1,254  2, 254 1    et  cc 0   c 0  1  5 cc  1
 
 f c0 f c0    f c0 
La loi de comportement du béton confiné est celle proposée par Bouafia, Y. et Iddir, A. et
Kachi, M.S. [34], [35], [36], [37].
modifiée :
74
_2
 
_
kc   c  k  1   c
'
 cc  c
Pour : 0 ≤ εc ≤ εcc0 (2.92)
_ _2
1  k c  2   c  k   c '
c
Avec :
  c2 σ
c  et σ  b2
 cc 0 cc
f c2
Eb 0   cc 0
kc 
f cc
Ebc 0  11000 3 f cc
k c  kc  1
'
 cc  f cc  E s  c   cc 0  Pour : εcc0 < εc ≤ ε65 (2.93)
Avec :
6  f c20
Es 
0.35  f cc
 65    cc
Es
 cc  0,65 f cc Pour : Pour : ε65 ≤ εc ≤ εccu (2.94)
75
II.7.4.2. Comportement du béton en traction dans la direction D1
La loi contrainte-déformation du béton en traction est linéaire avant la fissuration du béton.

Au-delà, la contrainte de traction diminue avec l’augmentation de la déformation moyenne de
traction laquelle inclut l’effet de l’ouverture des fissures (figure II.17).
1
ft
Eb0
ft 1
Figure II. 26: Comportement du béton en traction.
Avant fissuration du béton :
1   ft 1 = Ebo . 1 (2.95)
Pour décrire la branche décroissante, après fissuration, qui traduit l’influence du

béton tendu entre les fissures sur la déformation moyenne, on adopte la relation proposée par
Belarbi et Hsu (1994) [64] :
0,4
  (2.95 bis)
 b 1  f t  ft 
 1 
Où ft = ft
Eb 0
ft est la résistance du béton à la traction et ft est la déformation de fissuration du béton.
Contrairement à une idée assez répandue (Grelat 1978, par exemple) la contribution de
l’adhérence, donc du béton tendu entre les fissures, ne s’annule jamais totalement (Fouré,
2003) [65]. Par conséquent, cette loi (2.95) couvre l’ensemble du champ des déformations
jusqu’à la rupture des armatures. Elle s’applique aussi dans le cas particulier du
76
comportement uni-axial vis-à-vis des contraintes longitudinales, en remplaçant b1 par bx et
1 par x.
II.7.4.3. Comportement de l’acier
Les comportements des aciers passifs et de précontrainte sont caractérisés par les
types de relations admises par les règles BPEL 1999 :
- armatures longitudinales ax = Ea ( x ) . x (2.96)

- armatures transversales ay = Eat ( y ) . y (2.97)
- armatures de précontrainte px = Eap ( x ) . x (2.98)
II.7.5. Calcul de la distribution des contraintes tangentielles (y)

Pour calculer la distribution des contraintes tangentielles sur la hauteur d’une
section, on associe à chaque section de calcul une section voisine distante de d. Les
contraintes tangentielles sont alors déduites de l’équilibre des couches comprises entre la
section de calcul sous l’état de sollicitation (N, M, V) et la section voisine soumise aux efforts
(N1, M1, V1) tels que :
N1 = N
M1 = M – V.d (2.99)
V1 = V
La distance entre les deux sections est généralement prise égale à h/6 (où h est la hauteur
totale de la poutre) (Vecchio et Collins, 1988)[66]. Une étude paramétrique permet de
constater que la différence entre la valeur de l’effort tranchant extérieur et la résultante des
(y) calculée avec cette distance reste inférieure à 1 %, ce qui est négligeable.
Les deux sections sont analysées de façon à satisfaire chacune aux deux équations
d’équilibre. En appliquant la méthode classique de la RDM, on étudie l’équilibre de
l’ensemble des forces agissant sur une couche d’ordre k présentée à la figure II.18 :
77
Vk G k-1
Fk1 Fk hk
P Vk
Gk
d
Figure II.27:Forces agissant sur une couche de béton entre deux sections [55].
Où Fk et Fk1 désignent les efforts normaux agissant sur la couche k dans les deux sections,
Gk-1 et Gk sont les efforts de glissement résultants sur les faces inférieure et supérieure de la
couche k et Vk est la part de l’effort tranchant équilibrée par la couche k.
Ces efforts Fk sont déterminés à partir des contraintes normales longitudinales moyennes
dans le béton et dans les armatures longitudinales s’il en existe dans la couche. Dans une
couche quelconque i, on a :
Fi = bxi bi hi + axi Aaxi
Les efforts horizontaux agissant sur les faces supérieure et inférieure de la couche
de béton d’ordre k sont alors calculés comme suit :
k 1
G k -1   ( Fi1 - Fi )
i 1
Gk = Gk-1 + Fk1 – Fk
L’équilibre du tronçon permet de calculer la valeur de Vk, en écrivant que le moment des
forces est nul au point P :
( G k -1  G k ) h k
Vk 
2 d
La valeur moyenne de la contrainte tangentielle de la couche k est alors donnée par :
( G k -1  G k )
k  Vk  (2.100)
bk h k 2 d bk
78
II.7.6. Calcul de la distorsion moyenne
La distorsion varie d’une couche de béton à une autre. On définit une distorsion moyenne
de la section notée moy. Son calcul est effectué en remarquant, d’après le théorème des
travaux virtuels, que l’accroissement We du travail de l’effort tranchant extérieur est égal à
l’accroissement Wi du travail des contraintes tangentielles i. Ces derniers peuvent s’écrire :
We = V . moy
Wi = bi . hi . i . i
On doit avoir, m étant le nombre des couches de béton :

m
We = w i
1
m
i bi hi i
D’où :  moy  1 V
(2.101)
II.7.7. Résolution du problème
II.7.7.1. Transformation des équations générales
Pour la résolution du problème on explicite certaines grandeurs à partir des équations

générales précédentes. On tire de (2.75) :
x 2 (2.102)
1  
sin 2  tg 2 
On réécrit aussi les équations (2.80), (2.81) et (2.82) entre les contraintes du béton
sous la forme suivante (en tenant compte de (2.83):
b x = b2 + b1 - by (2.103)

 = tg-1 ( b2 bx )1/2 (2.104)
b2 'by
b2 = b y - τ (2.105)
tgθ
79
II.7.7.2. Algorithme et organigramme de calcul
L’étude de l’équilibre d’une couche de béton armé où x et  sont connus fait appel à un
système de 10 équations (2.74) (2.76) (2.77) (2.90) [(2.91) ou (2.92)] (2.93) (2.102) (2.103)
(2.104) (2.105) à 10 inconnues : y , 1 , 2 ,  ,  , ay , bx , by , b1 , b2 .Pour le résoudre,
on fait appel à une méthode itérative. Pour une distribution donnée des déformations
longitudinales, on suppose connue la déformation principale 2 et on cherche la valeur de
l’angle  qui permet de satisfaire les conditions de compatibilité et d’équilibre de la couche.
Les contraintes tangentielles (y) sont calculées par l’équilibre de deux sections voisines
(équation (2.100)). La déformation longitudinale x est tirée de Navier (2.73), dans laquelle u
et w résultent de l’équilibre global en flexion composée sous N, M selon (2.87) et (2.88),
avec un couplage à l’effort tranchant V. Enfin, la distorsion moyenne de la section est
calculée par (2.101). La première approximation de la distribution des (y) est obtenue à partir
de l’équilibre des deux sections voisines en flexion composée sans interaction de l’effort
tranchant. Après le calcul de l’équilibre des couches, on corrige cette distribution par
itérations, le paramètre de contrôle étant la valeur maximale de .
La résolution complète au niveau d’une section est décrite par les deux
organigrammes ci-joints : l’organigramme général (figure II.19) et le détail de l’étude d’une
couche (figure II.20).
80
Résultats de l’étape n [Ks]n
Accroissement du chargement P
Accroissement des efforts interne N , M , V

Section étudiée Section voisine [2.99]
N , M, V N , M1 = M – d.V , V
Equilibre en flexion composée N , M sans interaction de l’effort V
Choisir une valeur de la courbure w
Choisir une valeur de la déformation u

Calculer la déformation longitudinale x par la relation (2.73)
Calculer les contraintes bx et ax par les relations (2.90) uni-axiale, (2.95) et
1 (2.96) (2.97)
2 Calculer l’effort normal interne par la relation (2.87)
non
Nint = N
oui
Calculer le moment interne par la relation (2.88)
non
Mint = M
oui
Calculer la distribution (y) selon la relation [2.100] qui satisfait à Vint = V
u , w connus donc x(y) connu, (y) connu
Equilibre d’une couche yi à x et  connus (voir organigramme 2)
3
Boucle sur les couches
4 Equilibre en flexion composée avec interaction de l’effort V
5 Calculer les contraintes ax par les relations (2.96) (2.98) ,
bx est connue d’après l’équilibre des couches
6
Itération sur u et w pour obtenir Nint = N (2.87), et Mint = M (ou M1) (2.88)
Recalculer les contraintes tangentielles avec la relation (2.100) qui satisfait à Vint = V
non Vérifier la convergence des

contraintes tangentielles
oui
oui de la section moy par la relation (2.101)
Calculer la distorsion moyenne
Calculer les rigidités de la section N/Δδu ….. ΔV/Δγmoy
[Ks]n+1 de l’étape n+1
Figure II. 28:Organigramme général de l’étude d’une section.
81
Choisir une valeur de la déformation principale 2
Choisir une valeur de l’angle 
Continuité des déformations : Calculer 1 par (2.102) et y par (2.74)
Lois de comportement : Calculer b2 , b1 et ay par (2.90), [(2.95) ou (2.95 bis)] , (2.96) ou (2.97)
Equilibre sens y : calculer by par (2.74)
Continuité des contraintes du béton : Calculer bx par (2.103)
Calculer la nouvelle valeur ’ par la relation (2.104)
non
’ = 
oui
Calculer la nouvelle valeur de ’b2 par (2.105)
non ’b2 = b2
oui
Calculer la distorsion  par (2.77)
Figure II.29:Equilibre d’une couche à x et  connus.
Dans son état actuel, le logiciel n’effectue pas la vérification de convergence des contraintes
tangentielles, c’est-à-dire qu’il n’effectue pas de correction des (y) calculés initialement.
L’équilibre de flexion composée est calculé par la méthode classique des sections planes en
sollicitations uni-axiales, c’est-à-dire la première partie de l’organigramme 1. Le couplage
entre les contraintes normales x et les cisaillements  n’intervient donc que dans le calcul de
l’équilibre des couches, donc dans moy (Figure II.21).
Résultats de l’étape n [Ks]n
Accroissement des efforts N , M , V dans la section de calcul
Equilibre en flexion composée N , M sans interaction de l’effort V : calcul

des rigidités N/u , N/w, M/w et M/u
Calcul de la rigidité à l’effort tranchant de la section V/moy (voir

organigramme 1 en supprimant la boucle 6)
[Ks]n+1 Résultats de l’étape n+1
Figure II.30:Organigramme de la méthode simplifiée de calcul d’une section
82
II.8.Expression de la matrice de rigidité de l’élément dans son repère intrinsèque

Soit un élément poutre quelconque de la structure de nœud origine Io et de nœud
extrémité Jo, et xoyozo le repère local de référence lié à la position initiale de l’élément (figure
II.22). L’origine du système des axes est en Io
Sous l’action des efforts Fixo , Fiso , Mizo , Fjxo Fjyo , Mjzo , les nœuds se sont déplacés de
(ui , vi ) et (uj , vj ) dans le repère local.
Pour tenir compte, dans la matrice de rigidité de l’élément, des effets des non-
linéarités dues au comportement des matériaux constitutifs et de leur association, et des effets
du second ordre dus aux déplacements des nœuds, l’état de déformation de l’élément est
défini dans le repère intrinsèque xyz lié à la position déformée de l’élément, dans l’hypothèse
des petites déformations, des petits déplacements et en négligeant la déformation propre de
l’élément supposé suffisamment court. Les effets du second ordre sont introduits dans le
passage du repère intrinsèque xyz au repère xoyozo défini précédemment.
Soit un élément quelconque chargé en ces extrémités dans son repère intrinsèque lié
aux extrémités de cet élément (figures II.22, et II.23). Cet élément est supposé suffisamment
court en sorte qu’on puisse négliger les effets de second ordre dus à ses déformations propres
dans ce repère.
Figure II.31: Repères d’un élément poutre
83
Figure II. 32: Elément poutre déformé.
Le vecteur contenant ces charges, noté Fx, est exprimé par :

Fx  N i Tiy M iz N j Tjy M jz 
T
(2.106)
Soit Fn le vecteur colonne (3x1) contenant les efforts indépendants dans le repère
intrinsèque.

Fn =  N j , M iz , M jz T (2.107)
La condition d’équilibre de l’élément dans le repère intrinsèque nous permet d’écrire :
 
FX  T  Fn
T
(2.108)
Avec
-1 0 0 1 0 0
 
T   0 1/ L 1 0 1/ L 0 (2.109)
 
 0 1/ L 0 0 1/ L 1 

On désignera par Sn le vecteur contenant les déplacements dans ce repère, soit:

Sn =  e , iz ,  jz T (2.110)
84
Où e représente la variation de longueur de l’élément.
e = L - L0
D’autre part et compte tenu du fait que l’élément est supposé suffisamment court, pour
que les moments dus a l’interaction effort normal-flèche soient négligeables dans le repère
intrinsèque, les efforts dans une section d’abscisse x s’écrivent :
N(x) = Nj
Mz(x) = x/L Mjz - (1-x/L)Miz (2.111)
Ty(x) = (Miz + Mjz )/L
On peut écrire comme suit sous forme matricielle la relation liant les sollicitations
dans une section aux efforts aux extrémités de l’élément dans le repère intrinsèque :
 
FS = L(x) . Fn (2.112)
Avec :

FS = N(x) , M(x) ,T(x) 
T
(2.113)
et
-1 0 0 
 
L(x)   0 -(1-x/L) x/L  (2.114)
 
 0 1/L 1/L 
En négligeant les variations de longueur dans le calcul des variations de sollicitations,

on peut écrire, compte tenu de la relation (2.112) :
 
FS = L(x) . Fn (2.115)
Le théorème des travaux virtuels nous permet d’écrire :
 L 
L(x)
T
Sn  (x) dx (2.116)
0
85
Où
 (x) 
  
(x)  z (x)
 
  y(x) 
En remplaçant ces relations on obtient la matrice de souplesse de l’élément

K n 1 Telle que :
 
Sn  K n  1  
Fn  G n (2.117)
Avec :
L
K n 1  L(x)T K s 1 L(x) dx
0
 L 
L(x)K s  Ps dx
1
Gn  (2.118)
0
L’intégration est effectuée par la règle de SIMPSON, en considérant pour chaque

élément un certain nombre de sections de calcul dans lesquelles les matrices K s  sont
déterminées a chaque étape.
Finalement, la relation liant l’accroissement des déplacements et des efforts aux
extrémités de l’élément s’écrit :
  
Fn  Pn  K n  Sn (2.119)

Où K n  est la matrice de rigidité corde de l’élément dans son repère intrinsèque et où Pn est
donné par :
 
Pn  K n  G n (2.120)
II.9. Expression de la matrice de rigidité de l’élément dans le repère local

Dans ce paragraphe, on établit la matrice de rigidité de l'élément dans le repère local,
lié à la position initiale de l'élément, en tenant compte des effets du second ordre dus aux
déplacements des nœuds. Ces effets s'introduisent dans le passage des efforts et déplacements
86
du repère intrinsèque xyz (indice n) au repère local xoyozo (indice L) lié à la position initiale
de l’élément.
Nous noterons FL et SL les vecteurs contenant respectivement les sollicitations et
les déplacements aux extrémités de l’élément dans le repère local.

FL=Fixo , Fiyo, Mizo, Fjxo, Fjyo, Mjzo
T
(2.121)

SL = u i , v i , izo , u j , v j ,  jzo 
T
(2.122)
Sous l’action des forces Fjxo , Fjyo , Mizo , Mjzo les nœuds se sont déplacés de ( ui , vi ) et
(uj , vj ) dans le repère local. Dans une phase intermédiaire de calcul, on utilise un repère
intermédiaire du système d’axes xo’yo’zo’. Ces vecteurs s’écrivent alors sous la forme :

FU = Fjxo , Fjyo , M izo , M jzo 
T
(2.123)

SU = u , v , izo ,  jzo 
T
(2.124)
Où :
u = uj - ui
v = vj - vi (2.125)
La relation liant les déplacements dans ces deux repères est donnée par :
 
SU = TO . SL (2.126)
De même les vecteurs efforts sont liés par :
 
FL  T0  FU
T
(2.127)
Avec :
-1 0 0 1 0 0
 0 -1 0 0 1 0 
TO = 
 0 0 1 0 0 0
 (2.128)
 
 0 0 0 0 0 1
Les déplacements de l’élément dans le repère intrinsèque sont donnés par :
e  L - L0
iz  izo -  (2.129)
87
 jz   jzo - 
  Arcsin(v/L)
Avec : L  ( L0  u ) 2  v 2
Le déplacement axial e du nœud J est fonction des déplacements relatifs u et v des

nœuds I et J dans le repère intermédiaire. En effet :
e  ( L0  u ) 2  v 2 - L0
Où :
u et v sont définis par la relation (2.125).
La différenciation des expressions (2.129) permet d’établir la relation entre les


accroissements de déplacements nodaux n dans le repère intrinsèque, et les accroissements

des déplacements nodaux  U dans le repère intermédiaire.
 
 n  B  U (2.130)
Où [B] est la matrice de transformation géométrique qui s’écrit :
 e, u e, v 0 0
 
B   -, u -, v 1 0 (2.131)
 
-, u -, v 0 1
En ne tenant compte que des termes du premier ordre, dans le développement en séries
des dérivées partielles, cette matrice s’écrit alors :
 1 v/L 0 0
 
B  v / L 1/ L0  u/L 1 0
 2 2 (2.132)
 
 v / L2 1/ L 0  u/L2 0 1 
 
88
Selon le théorème des travaux virtuels, le travail virtuel des forces appliquées aux
 
nœuds Fn dans le repère intrinsèque dû à un déplacement virtuel Sn est égal au travail des
 
forces appliquées aux nœuds FU du à un déplacement SU , dans le repère local, c'est à dire :
   
FnT Sn  FUT SU (2.133)
Compte tenu de la relation (2.120), nous pouvons écrire :
   
Fn B SU  FUT SU (2.134)

La variation de Sn est arbitraire, d’où :
 
FU  B Fn
T
(2.135)
La différenciation de l’équation (2.135) nous conduit à :
  
FU  B Fn  B Fn
T T
(2.136)
Compte tenu des relations (2.130) et (2.120), nous avons :
   
FU  PU  B K n B SU  B Fn
T T
(2.137)

Où PU est le vecteur contenant l’action de la précontrainte à la mise en tension dans le repère
intermédiaire. Il est donné par :
 
PU  B Pn
T
(2.138)
D’autre part, les termes de la matrice B ne dépendent que des dérivées partielles des
fonctions e et, qui sont fonction des déplacements relatifs u et v des nœuds I et J
respectivement. On a alors :
89
B  B, u u  B, v v

T T T
(2.139)
On peut donc écrire :
 
B Fn  D SU
T
(2.140)
Avec :
 d11 d12 0 0
 
d 21 d 22 0 0
D    (2.141)
0 0 0 0
 
 0 0 0 0
Avec :
d11  N j e, uu - (M iz  M jz ) , uu
d12  d 21  N j e, uv - (M iz  M jz ) , uv (2.142)
d 22  N j e, vv - (M iz  M jz ) , vv
En ne tenant compte que des termes du premier ordre dans le développement des
termes de la matrice D on obtient :
d11  - 2 v (Miz  M jz ) /L30

d12  d 21  - v N j /L20 - (2 u - L0 ) (Miz  M jz ) /L30
d 22  (L0 - u) N j /L20  2 v (Miz  M jz ) /L30 (2.143)
L’équation (2.138) nous permet d’exprimer les efforts appliqués aux nœuds de
l’élément dans le repère intrinsèque en fonction des efforts appliqués aux extrémités de
l’élément dans le repère intermédiaire. En effet :
Fjxo   Nj 
   C   
Fjyo  - ( M iz  M jz ) /L 0 
90
Avec :
 1 -v / L0 
C   
 v / L0 1 - u / L0 
On a également
 Nj  Fjxo 
  C
1
  
- ( M iz  M jz ) /L 0  Fjyo 
Si bien que finalement
 Nj  L 0 (L 0 - u) / B1 v L0 / B1  Fjxo 
      
- ( M iz  M jz ) /L 0    v L 0 / B1 L20 / B1  Fjyo 
Avec
B1  L20 - u L 0  v 2
Alors les termes de la matrice D deviennent :
 d11  - 2 v 2 / L 0 B1 2 v / B1 
    Fjxo 
 d12    u v / L0 B1 B2 / L0 B1   
    Fjyo 
d 22   B3 / L0 B1 - (v L 0  u v) / L 0 B1 
Avec
d12 = d21
B 2  - L20  2 u L 0 - v 2
B3  ( L 0  u ) 2  2 v 2
Si on néglige les termes non linéaires en v et la contribution du déplacement longitudinal, on
obtient :
 d11   0 2 v / L20 
    Fjxo 
d
 12   0 -1 / L 0    (2.144)
    Fjyo 
d 22  1 / L0  v / L20 
91
Et par conséquent l’équation (2.127) peut s’écrire :
   
FU  PU  B Fn  D SU
T
(2.145)
En remplaçant les équations (2.119) et (2.130) dans l’équation (2.145), on obtient la

relation liant les accroissements des efforts à ceux des déplacements dans le repère
intermédiaire. Son expression est alors :
  
FU  PU  ( B K n B  D ) SU
T
(2.146)
Finalement, en substituant les équations (2.119), (2.126) et (2.127) dans l’équation

(2.146), on arrive à l’expression liant les accroissements des efforts et des déplacements des
nœuds de l’élément dans son repère local, lié à sa position initiale, qui fait intervenir la
matrice de rigidité de l’élément dans le repère local notée K L  sous la forme :
  
FL  PL  K L  SL (2.147)
L’expression de la matrice de rigidité de l’élément dans son repère local est donnée par :
K L   To T ( BT K n B  D ) To  (2.148)
Et
 
PL  To  B Pn
T T
(2.149)
II.10. Expression de la matrice de rigidité de l’élément dans le repère global lié à la

structure
La matrice de rigidité de la structure, notée R  , est formée à partir des matrices de
rigidité des éléments définissant la structure dans le repère global lié à la structure XYZ. Pour
permettre l’assemblage de la matrice de rigidité de la structure, il est nécessaire de former les
matrices de rigidité des différents éléments dans le repère global.

On notera FG le vecteur contenant les efforts appliqués aux extrémités d’un élément,

et SG celui contenant les déplacements des nœuds dans le repère global XYZ (indice G). On
a:
92

FG  FiX FiY M iZ FjX FjY M jZ  (2.150)

SG  U i Vi  iZ U j Vj  jZ  (2.151)
La relation liant ces grandeurs (et leur accroissement) dans les deux systèmes s’obtient
au moyen d’une matrice de rotation TG  dont les termes ne dépendent que de l’orientation
initiale des éléments. La relation liant les déplacements des nœuds dans le repère local et les
déplacements des nœuds dans le repère global s’écrit :
 
SL  TG  SG (2.152)
De même celle liant les efforts s’écrit :
 
FG  TG  FL
T
(2.153)
Avec
r0 
TG    (2.154)
  r0 
Où
 cos sin 0
 
r0  -sin cos 0 (2.155)
 
 0 0 1
L’équation (2.147) s’écrit alors en faisant intervenir la matrice de rigidité dans le

repère global :
  
FG  PG  K G  SG (2.156)
Avec
K G   TG T K L TG  (2.157)
 
PG  TG  PL
T
(2.158)
93
Où K G  représente la matrice de rigidité de l’élément dans le repère global lié à la structure

XYZ.
II.11. Expression des accroissements des efforts en fonction des accroissements des
déplacements dans le repère global
La relation liant les accroissements des efforts à ceux des déplacements dans le repère
global, lié à la structure, s’obtient au moyen de la matrice de rotation TG1 définie par :
TG1  r01 0  (2.159)

 0 r01 
Où
 cos sin 0
r01  sin cos 0 (2.160)
 0 0 0
Finalement, l’équation d’équilibre de l’élément câble, dans le repère global lié à la

structure s’écrit :
  
FG  PG  K G1 SG (2.161)
Avec
K G1  TG1
T
K L1 TG1
 
PG  TG1
T
PL (2.162)
II.12. Conclusion
Une méthode de calcul des éléments de section circulaire est développée dans ce chapitre
en y introduisant différentes formes de confinement interne et externe. L’étude est menée en
élasticité non linéaire, il est tenu compte du comportement non linéaire des différents
matériaux constituant la section transversale. Le confinement du béton est introduit par le
biais d’une loi de comportement également non linéaire, du béton confiné et du béton armé
confiné tenant compte du niveau de confinement, du type de confinement et de l’ensemble des
sollicitations (effort normal, moment fléchissant et effort tranchant). Enfin, la méthode de
calcul est injectée dans un logiciel de calcul non linéaire des éléments de structures afin
d’étudier la contribution et l’influence du confinement du béton sur le comportement global
des éléments de structures de section circulaire. Les résultats de la confrontation feront l’objet
du chapitre III.
94
Chapitre III Application du modèle et résultats de calcul
CHAPITRE III
APPLICATION DU MODELE ET
CONFRONTATION DES RESULTATS
III.1. Applications et confrontations dans le cas des sections en béton armé
III.1.1. cas des sections circulaires

Cette application est effectuée sur une série d’essais effectués par J .B. Mander et al. [23-
24] dont les propriétés des sections d’étude et des données expérimentales sont résumées dans
le tableau III.1:
Tableau III.1:Caractéristiques des données expérimentales des sections circulaires en béton.
fc0 ds s Asl Asp s
N° Essai
(MPa) (cm) (cm) (cm2) (cm2) ( %)
1 29 43 ,8 4,1 24,12 1,13 2,5
2 29 43 ,8 6,9 24,12 1,13 1,5
3 29 43 ,8 10,3 24,12 1,13 1,0
4 29 44,0 11,9 24,12 0,79 0,6
5 29 44,0 3,6 24,12 0,79 2,0
6 29 43,4 9,3 24,12 2,01 2,0
7 32 43 ,8 5,2 49,24 1,13 1,98
8 30 43 ,8 5,2 49,74 1,13 1,98
9 32 43 ,8 5,2 50,24 1,13 1,98
10 30 43 ,8 5,2 48,23 1,13 1,98
Les confrontations des courbes obtenues par le modèle et les courbes expérimentales
proposées par J .B.Mander et al. [23-24] pour le cas d’une section circulaire sont représentées
sur les figures suivantes :
95
courbe contrainte-deformation (essai N° 1) courbe contrainte-déformation (essai N° 2)

60
50
50
40
c ontrainte (M P a)
c ontrainte (M P a)
40
30
30
20 20
expérimentale
10 expérimentale
10 Modèle
Modèle
0 0
0 0.02 0.04 0.06 0.08 0 0.02 0.04 0.06
déformation déformation
Figure III.33:Confrontation de la courbe théorique avec la courbe

Expérimentale (essais N° 1 et essais N°2)
courbe contrainte-déformation (essai N° 4)

courbe contrainte-déformation (essai N° 3) 40
expérimentale
50
Modèle
40
30
c o ntrainte (M P a)
c ontrainte (M Pa)
30 20
20
expérimentale 10
10
Modèle
0 0
0 0.02 0.04 0.06 0.00 0.02 0.04 0.06
Figure III. 34:Confrontation de la courbe théorique avec la courbe

courbe contrainte-déformation (essai N° 5) courbe contrainte-déformation (essai N° 6)

60 60
50 50
contrainte (MPa)
contrainte (MPa)
40 40
30 30
20 20
expérimentale
10 10 expérimentale
Modèle Modèle
0 0
0 0.02 0.04 0.06 0.08 0 0.02 0.04 0.06 0.08
déformation déformaton
Figure III. 35 : Confrontation de la courbe théorique avec la courbe

96
courbe contrainte-déformation (essai N° 7) courbe contrinte-déformation (essai N° 8)

60 60
50 50
c ontrainte (M P a)
c ontrainte (M Pa)
40 40
30 30
20 20
expérimentale expérimentale
10 10
Modèle Modèle
0 0
0 0.02 0.04 0.06 0.08 0 0.02 0.04 0.06 0.08
Figure III.36: Confrontation de la courbe théorique avec la courbe

courbe contrainte-déformation (essai N° 9) courbe contrainte-déformation (essai N° 10)

60
60
50
50
c ontrainte (M P a)
40
c ontrainte (M P a)
40
30 30
20 20
expérimentale
10 expérimentale
10
Modèle
Modèle
0 0
0 0.02 0.04 0.06 0 0.02 0.04 0.06

Remarque
On remarque que les courbes calculées à l’aide du modèle actuel sont proches des courbes
expérimentales : le modèle est assez représentatif. Il faut noter que, dans le cas des sections
étudiées, le pourcentage d’acier de confinement varie de 0,6 à 2,5 %.
Le tableau III.2 résume les résultats de calcul de la résistance à la compression des sections
circulaires des essais menés par J .B.Mander et al. [23-24] en fonction du pourcentage de
confinement.
97
Tableau III.2: Caractéristiques des données expérimentales des sections circulaires en béton.
N° Essai 1 2 3 4 5
s
2,5 1,5 1,0 0,6 2,0
( %)
fcc
51,88 43,35 38,50 34,51 47,99
(MPa)
La figure III.6 nous montre la superposition des courbes calculées par le modèle proposé
en termes de contrainte-déformation des sections circulaires des essais menés par J .B.Mander
et al. [23-24] .On remarque que la contrainte maximale est proportionnelle au pourcentage de
confinement.
Supperpotion des courbes théoriques  = f(e)

60
50
40
contraintes
30
20
(ρs-fcc)=(2,5-51,88)
(ρs-fcc)=(2,0-47,99)
10 (ρs-fcc)=(1,5-43,35)
(ρs-fcc)=(1,0-38,50)
(ρs-fcc)=(0,6-34,51)
0
0.000 0.020 0.040 0.060 0.080
déformations
Figure III.38 : Superpositions des courbes théoriques  = f(e)
On observe clairement que les courbes de comportement se superposent graduellement

avec le pourcentage de confinement ainsi que la contrainte maximale et la déformation à la
rupture .Ce qui signifie que le comportement des sections est proportionnel au pourcentage de
confinement.
98
III.1.2. cas d’une section rectangulaire

Les données expérimentales et propriétés des sections rectangulaires d’étude sont
résumées dans le tableau III.3 :
Tableau III.3: Caractéristiques des données expérimentales des sections rectangulaires en béton.
N° fc0 Asl Asp dc bc w s s’ ρs

ρcc ρsy
essais MPa cm2 cm2 cm cm cm cm cm (%)
1 26 18,09 0,28 9,40 7,74 2,50 1,90 0,0229 0,0238 0,0139 3,77
2 29 18,09 0,28 9,40 7,74 2,50 1,90 0,0229 0,0238 0,0139 3,77
3 26 18,09 0,28 9,40 7,74 5,00 4,40 0,0229 0,0119 0,0122 2,41
4 26 18,09 0,28 9,40 7,74 2,50 1,90 0,0229 0,0238 0,0139 3,77
10 43 32,15 0,28 9,00 7,28 2,50 1,90 0,0558 0,0249 0,0140 3,89
11 43 32,15 0,28 9,40 7,28 5,00 4,40 0,0531 0,0119 0,0122 2,41
12 43 32,15 0,79 9,00 7,17 4,20 3,20 0,0558 0,0418 0,0235 6,53
13 43 11,30 0,79 9,00 14,50 3,00 2,00 0,0196 0,0558 0,0206 7,91
Les confrontations des courbes obtenues par le modèle et les courbes expérimentales
proposées par J .B.Mander et al. Pour le cas d’une section rectangulaire sont représentées sur
les figures suivantes :
courbe contrainte -déformation (essai N°1) courbe contrainte -déformation (essai N° 2)

60 60
50 50
contrainte (M Pa)
c ontrainte (M P a)
40 0.65 fcc 40
30 30
0.65 fcc
20 20
expérimentale
10 expérimentale 10
Modèle
Modèle
0 0
0 0.02 0.04 0.06 0 0.02 0.04 0.06
Figure III.39 : Confrontation de la courbe théorique avec la courbe

Expérimentale (essais N°1 et essais N°2)
99
courbe contrainte -déformation (essai N°3) courbe contrainte -déformation (essai N°4)
50 60
40 50
contrainte (MPa)
contrainte (MPa)
40
30 0.65 fcc
0.65 fcc 30
20 20
expérimentale
Modèle
Modèle
0 0
0 0.02 0.04 0.06 0 0.02 0.04 0.06 0.08

courbe contrainte -déformation (essai N° 10) courbe contrainte -déformation (essai N° 11)
80 70
70 60
60
contrainte (MPa)
50
contrainte (MPa)
50
0.65 fcc 40
40 0.65 fcc
30 30
expérimentale
10 Modèle 10
Modèle
0 0
0 0.02 0.04 0.06 0 0.01 0.02 0.03 0.04

courbe contrainte -déformation (essai N° 12) courbe contrainte -déformation (essai N°13)
90 80
80 70
70 60
contrainte (MPa)
contrainte (MPa)
60 50
50
40
40
30
30
20 expérimentale 20 experimentale
10 Modèle
10 modèle
0 0
0 0.02 0.04 0.06 0.08 0 0.02 0.04 0.06

100
Remarque :
On remarque que les courbes calculées à l’aide du modèle actuel sont proches des
courbes expérimentales : le modèle est assez représentatif. Il faut noter que, dans le cas des
sections étudiées, le pourcentage d’acier de confinement varie de24 à 79 ‰24 à 79 ‰.
III.1.3. cas d’une section béton enrobées de tube en aciers (CFST : Concrete filled steel
tubes)
Cette application est effectuée sur une série d’essais proposés par Lahlou [67] .les sections
sont en béton enrobées de tube en aciers(CFST) dont Les propriétés des données
expérimentales sont résumées dans le tableau III.4:
Tableau III.4: Caractéristiques des données expérimentales des sections circulaires en béton
fc0 ds tj Fuj
N° Essai
(MPa) (mm) (mm) (MPa)
1 44 150 3.2 457
2 95 150 6,4 321
3 131 150 6.4 321
Les confrontations des courbes obtenues par le modèle, les courbes expérimentales
proposées par Lahlou et al. [67] pour le cas d’une section circulaire sont représentées sur les
figures Suivantes III.11, III.12 et III.13 :
140
120
100
Contrainte [MPa]
80
60
Lahlou et al.
40
Baris Binici
20
modele proposé
0
0 0.02 0.04 0.06 0.08 0.1 0.12
déformation
Figure III.43 : Confrontation des courbes théoriques avec les courbes

Expérimentales proposées par Lahlou (essais N° 1)
101
250
200
Contrainte [MPa]
150
100
Lahlou et al.
50 Baris Binici
modele proposé
0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07
déformation
Figure III.44: Confrontation des courbes théoriques avec les courbes

300
250
200
Contrainte [MPa]
150
100 Lahlou et al.

Baris Binici
50
modele proposé
0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08
déformation
Figure III.45: Confrontation des courbes théoriques avec les courbes

La confrontation des courbes obtenues par le modèle de comportement du béton confiné

aux résultats expérimentaux donne une bonne concordance.
102
On remarque, dans cette série d’essais, qu’il y a absence de palier à 0.65fcc ; la déformation
à la rupture pour ces données expérimentales est atteinte dans la zone descendante (voir le
calcul de la déformation eccu formule 2.33 Bis §II.3.3).
Le modèle tient compte des paramètres suivant :
-La résistance à la compression du béton.
-Le diamètre, l’espacement, et le pourcentage des armatures transversales.
-Le diamètre des armatures longitudinales.
-La géométrie de la section en béton
-La limite élastique des armatures de confinement
-l’épaisseur du tube en acier dans le cas d’une section confiné par tube (CFST)
Son application, pour des pourcentages de confinement allant de 6 à 25 ‰ dans le cas des
sections circulaires, et allant de 24 à 79 ‰ dans le cas des sections rectangulaires, est
satisfaisante. L’application du modèle actuel pour des sections confinées par tube en acier est
satisfaisante
III.1.4.Influence des paramètres de confinement sur le calcul des sections
La simulation numérique, concerne une section circulaire en béton confiné soumise à

une flexion composée plane. Le calcul est effectué en élasticité non linéaire jusqu'à rupture,
les algorithmes développés ont été repris pour le calcul des déformations internes à partir
efforts appliqués.
Des applications ont été faites pour des sections en béton armé confiné avec des armatures en
aciers puis une étude comparative a été effectuée.
III.1.4.1. Influence des armatures transversales

Le calcul est effectué pour une section circulaire de diamètre 40 cm. Elle est armée de
barres d’acier de diamètre 12 mm et armatures transversales de diamètre 8mm pour la section
confinée et sans armatures transversales pour la section non confinée [68-69].
Les caractéristiques des matériaux et les courbes de chargement ainsi que les courbes
d’interactions (effort normal N en fonction du moment fléchissant M) obtenues, en utilisant
les lois de comportement réelles des matériaux, sont représentées sur les tableaux et figures
suivantes :
103
Caractéristiques de la section :
Tableau III.5: Caractéristiques du béton section 40cm armé avec 8T12 et des armatures
transversales T8 avec un espacement constant e = 10 cm.
sections
Fcc(MPa) Ftj(MPa) Ebc0(MPa) Kc Kc’ rt cc0 ccu Es(MPa)
Confinée
31.31 2.48 34671 5.01 4.01 -0.002 0.0045 0.05 1874.6
Non
25 2.1 32164.200 2.316 1.316 -0.002 0.002 0.0035 /
confiné
Tableau III. 6: Caractéristiques des aciers.
e(Mpa) Ea (MPa) u r (MPa)

400 200000 0.06 440
La confrontation des courbes moment-courbure avec et sans tenir en compte des armatures
transversales de confinement est représentée sur les figure III.14.
Courbe Moment-Courbure
100
90
80
70
Moment [KN.m]
60
50
40
30
sans conf inement
20
10 avec conf inement
0
0 0.01 0.02 0.03
Courbure (1/m)
Figure III.46: Confrontation des courbes Moment-Courbure calculées avec et

sans confinement
Le moment ultime calculé en tenant compte des armatures transversales de confinement est
égal à 94.63 kN.m, et sans armatures transversales le moment ultime calculé est égal à 88.8
kN .m ; soit une augmentation de 6%.
104
En terme d’effort normal les courbes efforts normal-déformation et la confrontation des

courbes avec et sans confinement sont représentées sur les figures III .15.
Courbe Effort Normal-Déformation
4500
4000
3500
Effort Normal[KN]
3000
2500
2000
1500
1000
avec conf inement
500
sans conf inement
0
0 0.001 0.002 0.003 0.004
Déformation
Figure III.47: Confrontation des courbes Effort Normal -Déformation calculée

sans et avec confinement
On voit bien qu’il y a une augmentation de l’effort normal ultime calculé en tenant
compte des armatures transversales de confinement. On passe d’un effort normal ultime égale
à 3370 kN sans armatures transversales à un effort normal ultime égale à 4175.25 kN avec
armatures transversale. Soit une augmentation de l’ordre de 20 %.
Pour la déformation maximale on remarque un grand écart observé, on passe d’une valeur
de 2.8 % sans armatures transversales à une valeur de 3 .5 % avec armatures transversales.
La confrontation des courbes d’interactions effort normal-moment avec et sans

confinement est représentée sur la figure III.16
Courbe d'intéraction N=F(M)

5000
Sans conf inement
4000 Avec conf inement
3000
Effort Normal [KN]
2000
1000
0
-50 0 50 100 150 200 250
-1000
Moment [KN.m]
Figure III.48 : Confrontation des courbes frontières efforts normal -moment calculée
On remarque qu’il y un gain de capacité portante considérable soit un gain de 20%.
avec et sans confinement
105
III.1.4.2. Influence du diamètre des armatures transversales

Le calcul est effectué pour une section circulaire de diamètre 40 cm. Elle est armée de
barres d’acier de diamètre 12 mm et armatures transversales de diamètre 8mm, 10mm.
suivantes :
Tableau III. 7: Caractéristiques du béton section 40cm armé avec 8T12 et des armatures
transversales T8, T10 avec un espacement constant e = 10 cm.
sections
Fcc(MPa) Ftj(MPa) Ebc0(MPa) Kc Kc’ rt cc0 ccu Es(MPa)
Confinée
31.31 2.48 34671 5.01 4.01 -0.002 0.0045 0.05 1874.6
T10
Confinée
34.59 2.87 35840 6.01 5.01 -0.002 0.0058 0.05 1180.7
T8
Tableau III.8 : Caractéristiques des aciers.

400 200000 0.06 440
La confrontation des deux courbes (confinement en T8 et T10) est représentée sur la figure
III.17 suivante:
Courbe Moment - Courbure

120
100
80
Moment [kN]
60
40
20 Asp:T8
Asp:T10
0
0 0.01 0.02 0.03
Courbure (1/m)
Figure III.49: Confrontation des courbes Moment-Courbure confinée en T8 et T10
106
On voit bien qu’il y a une augmentation du moment ultime avec des armatures transversales
de diamètre 10mm par rapport au diamètre 8mm. On passe d’un moment ultime égale à 94.63
kN .m à un moment ultime égale à 100.49 kN .m, soit une augmentation de 5.8 %.
la confrontation la confrontation des courbes de chargement en termes d’effort normal-
déformation avec un confinement en T8 et T10 sont représentées sur les figures III.18.
Courbe EffortNormal-Déformation
5000
4500
4000
Effort Normal [KN]
3500
3000
2500
2000
1500
1000 Asp:T8
500 Asp:T10
0
0.000 0.001 0.002 0.003 0.004 0.005
Déformation
Figure III.50 : Confrontation des courbes Effort Normal –Déformation confinée

en T8 et T10
Il y a une augmentation de l’effort normal ultime calculé avec des armatures transversales
de diamètre 10mm par rapport au diamètre 8mm. On passe d’un effort normal ultime égale à
4597.65 kN avec des armatures transversales de 10mm à un effort normal ultime égale à
4175.25 kN avec armatures transversale de8mm. Soit une augmentation de l’ordre de 10
% .en terme de déformation ultime soit un écart de l’ordre de 21%.
III.1.4.3. Influence de l’espacement des armatures transversales

Le calcul est effectué pour la même section, de diamètre 40 cm armée de barres de diamètre
12 mm, et des aciers transversale de diamètre 8mm .en faisant varier l’espacement des aciers
transversaux de confinement de 5 à 20 cm, Les caractéristiques du béton et des aciers sont
données aux tableaux III.9, III.10. Les courbes de chargement, les courbes d’interaction
(effort normal N en fonction du moment fléchissant M), et la confrontation en utilisant les lois
de comportement réelles des matériaux, sont représentées sur les figures III.19, III.20 et
III.21 :
107
Tableau III. 9 : Caractéristiques du béton section 40cm armé avec 8T12 et des armatures
transversales T8 avec un espacement variable.
Fcc(MPa) Ftj(MPa) Ebc0(MPa) Kc Kc’ rt cc0 ccu Es(MPa) esp

37.53 2.85 36830 6.88 5.88 -0.002 0 .007 0.05 865.25 5
31.31 2.48 34671 5.01 4.01 -0.002 0.0045 0.05 1874.6 10
28.79 2.34 33785 4.18 3.18 -0.002 0.0036 0.05 3093 15
27.74 2.26 33300 3.72 2.72 -0.002 0.0031 0.05 4560.9 20
Tableau III.10: Caractéristiques des aciers.

400 200000 0.06 440
120
100
80
Moment [KN.m]
60
40 Esp = 5cm
Esp = 10 cm
20 Esp = 15 cm
Esp = 20 cm
0
0 0.01 0.02 0.03
Courbure (1/m)
Figure III.51: Confrontation des Courbes Moment-Courbure confinée avec

espacement variable
On remarque que le moment ultime est plus important dans le cas d’un espacement de 5 et 10
cm par rapport à la valeur du moment ultime pour un espacement de 15 et 20cm. soit une
variation de 8% ; l’espacement est inversement proportionnel à la valeur du moment ultime.
108
Courbe Effort Normal Déformation

5000
4500
4000
Effort Normal [KN.m]

3500
3000
2500
2000
esp=5cm
1500
esp=10cm
1000
esp=15cm
500
esp=20cm
0
-0.001 0 0.001 0.002 0.003 0.004 0.005
Déformation
Figure III.52: Confrontation des Courbes Effort Normal –Déformation
On passant d’un espacement de 5 à 20 cm, on a un gain de déformation et d’effort normal soit

de 0.0024 à 0.0043 pour la déformation soit 44% et de 3389.44 KN à 4595.4 KN en terme
d’effort normal ultime soit une augmentation de 26%.

5000
esp=5cm
esp=10cm
4000 esp=15cm
esp=20cm
3000
Effort Normal [KN]
2000
1000
0
-50 0 50 100 150 200 250
-1000
Moment [KN.m]
Figure III.53:Confrontation des courbes d’interactions effort normal -moment

confinée avec espacement variable
Plus l’espacement des armatures transversales est réduit, on à un gain important de la

capacité portante.
109
III.1.4.4. Influence du diamètre des armatures longitudinales

Le calcul est effectué pour une section circulaire de diamètre 40 cm. Elle est armée de barres
d’acier de diamètre 12 mm, 14mm, 16mm et d’armatures transversales de diamètre 8mm.
suivantes :
Tableau III.11 : Caractéristiques du béton (section 40cm) armé avec 8T12 ,8T14 ,8T16 et des
armatures transversales T8
Fcc(MPa) Ftj(MPa) Ebc0(MPa) Kc Kc’ rt cc0 ccu Dia(mm)

31.31 2.48 34671 5.01 4.01 -0.002 0.0045 0.05 12
31.33 2.48 34679 5.02 4.02 -0.002 0.0045 0.05 14
31.35 2.48 34688 5.03 4.03 -0.002 0.0045 0.05 16
Courbe Effort Normal-Déforamtion

5000
4500
4000
3500
Effort Normal[KN]
3000
2500
2000
1500
1000 Dia=16cm
Dia=14cm
500
Dia=12cm
0
-0.001 0 0.001 0.002 0.003 0.004
Déformation
Figure III. 54 : Confrontation des courbes Effort Normal -Déformation calculée

avec les différents diamètres
On remarque qu’il n’ y’a pas une grande influence du diamètre des aciers longitudinaux
sur la valeur de l’effort normal et la déformation de rupture.
110

5000
Diam=16cm
Diam=14cm
4000
Diam=12cm
3000
Effort Normal [KN]
2000
1000
0
-50 0 50 100 150 200 250
-1000
-2000
Moment [KN.m]
Figure III. 55: Confrontation des courbes d’interaction avec les différents
diamètres
Il n y a pas une grande influence du diamètre des aciers longitudinaux sur la capacité
portante de la section.
III.1.4.5. Influence de la résistance à la compression du béton

Le calcul est effectué pour une section circulaire de diamètre 40 cm de classe de béton
variable (C25, C30, C40). Elle est armée de barres d’acier de diamètre 12 mm et d’armatures
transversales de diamètre 8mm. Les caractéristiques des matériaux et les courbes de
chargement ainsi que les courbes d’interactions (effort normal N en fonction du moment
fléchissant M) obtenues, en utilisant les lois de comportement réelles des matériaux, sont
représentées sur les tableaux et figures suivantes :
Tableau III.12: Caractéristiques de la section en béton (section 40cm) avec les différentes
classes
Fcc(MPa) Ftj(MPa) Ebc0(MPa) Kc Kc’ rt cc0 ccu fc(MPa)

31.31 2.48 34671 5.01 4.01 -0.002 0.0045 0.05 25
36.40 2.78 36456 4.14 3.14 -0.002 0.0041 0.05 30
39.52 3.39 39562 3.09 2.09 -0.002 0.0036 0.05 40
111
120
100
Moment [KN.m]
80
60
40
fc = 40MPa
20 fc = 30MPa
fc = 25MPa
0
0 0.01 0.02 0.03
Courbure (1/m)
Figure III.56: Confrontation des courbes Moment -Courbure calculées avec les
différentes classes du béton
Le moment ultime croit avec l’accroissement de la résistance à la compression du béton.
Courbe Effort Normal-Déforamtion

6000
5000
4000
Effort Normal[KN]
3000
2000
fc=40MPa
1000 fc=30MPa
fc=25MPa
0
-0.001 0 0.001 0.002 0.003 0.004
Déformation
Figure III.57: Confrontation des courbes Effort Normal -Déformation calculée

avec les différentes classes du béton
112
Plus la résistance à la compression du béton est importante plus l’effort normal ultime est
ultime est important.
Courbe d'intéraction N=f(M)

6000
fc=40MPa
fc=30MPa
5000 fc=25MPa
4000
Effort Normal [KN]
3000
2000
1000
0
-50 0 50 100 150 200 250
-1000
Moment[KN.m]
Figure III.58: Confrontation des courbes d’interactions calculées avec les

différentes classes du béton
La capacité portante devient importante avec l’accroissement de la résistance à la

n compression du béton.
III.1.4.6. Influence du diamètre de la section du béton

Le calcul est effectué pour une section circulaire de diamètre variable. Elle est armée de
barres d’acier de diamètre 12 mm et d’armatures transversales de diamètre 8mm. Les
caractéristiques des matériaux et les courbes obtenues, en utilisant les lois de comportement
réelles des matériaux, sont représentées sur le tableau III.13 et figure III.27 :
Tableau III.13: Caractéristiques du béton avec les différents diamètres de la section
Fcc(MPa) Ftj(MPa) Ebc0(MPa) Kc Kc’ rt cc0 ccu Dsec(cm)

33.26 2.56 35376 5.6 4.6 -0.002 0.005 0.05 30
31.31 2.48 34671 5.01 4.01 -0.002 0.0045 0.05 40
30.09 2.4 34217 3.6 2.6 -0.002 0.004 0.05 50
113
Courbe d'intéraction M=f(N)

7000
6000 D50
D40
5000 D30
Effort Normal [KN]
4000
3000
2000
1000
0
0 100 200 300 400
-1000
Moment [KN.m]
Figure III.59: Confrontation des courbes d’interaction avec les différents

diamètres de la section
La capacité portante s’amplifie proportionnellement avec l’accroissement du diamètre de

la section du béton.
III.2. Calcul global des poutres armées avec prise en compte du confinement
Afin de valider le modèle proposé, d’estimer la cohérence de son approche et de voir
l’apport du confinement sur le calcul des structures en béton armé confiné on a effectué le
calcul de plusieurs poutres en béton armé et précontraint qui ont fait l’objet de plusieurs
essais.
III.2.1 Essais du CEBTP – Poutres hyperstatiques HZ4
Il s’agit d’un programme d’essais réalisé au CEBTP (Trinh, 1995) [70]. Ce programme
comporte quatre poutres différenciées principalement par le degré de précontrainte qui varie
de 100 % (précontrainte totale) à 0 % (béton armé) et le pourcentage des armatures
transversales. Les quatre poutres sont identiques en forme et dimensions. Elles présentent une
section transversale en double T et une longueur totale de 10.40 m constituée de deux travées
de 5 m de portée chacune. Elles sont soumises à un chargement croissant jusqu’à rupture, par
114
deux charges concentrées d’intensité égale, Q, appliquées au milieu des travées. La figure
III.28 montre les caractéristiques géométriques et le détail de ferraillage de ces poutres.
Figure III.60: Dimensions et détails de ferraillage des poutres HZ [70].
La figure III.29 montre la disposition assez complexe des armatures passives

longitudinales dans les membrures inférieure et supérieure des différentes poutres HZ2, 3, 4.
La poutre HZ1, la plus précontrainte, est néanmoins munie d’un minimum d’armatures
passives, filantes sur toute la longueur de la poutre, en barres lisses  6 mm (en acier FeE 22).
Il n’y a pas d’armatures longitudinales dans l’âme. Chaque corps d’épreuve comporte une
armature transversale disposée en cours successifs dans des sections droites avec un
espacement t constant (voir tableau dans la figure III.28).
Figure III.61 : Disposition des armatures longitudinales pour les poutres HZ [70].
115
Les caractéristiques mécaniques des matériaux sont résumées dans tableau III.14.
Tableau III.14 : Caractéristiques du béton et des armatures passives [70].
Béton Armatures passives

Repère fcj(MPa) ftj(MPa) Eij(MPa) Barres fe Ea (GPa) fu (MPa)
HZ1 39 3.4 35400 6 (MPa)
340 206 435
HZ2 33 3.0 31250 HA10 428 198 545
HZ3 34 3.4 32080 HA16 430 213 526
HZ4 32 3.3 32000 HA20 424 195 543
HA25 450 230 558
La précontrainte des poutres est réalisée en post-tension. Les câbles, un par poutre, sont
constitués d’un nombre variable de fils 7 : 16, 10 et 7 respectivement pour les poutres HZ1,
HZ2 et HZ3. Le fil 7 est caractérisé par une limite d’élasticité de 1540 MPa, un module
élastique de 175 GPa et une contrainte de rupture de 1970 MPa.
Pour notre cas l’application du modèle avec prise en compte du confinement a été faite
uniquement pour la poutre HZ4 qui est en béton armé. La figure III.30 montre l’évolution de
la flèche, au droit de la charge, en fonction de la charge appliquée pour la poutre HZ4. Le
calcul non linéaire en flexion est effectué d’une part, en considérant comme rigidité à l’effort
tranchant la valeur de l’élasticité linéaire, et d’autre part dans un deuxième calcul, en
considérant la rigidité à l’effort tranchant et tout en tenant compte du confinement.
116
Charge Q(kN)
800
600
Valeures experimentales
400
Rigidité tangente solution de

l'élasticité linéaire
200
valeures calculées avec
modele proposé
0
0 5 10 15 20 25
Poutre HZ4 Flèche (mm)
Figure III.62 : Evolution de la flèche en fonction de la charge pour les poutres HZ.
Le tableau III.15 montre que les valeurs ultimes de la charge Q calculées par le présent
modèle sont proches de celles relevées expérimentalement.

Pour la poutre HZ4.
Valeur ultime de l’effort Qu kN

Calculée Calculée
Exp/présente
Poutre Expérimentale Kachi et Présente Exp/Cal kachi
étude
al.2006 étude
HZ4 590 570 580.023 1.03 1.01
On peut noter que le modèle développé estime de façon satisfaisante la contribution du

confinement à la charge ultime des poutres en béton armé. Il permet donc de bien simuler le
comportement, jusqu’à rupture, des poutres soumises à la flexion avec effort tranchant en
tenant compte des armatures transversales de confinement. Les efforts maximaux et les
déplacements à la rupture des poutres sont correctement estimés. Le logiciel développé
permet ainsi de prédire le mode de rupture réel des poutres, notamment dans le cas d’une
rupture à l’effort tranchant. En effet, dans un calcul classique des poutres, la rigidité de
cisaillement reste constante ; elle est solution de l’élasticité linéaire. Le calcul ne prédit alors
117
que les ruptures en flexion ; la résistance à l’effort tranchant doit être vérifiée
indépendamment par des formules empiriques.
III.2.2.Etude paramétrique
Le calcul a été fait pour la poutre HZ4 de l’essai réalisé au CEBTP [70].la poutres est en
béton armé dont les caractéristiques géométriques et le détail de ferraillage est résumé dans le
tableau III.16 :
Tableau III.16. Caractéristiques du béton et des armatures passives poutre HZ4 [70].
Béton Armatures logitudinales Armatures transversales

Repè fcj ftj Eij Barres fe Ea fu Barres fe Ea
re (MPa) (MPa (MPa) (MPa) (GPa) (MPa) (MPa (GPa
32 )
3.3 32000 HA20 424 195 543 HA10 )
424 )
195
HZ4
HA25 450 230 558 esp=20cm / /
Le calcul a été fait pour des espacements allant de 10 cm à 30 cm.

Le tableau III.17 montre que les valeurs ultimes de la charge Q et de la flèche ultime
calculée par le présent modèle pour les différents espacements.
Tableau III. 17: Comparaison des valeurs ultimes expérimentales et calculées de la charge Q
Espacement Charge ultime Flèche ultime Exp/présente

Poutre
Cm Qu kN (mm) étude
10 590.0253 17.26 0.99
13 590.0247 17.20 0.99
15 590.0246 17.00 0.99

HZ4
20 580.0232 16.37 1.017
25 525.0176 13.02 1.124
30 525.0176 13.02 1.124
On observe que la charge ultime et la flèche ultime augmente avec l’augmentation du

confinement dans le cas des poutres en béton armé .On remarque qu’il y’a un apport de
11,02% en terme de charge ultime pour un espacement de 10 à 30 cm et un apport de 24,56
% flèche ultime.
118
Charge Q (kN)
800
(590.0253, 17.26)
(590.0246, 17.00)
600
(590.0232, 16.37)
(525.0176, 13.02)
(590.0246, 17.00)
400
200
esp=10cm
esp=15cm
esp=20cm
esp=30cm
0
0 2 4 6 8 10 12 14 16 18 20
Poutre HZ4 Flèche (mm)
Figure III.63: Évolution de la flèche en fonction de la charge pour les poutres HZ4 pour des
espacements 10, 15,20 et 30cm
119
III.3.Calcul en cisaillement d’une section de poutre
Un programme de calcul en cisaillement-flexion des sections de poutres a été élaboré. Afin

de vérifier l’influence des simplifications admises dans cette étude ainsi que la cohérence de
son approche, le confinement de ces sections a été introduit dans le calcul Nous avons
effectué le calcul de plusieurs sections de poutre avec des formes et des détails de ferraillage
variés.
III.3.1. Essais de l’Université de Toronto –Comportement de section
Le modèle a été utilisé pour simuler le comportement de plusieurs poutres testées par
différents chercheurs à l’Université de Toronto (Vecchio Collins, [71], [63], [66]). Les
poutres ont une section transversale rectangulaire pleine ou creuse. Elles sont en béton armé
ou précontraint. La précontrainte est soit adhérente soit non adhérente (dans ce dernier cas,
elle est équivalente à une force axiale extérieure). Pour le calcul, les armatures passives sont
concentrées en deux lits inférieur et supérieur ; il n’y a pas d’armatures longitudinales
réparties. Différents modes de chargement ont étés appliqués ; cependant, dans tous les cas, le
moment fléchissant est nul dans la zone de poutre où se produit la rupture par effort tranchant
et où on a mesuré la distorsion. Les caractéristiques géométriques et les propriétés des
différents matériaux sont données au tableau (III.18).
120
Tableau III.18 : Propriétés des différentes poutres [71],[63],[66].
Dimensions Béton Armatures

Armatures Armatures de
Repèr transversales
longitudinale précontrainte
e de la s
poutre EXT INT fc b0 Barre st fe Nb.x fe Ap fpe p
Mm mm MPa 10-3 s  mm Mpa  Mpa mm2 Mpa 10-3
(m) (mm)
12x 29 345 _ _ _
SA3 305x61 152x40 40.0 2.8 9.5 72 373 4 x 22 462
0 6
12x 29 345 _ _ _
SA4 305x61 152x40 40.0 2.8 9.5 72 373 4 x 22 462
0 6
SK1 305x61 - 26.9 2.2 9.5 100 400 8 x 25 442 154 1450 4.8
0 5 0 2
SK2 305x61 121x38 26.9 2.2 9.5 100 400 8 x 25 442 154 1450 4.8
0 1 5 0 2
SK3 305x61 - 28.2 2.2 9.5 100 400 8 x 25 442 _ _ _
0
SK4 305x61 121x38 28.2 2.2 9.5 100 400 16x 25 442 - - -
0 1
SP0 305x61 152x40 25.0 2.3 9.5 150 373 16x 25 421 _ _ _
0 6
SP1 305x61 152x40 33.5 2.3 9.5 150 373 12x 22 421 510 1450 4.2
0 6 1
SP2 305x61 152x40 32.0 2.0 9.5 150 373 12x 22 421 101 1450 4.1
0 6 0 1
SP3 305x61 152x40 32.2 2.0 9.5 150 373 12x 22 421 152 1450 4.2
0 6 0 6
SM1 305x61 152x40 29.0 2.4 9.5 175 424 12x 22 452 -_ -_ -_
0 6
CF1 305x61 152x40 38.6 3.0 9.5 150 367 6 x 9.5 367 930 1450 5.1
0 6 7
Les résultats obtenus pour les poutres SA3, SA4, SK1, SK2, SK3 sont présentés à la figure
III.32 sous forme de courbes comparatives avec les résultats expérimentaux d’une part, et les
résultats obtenus par Vecchio et Collins [71], [63], [66] ainsi que ceux obtenus par Kachi et
al. [52-60]. d’autre part. Il faut bien noter qu’il s’agit de la relation effort tranchant V –
distorsion moy d’une zone de poutre soumise au cisaillement (éventuellement combiné à un
effort normal) sans gradient des déformations longitudinales et en tenant compte du
confinement.
121
V(kN)
V(kN)
900
800
700
600
500
Experimentale
400
Collins ET Vecchio
300 Kachi2006
Kachi et al. 2006
200 Présente étude
100
0
0 2 4 6 8
Poutre SA3 moy (10-3)
V(kN)
600

500
400 V(kN)
300
Experimentalle
Collins et Vecchio
200
Kachi et al. 2006
Kachi 2006
100 Présente étude
0
0 2 4 6 8
122
V(kN)
900
800
700
600
500
experimentale
400
Collins et Vecchio
300
Kachi2006
Kachi et al. 2006
200
présente étude
100
0
0 2 4 6 8 10 12
Poutre SK1 moy (10-3)
V(kN)
800
700
600
500
400
Experimentalle
300 Collins et Vecchio
200 Kachi et al. 2006
Kachi2006
100 Presente étude
0
0 2 4 6 8
PoutreSK2 moy (10-3)
123
V(kN)
900
800
700
600
500
Kachi et al. 2006
Kachi 2006
400 Experimentalle
Collins et Vecchio
300
Presente étude
200
100
0
0 2 4 6 8 10 12
PoutreSK3 moy (10-3)
Figure III.64 : Comparaison des résultats du calcul aux résultats expérimentaux et à ceux
obtenus par le modèle de Vecchio et Collins et le modèle de Kachi et al.2006.
Le tableau III.19 compare les valeurs ultimes de l’effort tranchant calculées par le présent
modèle à celles relevées expérimentalement.
Tableau III. 19: Comparaison des valeurs ultimes expérimentales et calculées de l’effort
tranchant.
Valeur ultime de l’effort tranchant Vu Kn
Poutre Expérimentale Calculée Exp/Cal
SA3 716 771 0.93
SA4 534 569.75 0.94
SK1 672 818.5 0.82
SK2 530 609.75 0.87
SK3 725 792.25 0.91
Moyenne des rapports (Exp/Cal) de l’ensemble des poutres 0,894
124
On peut noter que les résultats du calcul, comparés aux résultats expérimentaux,
montrent que le modèle approche bien le comportement réel de l’ensemble des poutres
étudiées, aussi bien avant qu’après la fissuration du béton. L’erreur commise sur les valeurs
ultimes de la charge varie entre 0 et 18 %. Le calcul prévoit assez correctement l’évolution de
moy avec V, donc la rigidité à l’effort tranchant en tenant compte du confinement.
Le terme de cisaillement (rigidité ∆v/∆moy) est pratiquement le même pour le béton
confiné et non confiné. Mais une fois la limite des aciers est atteinte, on remarque un léger
accroissement de la rigidité au cisaillement dans le cas du béton confiné. La déformation et
l’effort à la rupture augmentent avec le confinement.
Dans le cas du béton précontraint on remarque que le confinement augmente la rigidité
au cisaillement correspondante à la déformation et l’effort à la rupture.
III.4.Etude de la fissuration
Le calcul a été fait pour la poutre HZ4 de l’essai réalisé au CEBTP [70].la poutre est en béton
armé dont les caractéristiques sont donnés dans figure III.33 et le détail de ferraillage est
donné dans le tableau III.16.
Figure III. 65 : caractéristiques géométriques et détails de ferraillage de la poutre HZ4
Les figures III.34, III.35, III36 montrent la confrontation des résultats obtenus de la hauteur
décomprimée en fonction de la déformation ex ; avant la fissuration, au début de la fissure et
à la rupture pour la poutre HZ4 par la présente étude en tenant compte du confinement et les
résultats obtenus par Kachi et al. En 2006 pour la poutre HZ4 sans tenir compte du
confinement.
125
Y (m)
0.7
0.6
0.5
0.4
0.3
0.2 Kachi 2006 (non conﬁné)

Présente étude (confiné)
0.1
0
-1.00E-03 -5.00E-04 0.00E+00 5.00E-04 1.00E-03 1.50E-03 x

section avant fissuration de la poutre HZ4.
Y (m)
x

section au début de la fissuration de la poutre HZ4.
126
Y (m)
0.7
0.6
0.5
0.4
0.3
0.2 Kachi 2006 (non conﬁné)

Présente étude (confiné)
0.1
0
-2.50E-03 -1.50E-03 -5.00E-04 5.00E-04 1.50E-03 x

section à la rupture de la poutre HZ4.
On remarque bien d’après ces figures que le confinement retarde la fissuration et la

rupture de la poutre. La hauteur comprimée à la rupture est de 17cm sans confinement et de
5.1cm en tenant compte du confinement par conséquent la hauteur décomprimée est de 43 cm
sans confinement et de 54.9 cm en tenant compte du confinement
En faisant varier l’espacement de 5 cm à 25 cm. Les résultats de calcul pour la
Distribution des déformations longitudinales moyennes sur la hauteur de la section à la
rupture de la poutre HZ4 montre que la hauteur décomprimé est inversement proportionnelle
au pourcentage de confinement. La figure III.37 montre la variation de la hauteur
décomprimée en fonction de l’espacement :
127
Y (m)
hauteur décomprimée en fonction de l'espacement

0.60
h=f(esp)
0.58
0.55
0.53
0.50
0.00 5.00 10.00 15.00 20.00 25.00 30.00 esp (cm)
Figure III.69: hauteur décomprimée de la section en fonction de l’espacement a la rupture de

la poutre HZ4.
Il y a une variation très forte de la hauteur décomprimée et ce jusqu’à un espacement de

10 cm, au-delà de cette valeur, la hauteur décomprimée varie très peu. Le confinement se fait
sentir pour des valeurs d’espacement inférieures à 10 cm.
En faisant varier l’espacement de 7 cm à 25 cm. Les résultats de calcul pour les efforts
internes effort tranchant, Moment de flexion et effort normal au début des premières fissures
et à la rupture de la poutre HZ4 sont montés dans le tableau III.20:
Tableau III.20 : Résultats calculés des efforts au début des premières fissures de la Poutre
HZ4.
Espacement Effort tranchant Moment Effort normal

Poutre
Cm T kN M kN.m N x10-2 kN
HZ4 7 155.28 201.43 6.00
HZ4 10 154.68 198.42 6.61
HZ4 15 153.95 194.76 7.30
HZ4 17 153.71 193.57 7.30
HZ4 20 153.34 191.68 7.74
HZ4 25 153.25 191.22 7.99
128
Tableau III.21 : Résultats calculés des efforts à la rupture de la Poutre HZ4.
Espacement Effort tranchant Moment Effort normal

Poutre
Cm T kN M kN.m N x10-1 kN
HZ4 7 375.66 426.65 8.06
HZ4 10 368.88 419.25 8.09
HZ4 15 368.63 417.96 8.11
HZ4 17 368.55 417.58 8.19
HZ4 20 368.45 417.09 8.27
HZ4 25 335.11 316.90 8.99
On remarque que les efforts T et M sont proportionnel au pourcentage de confinement .on

voit que le confinement retarde légèrement la fissuration et la rupture.
III.5.Pieu en béton armé soumis à une flexion composée
Nous avons effectué le calcul avec le programme sur un pieu dont l’étude
expérimentale a été réalisée au CEBTP [72].
Le pieu à un diamètre de 500 mm et de 4m de long .Il est soumis à une flexion
composée. L’effort normal de compression est de 1370 kN (appliqué à l’aide d’une
précontrainte extérieure). Le pieu est renforcé longitudinalement par 5 barres d’acier HA de
diamètre 16mm (correspondant à 0,5 % en volume) selon le min du DTU 13.2 et
transversalement par des HA 6.les caractéristiques géométriques et le détail de ferraillage est
résumé dans le tableau(III.22) suivant :
Tableau III.22 : Caractéristiques du béton et des armatures passives [72].
Béton Armatures logitudinales Armatures transversales

Repè fcj ftj Eij Barre fe Ea fu barres fe Ea
re (MPa) (MPa (MPa) s (MPa (GPa) (MPa) (MPa (GPa
Pieu ) 5HA ) ) )
42.41 4.35 38059.47 475 20 475 HA6 428 200
zhan 16
La confrontation du calcul par la présente étude aux résultats expérimentaux en terme de charge
flèche est illustrée à la Figure III.38.
129
Charge Q(kN)
4.50E+02
4.00E+02
3.50E+02
3.00E+02
2.50E+02
2.00E+02 Zhan exp
1.50E+02 Présente étude
1.00E+02
5.00E+01
0.00E+00
0.00E+00 2.00E+00 4.00E+00 6.00E+00 8.00E+00 1.00E+01 1.20E+01 1.40E+01
-5.00E+01
Pieu Zhan Flèche (mm)
Figure III.70 : Évolution de la flèche en fonction de la charge pour le pieu de Zhan.
Le tableau III.23 montre que les valeurs ultimes de la charge Q calculées par le présent
modèle sont proches de celles relevées expérimentalement soit une erreur de 2%.

Au début de la première fissure et à la rupture.
Valeur ultime de l’effort Qu kN

Présente
Pieu Expérimentale Exp/présente étude
étude
Apparition
200 200.11 0.99
des 1 ères fissures
A la rupture 410 417 0.98
On peut noter que le modèle développé estime de façon satisfaisante la contribution du

confinement. Il permet donc de bien simuler le comportement, jusqu’à rupture, des poutres
soumises à la flexion composée avec effort tranchant en tenant compte des armatures
transversales de confinement (soit 98% de précision).
130
III.6. Conclusion
La méthode de calcul présentée dans le cadre de ce travail a pour objectif l’introduction

du confinement en tenant compte des différentes dispositions des armatures transversales ou
des tubes, ainsi que l’effet de l’effort tranchant dans le domaine élastique non linéaire dans un
logiciel de calcul des poutres hyperstatiques. Le confinement est introduit par une loi de
comportement également non linéaire qui permet de mieux simuler le comportement des
sections en béton armé confiné.
L’effet de l’effort tranchant dans le domaine non linéaire a été pris en compte dans le
modèle de calcul par une estimation du module de cisaillement réel d’une section de poutre.
La simulation du comportement des poutres hyperstatiques en béton armé avec prise en
compte du confinement, a été faite jusqu’à la rupture.

satisfaisante.
Le modèle développé estime de façon satisfaisante la contribution du confinement et

de l’effort tranchant à la déformabilité des poutres. Il permet donc de bien simuler le
comportement, jusqu’à rupture, des poutres soumises à la flexion composée avec effort
tranchant en tenant compte des armatures transversales de confinement.
Les résultats du calcul, comparés aux résultats expérimentaux, montrent que le
modèle approche bien le comportement réel de l’ensemble des poutres étudiées, aussi bien
avant qu’après la fissuration du béton. Le calcul prévoit assez correctement l’évolution de
moy avec V, donc la rigidité à l’effort tranchant en tenant compte du confinement.
131
Le confinement augmente la rigidité au cisaillement correspondante à la déformation et

l’effort à la rupture.
La hauteur décomprimée a une variation très forte et ce jusqu’à un espacement de 10
cm, au-delà de cette valeur, la hauteur décomprimée varie très peu. Le confinement se fait
Le modèle développé estime de façon satisfaisante la contribution du confinement. Il
permet donc de bien simuler le comportement, jusqu’à rupture, des poutres soumises à la
flexion composée avec effort tranchant en tenant compte des armatures transversales de
confinement.
132
CONCLUSION GENERALE
La prise en compte du confinement dans la modélisation et le calcul des structures en

béton armé et en béton précontraint à donner des résultats satisfaisants et à permis de mieux
approcher le comportement réel.
Vue la complexité des modèles existants et la difficulté de maîtriser certain de leur
paramètres, nous avons proposé dans le cadre de cette étude un modèle simple dont les
paramètres sont facilement maîtrisable.
Le confinement est introduit par un coefficient de confinement effectif, le pourcentage
d’armatures transversales et la limite élastique de l’acier dans le cas d’un béton confiné par
armatures transversales ainsi que l’épaisseur du tube dans le cas d’un béton confiné par
CFST.
L’effet de l’effort tranchant dans le domaine non linéaire a été pris en compte dans le
modèle de calcul par une estimation du module de cisaillement réel d’une section de poutre.
La simulation du comportement des poutres hyperstatiques en béton armé avec prise en
compte du confinement, a été faite jusqu’à la rupture.

satisfaisante.
Le modèle développé estime de façon satisfaisante la contribution du confinement et de
l’effort tranchant à la déformabilité des poutres. Il permet donc de bien simuler le
comportement, jusqu’à rupture, des poutres soumises à la flexion composée avec effort
tranchant en tenant compte des armatures transversales de confinement.
Les résultats du calcul, comparés aux résultats expérimentaux, montrent que le modèle
approche bien le comportement réel de l’ensemble des poutres étudiées, aussi bien avant
133
qu’après la fissuration du béton. Le calcul prévoit assez correctement l’évolution de moy
avec V, donc la rigidité à l’effort tranchant en tenant compte du confinement.
Le confinement augmente la rigidité au cisaillement ainsi que la déformation et l’effort

à la rupture.
La hauteur décomprimée a une variation très forte et ce jusqu’à un espacement de 10
cm, au-delà de cette valeur, la hauteur décomprimée varie très peu. Le confinement se fait
Le modèle développé estime de façon satisfaisante la contribution du confinement. Il
permet donc de bien simuler le comportement, jusqu’à rupture, des poutres soumises à la
flexion composée avec effort tranchant en tenant compte des armatures transversales de
confinement.
Le modèle analytique présenté repose sur la modélisation multicouche des sections.
L'équilibre global des sections est analysé sur la base de la planéité des sections dans le sens
longitudinale. Le modèle introduit un terme de cisaillement diagonal non linéaire dans la
matrice de rigidité de la section pour prédire le comportement global jusqu'à rupture des
poutres en flexion composée. Les résultats de la confrontation sont très satisfaisants.
Le confinement interne du béton est introduit et des confrontations ont été faites sur
des sections de poutre ainsi que sur des poutres en béton armé et en béton précontraint, les
résultats du calcul confrontés aux résultats expérimentaux sont très satisfaisants.
Dans le cadre de cette étude nous proposons une méthode de calcul des éléments de
section circulaire basée sur la méthode de calcul par les armatures discrètes est développée en
y introduisant différentes formes de confinement internes et externe. L’étude est menée en
élasticité non linéaire, il est tenu compte du comportement non linéaire des différents
matériaux constituant la section transversale. Le confinement du béton est introduit par le
biais d’une loi de comportement également non linéaire, du béton confiné et du béton armé
confiné tenant compte du niveau de confinement, du type de confinement et de l’ensemble
des sollicitations (effort normal, moment fléchissant et effort tranchant). Enfin, la méthode de
calcul est injectée dans un logiciel de calcul non linéaire des éléments de structures afin
d’étudier la contribution et l’influence du confinement du béton sur le comportement global
des éléments de structures de section quelconques. Les résultats de la confrontation sont très
satisfaisants.
134
REFERENCES BIBLIOGRAPHIQUES
[1] Richart, F.E., Brandtzaeg, A.,and Broun,R.L. ,« The failure of plain and spirally
reinforced concret in compression. » Bulletin 190, 1929, Univ.of Illinois, Engineering
Experimental Station, Champaing III.
[2] Belounis, R.et al., « Etude de la performance des poteaux Mixtes acier-beton sous
differentes Sollicitations », mémoire de magister, Université Mentouri Constantine, 2008.
[3] Y.C.Wang, «Tests on slender composite columns», Journal of Constructional Steel

Research 49 (1999) 25–41.
[4] Kefeng.T., John.M.N., Xincheng.P., « Mechanical properties of high strength concrete

filled steel tubular columns Part1: Concentrically loaded », ACI Journal 22 April 2003.
[5] M.Mouli, H.Khelafi, «Strength of short composite rectangular hollow section columns
filled with lightweight aggregate concrete », Engineering Structures (2006),
doi:10.1016/j.engstruct.2006.10.00.
[6] Dalin,Liu, CTests on high-strength rectangular concrete-filled steelhollow section stub

columns»,Journal of Constructional Steel. Research 61 (2005) 902–911.
[7] Ehab Ellobody,Ben Young, Dennis Lam, « Behaviour of normaland high strength
concrete-filled compact steel tube circular stub columns », Journal of Constructional Steel
Research 62 (2006) 706–715.
[8] M.Mouli, H.Khelafi, «Strength of short composite rectangular hollow section columns
filled with lightweight aggregate concrete », Engineering Structures (2006),
doi:10.1016/j.engstruct.2006.10.00.
[9] Dalin Liu., «Tests on high-strength rectangular concrete-filled steel hollow section stub
columns»,
Journal of Constructional Steel, Research 61 (2005) 902–911.
135
[10] Muhammad Naseem Baig, FAN Jiansheng, NIE Jianguo, «Strength of Concrete Filled
Steel Tubular Columns», ISSN 1007-0214 05/15 pp657-666 Volume 11, Number 6,
December 2006.
[11] Ever.J.Barbero, «Prediction of buckling mode interaction in composite columns», ACI

Journal December 10, 1999.
[12] A.A. Marinopoulou,V.D. Balopoulos, C.N. Kalfas, «Simulation of partially encased

composite steel concrete columns with steel columns»,Journal of Constructional Steel
Research (2006),doi:10.1016/j.jcsr.2006.11.003.
[13] S.Ramana Gopal, P. Devadas Manoharan,«Expérimental behaviour of eccentrically

loaded slender circular hollow steel columns in-filled with fibre reinforced concrete»,Journal
of Constructional Steel Research 62 (2006) 513–520.
[14] Eurocode 4 – Conception ET Dimensionnement des Structures Mixtes Acier –Béton.

Partie 1-1: règles générales et règles pour les bâtiments.
[15] Mohanad Mursi, Brian Uy “Strength of slender concrete filled high strength steel box
columns”. Journal of Constructional Steel Research 60 (2004) 1825–1848
[16] Richard Maranda “Analyse par éléments finis de poteaux mixtes avec section d’acier en I
de classe 4”.Mémoire présentée en vue de l’obtention du diplôme de maîtrise des sciences
appliquées (génie civil) décembre1998. Département des génies civil, géologique, et des
mines, Ecole polytechnique de MONTRÉAL.
[17] Mathias Johansson “Composite Action and Confinement Effects in Tubular Steel-
Concrete Columns”. Steel and Composite Structures.
[18]Hognestad E. “A study of combined bending and axial load in reinforced concrete

members”. Bulletin series no. 399, Engineering Experiment Station,Univ. of Illinois, Urbana
Ill; 1951.
[19] Bousalem B., Chikh N., « Development of a Confined Model for Rectangular Ordinary
Reinforced Concrete Columns » ;materials and structures, 2007, vol. 40, pp 605-613.
136
[20] Ahmad, S. M., and Shah, S. P. (1982). “Complete triaxial stress-strain curves of concrete
confined by spiral reinforcement,” J. Struct.Div. ASCE, 108(4), 728-742.
[21] Ahmad, S. M., and Shah, S. P. (1982). “Stress-strain curves of concrete confined by
spiral reinforcement,” ACI Struct. J., 79(6),484-490.
[22] Toutanji, Houssam A., (1999), “Stress-Strain Characteristics of concrete Columns

Externally Confined with Advanced Fiber Composite Sheets,” ACI Material Journal, V.96,
No. 3, May-June 1999, pp. 397-404.
[23] Mander, J. B., Priestley, M. J. N., and Park, R., Observed Stress-Strain Behavior of
Confined Concrete.” J. Struct. Engrg., ASCE, 114(8), 1827-1849,1988.
[24] Mander, J. B., Priestley, M. J. N., and Park, R., “Theoretical Stress-Strain Model for
Confined Concrete.” J. Struct. Engrg., ASCE, 114(8), 1804-1826,1988.
[25] saatcioglu M. and Razvi S.R., « Strength and Ductility of Confined Concrete », journal
of structural engineering , 1992, vol .118,No.6 pp 1590-1607.
[26] Hoshikuma J,K . Kawashima,K. Nayaga and A.W. Taylor, « Strass-Strain Model for
Confined RC in bridge Piers », Journal of structural engineering, 1994, vol .123,No.0.5 pp
624-633.
[27] Fardis, M.N., and Khalili, H. (1981) “Concrete encased in fiberglass-reinforced plastic.”
ACI J., 78(6), 44-446.
[28] Karbhari, V. M., and Gao, Y. (1997), “Composite jacketed concrete under uniaxial
compression –Verification of simple design equations.” J. Mater. Civ. Eng.,9(4), 185-193.
[29] Mirmiran, A., and Shahawy, M., (1997), “Behavior of Concrete Columns Confined by
Fiber Composites,” ASCE Journal of Structural Engineering, Vol. 123,No. 5, 583-590.
137
[30] Wu, H., and Xiao, Y., (2000), “Compressive Stress Strain Behavior of Concrete
Comfined by Carbon Fiber Jackets,” Proceedings of the Sixth ASCCS Conference on
Composite and Hybrid Structures, Y. Xiao, and S. A. Mahin,eds., Los Angeles, Calif., Mar.
22-24, 2000. pp.
[31] Harmon, T. G.; Ramakrishnan, S. and Wang, E. H., (1998), “Confined Concrete
Subjected to Uniaxial Monotonic Loading,” ASCE Journal of Structural Engineering, Vol.
124, No. 12, December, 1303-1309.
[32] Fujikake, Kazunori; Mindess, Sidney and Xu, Hanfeng (2004), “Analytical Model for
Concrete Confined with Fiber Reinforced Polymer Composite,” J.Compos. for Constr.,
Volume 8, Issue 4, pp. 341-351 July/August 2004.
[33] Becque, J., Patnaik, A. K., and Rizkalla, S. H. (2003), “Analytical models for concrete
confined with FRP tubes,” J. compos. Constr., 7(1), pp 31-38.
[34] Bouafia Y, Idir A, Kachi MS., « Influence of the taking account of the confined concrete
on the structures non linear calculation », ACMA 2010 May19-21 2010. Marrakech,
Morocco.
[35] Iddir. A., « influence de la prise en compte du béton confiné dans le calcul non lineaire
des sctions»., mémoire de magister, université mouloud mammeri, tizi ouzou, 2009.
[36] Iddir, A., Bouafia,Y., Kachi. Ms., « La prise en compte du confinement dans le calcul
non linéaire des sections», CMMS 08 Les 16, 17 & 18 November 2008, Université M. M. de
Tizi-Ouzou, Algérie.
[37] Bouafia Y, Idir A, Kachi MS.,Dumontet H. « stress- strain relationship for the confined
concrete », WCCM XI-ECCM V-ECFD VI, July 20 - 25, 2014, Barcelona, Spain.
[38] Sargin M.,« Stress –Strain Relationship for concrete and the analysis of structural
concrete sections » PHD thesis, University of Waterloo, Ontario, Canadi; Mar; 334 p,1968.
138
[39] Virlogeux M., « la précontrainte extérieure, Innovation dans le domaine des ouvrages
d’arts en béton précontraint », Annales de l’ITBTP, Décembre, 1983.
[40] Virlogeux M., « étude d’une section de poutre en élasticité non linéaire », Annales de
l’ITBTP, Mai, 1986.
[41] Virlogeux M., « Calcul des structures en élasticité non linéaire », Annales des ponts et
chaussées n°39 et 40, 1986.
[42] Virlogeux M., M’rad, A., « étude d’une section de poutre en élasticité non linéaire »,
Annales de l’ITBTP, n°444, Mai, 1986.
[43] Virlogeux M., « Calcul des structures en élasticité non linéaire », Annales des ponts et
chaussées, 3ème trim.1986.
[44] Virlogeux M., « la précontrainte extérieure au béton », X congrès International de la FIP,

New delhi, 1986.
[45] Virlogeux M. M’rad A., « Flexural behavior of externally prestressed structures for
ultimate loads, workshop behavior of external prestressing in structures », Rémy –lès –
Chevreuses, juin 1993.
[46] Virlogeux M., « External prestressing : from construction history to modern technique
and technology, External prestressing in bridges », ACI, SP-120, Ed Antoine Naamane and
John Breen, 1990.
[47] M’rad A., « Calcul à la rupture des ponts à précontrainte extérieure au béton »,Thèse de
doctorat de 3ème cycle, université pierre et Marie curie (Paris VI), 1984.
[48] Muller J., Gautier Y., « Ultimate béhavior of precast segmental box girders with external
tendons », External prestressing in bridges, ACI, SP-120, Ed Antoine Naamane and John
Breen, 1990.
[49] Kihg-Hwee.,Chee Khoon NG., « effect of behavior and tendon configuration on behavior
of externally prestressed beams »,ACI structural journal. V.94 N°1, janvier, fevrier 1997.
139
[50] Rezendi-Martin P.L., « Modélisation du comportement jusqu’à la rupture en flexion des
poutres en béton à précontrainte extérieure ou mixte », Thèse de doctorat à l’école centrale de
paris, 1989.
[51] Guriani E., « on the effective axial stiffness of a bar in cracked concrete Bond in
concrete », Ed. P. Bartos, Applied science publishers, Londres, 1982.
[52] Kachi M.S., « modélisation du comportement jusqu’à rupture des poutres à précontrainte
extérieure», Thèse de doctorat, Université de Tizi - ouzou, Algérie, février 2006.
[53] Nait-Rabah O., « Simulation numérique des ossatures spatiales », Thèse de doctorat,
Ecole Centrale de Paris, 1990.
[54] Garcia J.F., « Elaboration d’une méthode d’évaluation de la vulnérabilité sismique des
bâtiments », INSA Strasbourg, Spécialité- Juillet 2007.
[55] Kachi, MS., Fouré, B., Bouafia, Y., Muller, P., « L'effort tranchant dans la modélisation
du comportement jusqu'à rupture des poutres en béton armé et précontraint », Revue
Européenne de Génie Civil Volume 10, Décembre 2006, pages 1235-1264.
[56] Kachi, M S., Bouafia Y., Muller P., Fouré B ., « Modélisation du comportement jusqu’à
rupture des poutres en béton armé et précontraint en élasticité non linéaire », Colloque
international Caractérisation et Modélisation des matériaux et structures CMMS08 » Tizi
Ouzou – Algérie. Le 16, 17 et 18 Nov 2008.
[57] Kachi, MS., Bouafia Y., Muller, P., « Modélisation du comportement jusqu’à rupture en
cisaillement d’une section de poutre en béton armé et précontraint », Conférence
Internationale Maitrise et Management des Risques Industriels (M2RI’08) Oujda 24-25 avril
2008.
[58] Kachi, MS., Bouafia Y., Muller, P., « Modelization of Prestressed or Reinforced
Concrete beams until fracture in Non Linear Elasticity », 12th International Conference on
Fracture (ICF12). Ottawa, Canada, July 12-17, 2009.
140
[59] Bouafia, Y., Kachi, M S., Muller, P., « Model ling of Externally Prestressed Beams
Until Fracture on Non Linear Elasticity », International Conference on Structural analysis of
Advanced Materials. (IC SAAM-2009) Ecole Nationale D’Ingénieurs de Tarbes, France
September 7-10, 2009.
[60] Bouafia, Y., Fouré, B., Kachi, MS., « Shear strength of externally prestressed beams ».
II. International Symposium on Cement and Concrete Technology in the 2000s.Turkish
Cement Manufacturers’ Association with European Cement Association.Sept. 6 - 10, 2000,
Istanbul - Turkiye, ISBN 975-8136-09-7. Vol. 2, pp. 522-531.
[61] Grelat A., « Analyse non linéaire des ossatures hyperstatiques en béton armé », Thèse de
docteur Ingénieur, Université Paris VI, 1978.
[62] BPEL rules, « Technical design and costing and prestressed concrete structures
according to limit state method », DTU-BPEL Rules 91, Revised 99.
[63] Vecchio F.J., Collins M.P., « The modified compression field theory for reinforced
concrete elements subjected to shear », ACI Journal, March-April 1986, p. 219-231.
[64] Belarbi A; Thomas T.C.HSU, « Constitutive low of softened concrete in biaxial

tension-compression », ACI Structural Journal September-October 1995.
[65] Foure, B., «Déformations limites des armatures tendues et du béton comprimé pour le
calcul sismique des stuctures », VI Colloque National AFPS,juillet 2003,vol II ,pp 67-74.
[66] Vecchio F.J., Collins M.P., « Predicting the response of reinforced concrete beams
subjected to shear using modified compression field theory », ACI Structural Journal, May-
June 1988, p. 258-268.
[67] Lahlou K, Lachemi M, Aitcin PC., « Behavior of HSC filled tube columns under
dynamic compressive loading », ASCE, Journal of Structural Engineering
1999;125(10):1100–8.
141
[68] Bouafia, Y., kachi, M.S., Ait tahar, K., Ghazi, F., « Prédiction de la rupture des sections
circulaires en béton armé, et béton armé de fibres », JET 2006, association franco -
maghrébine de mécanique et des matériaux, 26-28 Avril 2006, Marrakech, Maroc.
[69] Ghazi F., « Evaluation de la Resistance à la Rupture en Flexion Composée des Sections
Circulaire en Béton Armé et en Béton Armé de Fibres », Mémoire de magistère, Université de
Tizi - ouzou, Algérie, 2005.
[70] Thrin J.L., « Précontrainte partielle essais de poutres continues »,Annales de

l’I.T.B.T.P.N°530,série béton 318,janvier 1995.
[71] Vecchio F.J., Collins M.P., « The response of reinforced concrete to in-plane shear and
normal stresses », University of Toronto, Department of civil Engineering, Publication No.82-
03, March 1982.
[72] Zhan Z.F., Foure B. Trinh L.J., « contribution au dimensionnement des pieux en béton de
fibres » .Doctorat, Université D’orleans C .E.B.T.P 1991.
[73] Kachi, Ms., Bouafia, Y., Saad, M., Bouhrat, S., « Contribution au calcul des sections
circulaires en béton armé par les armatures discrètes », ACMA 2012, May09-12 2012. Fès,
Morocco.
[74] Dall’Asta A., Zona A., « Three field mixed formulation for the non linear analysis of
composites beams with deformable shear connection », Finit element in analysis and design.
40. (2004).PP425-448.
[75] Dehousse N., « Fiber reinforced cement and concret », Ryhlem symposium, G.B.,
communication 4.1. The construction press L.T.D., ED. Adam Neuilly, 1975.
[76] Di Prisco M., Gombrova G., « compressive model for study of shear in thin-webbed RC
and PC beams », Journal of structural engineering. December 1995.
[77] Drizi S., « contribution à l’étude du comportement des voiles porteurs en béton vis-à-vis
des sollicitations sismiques », Thèse de doctorat, ParisVI, avril 1993.
142
[78] Dronnadulla V. and all., « Ultimate load behaviour of precast single-cell segmental box
bridges with externally post-tensioning », Workshop behavior of external prestressing in
structures.Saint Rémy-lés-Chevreuse, juin 1993.
[79] Edington J., Hannant D.J., « Stel fiber reinforced concrete, the effect of fiber oriontation
of compacting by vibration », Matériaux et construction 5, n°25, pp41-44, janvier –fevrier
1972.
[80] Espion B., « contribution à l’analyse non linéaire des ossatures planes, Application aux
structures planes en béton armé », Thèse de doctorat, Université libre de Bruxelles, 1986.
[81] Falkner H., Teutsch M., Huang Z., « prestressed segmental box girders winth unbonded
tendons under combined torsion bending and shear », Workshop behavior of external
prestressing in structures.Saint Rémy-lés-Chevreuse, juin 1993.
[82] Ferraro Maia A.C, Grelat A., Foure B., « Analyse non linéaire des ossatures en béton
armé ou précontraint compte tenu de retrait, du fluage et de la relaxation », Annales de
l’I.T.B.T.P. N°412, Février 1983.
[83] Foure B., « déformations limites des armatures tendues et du béton comprimé pour le
calcul sismique des structures », VI colloque national AFPS. juillet 2000, vol. II,pp 67-74.
[84] Foure B., « Béton à haute résistance », rapport interne OG, C.E.B.T.P., Nov.1985.
[85] Foure B., Rezende-Martins P.C. DE., Howang L.M., « Problème de sécurité à rupture et
de modelisation du comportement des poutres en béton à précontrainte extérieure », Annales
de l’I.T.B.T.P. Février 1991.
[86] Zhong.T, Lin.Hai.H, Zhi-Bin.Wang, « Experimental behaviour of stiffened concrete

filled thin walled hollow structural (HSS) stub column », Journal of Constructional Steel.
Research 61 (2005) 962–983.
[87] Mohanad Mursi, Brian Uy., « Strength of slender concrete filled high strength steel box
column », Journal of Constructional Steel Research 60 (2004) 1825–1848.
143
[88] Muhammad Naseem Baig, FAN Jiansheng, NIE Jianguo, «Strength of Concrete Filled
Steel Tubular Columns »,ISSN 1007-0214 05/15 pp657-666 Volume 11, Number 6,
December2006.
[89] Ever.J.Barbero, «Prediction of buckling mode interaction incomposite columns »,ACI

Journal December 10, 1999.
[90] A.A. Marinopoulou,V.D. Balopoulos, C.N. Kalfas, « Simulation of partially encased

composite steel–concrete columns with steel columns », Journal of Constructional Steel
Research (2006)
[91] N.E.Shanmugan, B.Lakshmi, «State of art report on steel –concrete composite

columns », Journal of Constructional Steel. Research 57 (2001) 1041–1080.
[92] Lin-Hai Han, You-Fu Yang, «Cyclic performance of concretefilled steel CHS columns
under flexural loading », Journal of Constructional Steel. Research 61 (2005) 423–452.
[93] Richard Maranda, «Analyse par éléments finis de poteaux mixtes avec section d’acier en
I de classe 4», Mémoire présentée en vue de l’obtention du diplôme de maîtrise des sciences
appliquées (génie civil) décembre1998. Département des génies civil, géologique, et des
mines, Ecole polytechnique deMONTRÉAL.
[94] Johansson M. and Gylltoft K., «Mechanical Behaviour of Circular Steel-Concrete

Composite Stub Columns», Journal of Structural Engineering, Vol. 128, No. 8, August
2002,ASCE, 1073-1081.
[95] C.C. Weng, S.I. Yen, « Comparisons of concrete-encased composite column strength
provisions of ACI code and AISC specification», Engineering Structures 24 (2002) 59–72.
[96] Kachi M.S., « Calcul non linéaire, jusqu’à rupture, d’une section en béton armé de fibres
métalliques », Thèse de magister, Université de Tizi - ouzou, Algérie, février 1997.
144
[97] Bouafia, Y., Kachi, MS., Fouré, B., « Relation contrainte - déformation dans le cas du
béton armé de fibres d'acier », Annales de l'ITB, Ed. ESKA, France, N°3, juin 2002, pp. 5-
17.
[98] Bouafia, Y., Kachi, MS., Fouré, B., « Numerical modeling of the behavior of steel fiber
reinforced concrete », II. International Symposium on Cement and Concrete Technology in
the 2000s. Turkish Cement Manufacturers’ Association – European Cement Association.
Sept. 6 - 10, 2000, Istanbul - Turkiye, ISBN 975-8136-09-7. Vol. 2, pp.582-591.
[99] Bouafia, Y., Fouré, B., Kachi, MS., « Relation effort - ouverture des fissures dans le cas
du béton armé de fibres d'acier », Annales de l'ITB, Ed. ESKA, France, N°2, avr. 1998, pp. 5-
14.
[100] Victor, Davidovici., «Formulaire du béton armé 1 », Edition du moniteur paris, 1996.
[101] Sheikh, S. A. and Uzumeri S., « Analytical Model for Concrete Confinement in Tied
Columns », journal of structural Division, ASCE,1982, vol .108,No.12 pp 2703-2722.
[102] Montoya E., Vecchio F.J and Sheikh S. A., « Compression Filed Modeling of Confined
Concrete », structural engineering mechanics, 2001, vol. 24,No.3 pp 231-248.
[103] Virlogeux anales l’I.T.B.T.P N°444, mai 1986
[104] Bouafia Y., Kachi MS., « calcul non linéaire jusqu’a rupture d’une section de poutre
soumise a une flexion plane », Journées scientifiques, CTC-Centre- Université M.M. de Tizi-
Ouzou, 25-26 juin 1997.
[105] Vecchio F.J., Collins M.P., « Predicting the response of reinforced concrete beams
subjected to shear using modified compression field theory», ACI Structural Journal, May-
June 1988, p. 258-268.
[106] Foure B., Rezende-Martins P.C. DE., « comportement en flexion jusqu'à rupture des
poutres à précontrainte extérieure », AFPCXI congrès de la FIP Hambourg, 1990.
145
[107] Foure B.,Bouafia Y., Soubret R., Thomas P., « shear test on keyed joints between
precast segments», Workshop behavior of external prestressing in structures.Saint Rémy-lés-
Chevreuse,France juin 1993.
[108] Garcia-Vargas JA., « la précontrainte extérieure appliquée à des poutres

hyperstatiques», Thèse de doctorat, Ecole centrale des Arts et Manufactures, Paris 1994.
[109] Grelat A., « comportement non linéaire et stabilité des ossatures en béton armé»,
Annales de l’I.T.B.T.P.N°336, Novembre 1978.
[110] Grelat A., « Analyse non linéaire des ossatures hyperstatiques en béton armé », Thèse
de doctorat Ingénieur, Université paris VI, 1978.
[111] Grelat A., Nait Rabah O., « simulation numérique du comportement des structures en
béton armé ou précontraint hors du domaine linéaire », Rapport du SES-CEBTP, Octobre
1988.
[112] R.Mahtout, A. Iddir,Y. Bouafia, Ms. Kachi., « la prise en compte du confinement dans
le calcul non linéaire des sections », J.S.UMMTO 20-21 November Université M. M. de Tizi-
Ouzou, Algérie.
[113] Gongchen D. Xuekang T., « contrainte ultime dans les câbles non adhérant de poutres
en béton à précontrainte partielle», Annales de l’I.T.B.T.P.N°412, Février 1988.
[114] Higer F.J., Mcgrath T.J., «Shear strength of pipe box section and other one-way flexural
members», ACI journal November-December 1982
[115] Hindi A., Kreger M. et Breen J., « Enhacing the strenth and ductility of post-tensioned
segmental box-grider bridges», Rapport de recherche n°1209-2F. Canter for transportation
research. Université du Texas à Austin, avril 1991.
[116] Hindi A., Kreger M. et Breen J., «Enhacing the strenth and ductility of post-tensioned
segmental box-grider bridges», Workshop behavior of external prestressing in structures.Saint
Rémy-lés-Chevreuse,France juin 1993.
146
[117] Hyo-Gyoung Kwak., Sun-Pil Kim., «Non linear analysis of RC beams based on
moment-curvature relation», Computers and structures.80,pp615-628,2002.
[118] Kachi M.S.,Bouafia Y.,Drizi S.,«Calcul non linéaire jusqu'à rupture d’une section en
béton armé, béton précontraint et section mixtes», Premier séminaire national en Génie civil.
Biskra Algérie. 8 et 9 décembre 1996.
[119] Adjrad A., Bouafia Y., Kachi MS., Dumontet H., « Modeling of Externally Prestressed
Beams until Fracture in Non Linear Elasticity » , Applied Mechanics and Materials, Vol. 749
(2015), pp 379-385.
[120] Kaverchin M., « Recherche sur les structures en béton », Annales de l’I.T.B.T.P.N°360,
Avril 1978.
[121] Kaverchin M., « Utilisation de la précontrainte extérieur » Annales de l’I.T.B.T.P. Juin

1983.
[122] Lim T.Y., Pramasivam P., Lee S., «Analytical model for tensile behavior of steel fiber
reinforced concrete J.ACI,July-August 1987.
[123] Macleod I.A.,Houmsi A., «shear strength of haunched beams without shear
reinforcement», ACI structural journal, Janury-Febrery 1994.
[124] Mansour M.Y., Dicleli M., Lee J.Y.,Zhang J., « Predecting the shear strength of
reinforced concrete beams using artificial neural networks», Engineering
Structures.26.(2004).PP781-799.
[125] Mahtout R., «Vulnérabilité des structures poteaux-poutres en béton armé vis-à-vis des
efforts sismiques», mémoire de magister, université mouloud mammeri, tizi ouzou, 2009.
[126] Richart FE, Brandtzaeg A, Brown RL., « A study of the failure of concrete under
combined compressive stresses », Engineering Experiment Station, Univ. ofIllinois, Urbana,
I11; 1928.
147
[127] Teng JG, Jiang T, Lam L, Luo YZ., « Refinement of a design-oriented stress–strain
model for FRP confined concrete », J Compos Constr ASCE 2009;13(4):268–78.
[128] Jiang T, Teng JG., « Analysis-oriented stress–strain models for FRP-confined

concrete », Eng. Struct 2007;29:2968–86.
[129] Teng JG, Huang YL, Lam L, Ye LP., « Theoretical model for fiber-reinforced polymer-
combined concrete », J Compos Constr 2007;11(2):201–10.
[130] Becque J, Patnaik AK, Rizkalla SH., « Analytical models for concrete confined with
FRP tubes ». J Compos Constr 2003;7(1):31–8.
[131] Harmon TG, Ramakrishnan S, Wang EH., « Confined concrete subjected to uniaxial
monotonic loading », J Eng Mech 1998;124(12):1303–8.
[132] Toutanji HA. , « Stress–strain characteristics of concrete columns externally confined

with advanced fiber composite sheets », ACI Mater J 1999;96(3):397–404.
[133] Saafi M, Toutanji HA, Li Z., « Behavior of concrete columns confined with fiber
reinforced polymer tubes », ACI Mater J 1999;96(4):500–9.
[134] Lillistone D, Jolly CK., « An innovative form of reinforcement for concrete columns
using advanced composites », Struct Eng 2000;78(23/24):20–8.
[135] Li YF, Lin CT, Sung YY., « A constitutive model for concrete confined with carbon
fiber reinforced plastics », Mech Mater 2003;35:603–19.
[136] Jin XN, Pan JL, Liu GY, Lai WH., « Research of stress–strain curve of concrete
confined by fiber reinforced plastics under axial compression ». J Build Struct2003;24(4):47–
53 [in Chinese].
[137] Lam L, Teng JG., « Design-oriented stress–strain model for FRP-confined concrete »,
Constr Build Mater 2003;17:471–89.
[138] Lam L, Teng JG., « Design-oriented stress–strain model for FRP-confined concrete in
rectangular columns », J Reinf Plast Compos 2003;22(13):1149–86.
148
[139] Harajli MH, Hantouche E, Soudki K. , « Stress–strain model for fiber-reinforced
polymer jacketed concrete columns », ACI Struct J 2006;103(5):672–82.
[140] Baris B., « An analytical model for stress–strain behavior of confined concrete »,
Engineering Structures 27 (2005) 1040–1051, February 2005.
[141] Orito Y, Sato T, Tanaka N, Watanabe Y., « Study on the unbonded steel tube composite
system. In: Composite construction in steel and concrete proceedings », ASCE engineering
foundation. 1987,p. 786–804.
[142] McAteer P, Bonacci JF, Lachemi M., « Composite response of highstrength concrete
confined by circular steel tube. ACI, Structural », Journal 2004;101(4):466–74.
[143] Saad, M., Kachi, MS.,Bouafia, Y., Muller, P., Fouré, B., « Influence du pourcentage
d'acier sur le comportement du béton tendu fissuré », European Journal of Environmental and
Civil EngineeringVolume 14, Issue 3, March 2010, pages 303-327.
[144] DiPrisco M .,Gombrova G.,« Comprehensive model for stady of shear in thin webbed
RC and PC beams », Journal of structural engineering. December 1995.
[145] Falkner H.,Teutsch M.,Huang Z., « Prestressed segmental box girders with unbounded
tendons under combined torsion bending and shear », workshop behavior of external
prestressing in structures . Saint Remy-les-Cheuvreuse,jouin 1993.
[146] Gregor T.,Collins M.P., « Test of large partially prestressed concrete girders », ACI
structural journal January-February 1995.
[147] Heger F. J.,Mcgrath T.J., « Shear strength of pipe box section and other one-way
flexural members », ACI Journal November-December 1982.
[148] HUANG J., EIBL J., « Design of segmental bridges under combined bending, shear
and torsion fe-study », Workshop Behavior of external prestressing in structures. Saint
Remy- Les-Chevreuse,juin 1993.
149
[149] Keun Hyeok Yang,Heon-Soo Chung,Eun-Taik Lee,Hee-Shang Eun, « shear
characteristics of high-strength concrete deep beams without shear
reinforcement »,Engineering Structures.25.(2003).PP1343-1352.
[150] Kishi N., Nikami H., Hatsuoka K.G., Ando T., « Impact Be havior of Shear-Failure-
type RC beams without shear bar », International Journal of impact Enginnering.27
(2002).PP955-968. 92.
[151] Mau S.T., Thomas T.C.Hsu, « Shear design and analysis of low-rise structural
walls », ACI Structural Journal March-April 1986.
[152] Ramirez G.,Macgregor R., Kreger M.E.,Robberts-Wollmann R.,Breen J., « Shear

strength of segmental structures», Workshop Behavior of external prestressing in structures.
Saint Remy-les-Cheuvreuse,juin 1993.
[153] Salandra M.A.,Shuaib H.A., « Shear capacity of reinforced lightweight high strength
concrete beams », ACI structural Journal November-December 1989.
[154] Stanciulescu-Panea I., « Modélisation de l’influence de l’effort tranchant dans le

comportement non linéaire des poutres en béton armé », mémoire de DEA de l’Ecole
Nationale des Ponts et Chaussées, Juillet 1996.
[155] Stroband J., « Shear capacity of high strength concrete beams with shear
reinforcement », Progress in concrete research vol num 5 1997.
[156] Thomas T.C.Hsu, « softened truss m odel theory for shear and torsion », Rechearch
Rapport UHCE-86-9 University Of Houston August 1986.
[157] Thomas T.C.Hsu, « Non-linear analysis of concrete membranes elements », ACI

Structural Journal September-October 1991.
[158] Thomas T.C.Hsu, Mo Y.L., « Softening concrete in low-rise shear walls », ACI
Structural Journal November-December 1985.
150
[159] Eurocod 2, Env 1992-1-1,Nf P 18 711, « Calcul des structures en béton,Partie 1-
1:Règles générale s et règles pour les bâtiments », Décembre 1992.
[160] Panagiotis Kotronis, « Cisaillement Dynamique De Murs En Beton Arme Modeles

Simplifies 2D et 3D », thèse de doctorat-Ecole Normale Supérieure de Cachan ; décembre
2000.
[161] Davidovici V. et al. , « Génie Parasismique », Presses de l'Ecole Nationale des Ponts et
Chaussées, 1985.
[162] Martin-Pérez B. & Pantazopoulou S. J ., « Mechanics of concrete participation in
cyclic shear resistance of RC », Journal of Structural Engineering, Vol. 124, No 6, June, pp.
633-641, 1998.
[163] Lampert P. & Thurlimann B., « Ultimate strength and design of reinforced concrete
beams in torsion and bending », IABSE, No. 31-I, 107,131, 1971.
[164] Moy.L. & Rothert. H., « Effect of softening m odels on behavior of reinforce concrete
framed shearwalls », ACI Structural journal, V.94, No.6, Novem ber- December, 1997.
[165] Davenne L., «Modélisation de l'influence des armatures transversales sur le

comportement non linéaire d'éléments de béton armé», Thèse Université Pierre et Marie
Curie, Paris 6, 1990.
[166] Baby F, Marchand P, Atrach M, Toutlemonde F., « Analysis of flexure-shear behavior

of UHPFRG beams based on stress field approach»,Engineering structures 56 (2013) pp 194-
206.
[167] Iddir, A., Bouafia,Y., Kachi. Ms., « Simulation of the influence of the confinement,
until fracture, on the shear stiffness of the reinforced or prestressed concrete», 20th European
Conference on Fracture (ECF20), Procedia Materials Science 3 ( 2014 ) 679 – 685.
[168] Bouafia Y.,« Résistance à l'effort tranchant des poutres à précontrainte extérieure -
Etude expérimentale et calcul à la rupture », Thèse de Docteur (Nouvelle thèse), Ecole
Centrale Paris, 288 p, Nov. 1991.
151
[169] Bouafia Y.,« Simulation numérique du comportement moyen, jusqu'à rupture, d'une
zone de poutre. Application au béton armé ou précontraint et béton de fibres», Mémoire
DEA, Université de Paris 6, 77 p, juin 1987.
[170] Adjrad A., Bouafia Y., Kachi MS, Dumontet H., « Modeling of internally or externally
prestressed concrete beams until fracture in non linear elasticity», IACSIT, International
Journal of Engineering and Technology (IJET), july, 2014.
[171] Saad M., Bouafia Y., Kachi M. S. « Contribution à l’évaluation d’ouverture des fissures
dans les éléments en béton armé », Annales du bâtiment et travaux publics, Ed. ESKA, n° 4,
Juin 2014.
[172] Adjrad A., Kachi, MS., Bouafia Y., Iguetoulène F., « Nonlinear modeling structures on
3D », The International Conference on Structural Analysis of Advanced Materials, ICSAAM
2011, 7-11 september, Sinaia, Romania, pp. 1 – 9.
[173] Adjrad A., Bouafia Y., Kachi MS., Dumontet H., « Modeling of Externally Prestressed
Beams until Fracture in Non Linear Elasticity », International Conference on Intelligent
Materials and Manufacturing Engineering, IMME 2015 - M060, Phuket Island, Thailand, 16 –
17 January 2015.
[174] Bouafia, Y.,Saad, M., Kachi, MS., « Relations contraintes – déformations moyennes et
ouvertures des fissures dans le cas d’un tirant en béton armé ». Société Marocaine des
Sciences Mécaniques - Réseau Universitaire de mécanique, El jadida - Maroc, 8 è Congrès de
Mécanique, 17 - 20 avril 2007.
[175] Bouafia, Y.,Saad, M., Kachi, MS., « Influence du pourcentage d’acier sur le
comportement du béton fissuré en traction ». Société Marocaine des Sciences Mécaniques -
Réseau Universitaire de mécanique, Maroc (Casablanca), 7 è Congrès de Mécanique, 19 - 22
avril 2005.
152
Abstract
In the civil engineering, the calculation of circular sections is not easy because there are many unknowns in the
equations of equilibrium. The required experience, advantages and disadvantages and problems to solve have
been widely discussed by many authors. The calculation of such sections is difficult to predict exactly there
behavior. The effect of the tensioned concrete, even after cracking, has a great importance. The effect of the
confinement of concrete generally gives this one greater resistance and ductility according to the level of
confinement. In this study we propose a model which is based on the equilibrium, compatibility equations and
which are formulated by average stresses and average strains. According to the assumptions. This model is able
to analyze sections having unusual forms or reinforcing details, loaded in combined bending, axial load and
shear. In the case of the confined concrete, a stress-strain relation of this one is proposed. In this case, for the
purpose of investigating confinement effect, an effectively confined ratio was introduced and the effects
according to concrete compressive strength. Predictions of the model are compared with several calculated
unconfined reinforced or prestressed beams and are shown to estimate the effect of the confinement on the shear
stiffness of the sections in the case of the elastic linear domain, after concrete cracking and after reinforcement
yielding. The results of the comparison give a good satisfaction.
Keywords: nonlinear elasticity, confinement, circular section, shear force, shear stiffness, ties discrete,
equilibrium, compatibility, cracks.
Résumé
Dans le domaine de génie civil le calcul des sections circulaires n’est pas aisé compte tenu de la disposition des
armatures qui induit plusieurs inconnues dans les équations d’équilibre. L’expérience acquise, les avantages et
les inconvénients ainsi que les problèmes à résoudre ont été largement discutées par de nombreux auteurs. Le
calcul de telles sections présente des incertitudes sur certains aspects de leur comportement. L’effet de
l’influence du béton tendue même après fissuration est d’une grande importance. L’effet confinement du béton
en général lui confère une plus grande résistance et une plus grande ductilité en fonction du niveau du
confinement. Dans le cadre de cette étude nous proposons un modèle de calcul basé sur les équations
d’équilibres et de compatibilités qui sont formulées en termes de contrainte et déformation moyenne selon des
hypothèses. Ce modèle est capable d’analyser des sections quelconques avec tous les détails de ferraillage sous
l’ensemble des sollicitations (effort normal, moment fléchissant et effort tranchant). Une loi de comportement
contrainte- déformation du béton confiné a été proposée. Le confinement du béton sera introduit par le biais d’un
coefficient de confinement effectif, et d’une contrainte de confinement latérale. La méthode de calcul est injectée
dans un logiciel de calcul non linéaire des éléments de structures afin d’étudier la contribution et l’influence du
confinement du béton sur le comportement global des éléments de structures de section circulaire et sur la
rigidité au cisaillement dans le cas des poutres en béton armé et précontrainte. Enfin, le modèle de prise en
compte de l’effet de l’effort tranchant dans le domaine non linéaire a été également proposé, ce dernier permet
de prédire le comportement des poutres en cisaillement combiné à l’existence des contraintes normales, ce qui
permet alors d’estimer la rigidité réelle de cisaillement des poutres, pour un niveau de charge donné.
Mots-clés : Elasticité non linéaire, Confinement, section circulaire, effort tranchant, rigidité au cisaillement,
armatures discrètes, équations d’équilibres, équations de compatibilités, rupture.

These Iddir Abdelkader PDF

Transféré par

Droits d'auteur :

Formats disponibles

These Iddir Abdelkader PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

These Iddir Abdelkader PDF

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Université Mouloud Mammeri de TIZI-OUZOU

Faculté du Génie et de la Construction

Modélisation des éléments de structure de section

Devant le jury composé de :

Abdelkader TAHAKOURT Professeur, U. Béjaia Président

Je souhaite exprimer ma reconnaissance à mon co-directeur de thèse le Pr BOUAFIA Youcef

Je remercie le Pr TAHAKOURT Abdelkader qui me fait l’honneur de présider le jury, je

Je remercie mes proches pour leurs soutiens et leurs encouragements.

A la mémoire de mon grand-père et de ma mère

Tableau III.1:Caractéristiques des données expérimentales des sections circulaires en béton.

Aaj : section d’une barre d’ordre j

De plus, l’influence de la prise en compte du confinement dans la modélisation et le calcul

La présente étude est structurée comme suit :

- le chapitre 2 est consacré au développement du modèle de calcul non linéaire du

I-1.Première partie : Le confinement

I.1.2. Principe du confinement

I.1.3. Types de confinement

I.1.3.1. Le confinement interne

I.1.3.2. Le confinement externe

Figure I.3: béton armé confiné par FRP et tube en acier.

I.1.4. Poteaux mixtes

poteaux remplis de béton. Vu la souplesse de sa méthode Y.C.Wang [3] à conclu qu’elle

I.1.4.2. Effets de chargement, et de sollicitations extérieurs [2]

I.1.4.3. Tube en acier [2]

I.1.5. Quelques modèles de comportement du béton confiné en élasticité non linéaire

I.1.5.1. Modèle de Richart

I.1.5.2. Modèle d’Hognestad

Figure I.4 : Courbe contrainte-déformation du modèle Hognestad [18].

I.1.5.3. Modèle Bousalem et al .[19]

Figure I.5: Courbe contrainte-déformation du modèle Bousalem et al [19].

I.1.5.4. Modèle d’Ahmad et Shah

L’équation équivalente est:

'  f '  4 .255 f

I.1.5.5. Modèle de Toutanji

Où : fcom est la résistance à la traction du composite confiné.

I.1.5.6. Modèle Mander et al.

 7,94 f ' l 2 f ' l 

Figure I.6: Courbe contrainte-déformation du modèle Mander et al.

I.1.5.7. Modèle Saatcioglu et al. [25]

 cc  f cc  Z  c   cc 0   0.2 f cc (Partie descendante)

I.1.5.8. Modèle Hoshikuma et al. [26]

 cc  f cc  E des  c   cc   0.5 f cc (Partie descendante)

I.1.5.9. Modèle de Fardis et Khalili

Où Ecom est le module d'élasticité du confinement dans la direction transversale.

I.1.5.10. Modèle de Karbhari

I.1.5.11. Modèle de Mirmiran

  0.872 f  0.371 f  6.258 MPa 

I.1.5.12. Modèle de Xiao et Wu

I.1.5.13. Modèle de Harmon

   elastic   crack   void

I.1.5.14. Modèle Fujikake

  1   13 31   23 32   12 21   21 13 31   12 23 31

I.1.5.15. Modèle de Beque

Ce modèle dépend essentiellement du module Ks et du module de cisaillement Gs.

Les différentes confrontations effectuées par ces auteurs montrent clairement :

- Le confinement est introduit par un coefficient de confinement effectif, le pourcentage

- Le domaine d’application de ces modèles se limite à des pourcentages de confinement allant