Mécanique des structures

et résistance des matériaux

Christophe Fond
Université de Strasbourg

Laboratoire ICube
Préambule ............................................................................................................................... 3
1. Introduction ........................................................................................................................ 3
1.1. Nature de différents types de problème en mécanique des milieux continus ............. 3
1.2. Quelques considérations pour le génie civil ................................................................ 4
1.3. A propos de la théorie des poutres .............................................................................. 5
1.4. Principales notations ................................................................................................... 6
1.5. Déformations ............................................................................................................... 6
1.5.1. Admissibilité de la cinématique ........................................................................... 6
1.5.2. Allongement - raccourcissement .......................................................................... 7
1.5.3. Distorsion ............................................................................................................. 7
1.5.4. Tenseur des déformations..................................................................................... 7
1.5.5. Equations de compatibilité ................................................................................... 8
1.6. Contraintes .................................................................................................................. 8
1.6.1. Vecteur contrainte ................................................................................................ 8
1.6.2. Tenseur des contraintes ........................................................................................ 8
1.6.3. Réciprocité des contraintes tangentielles ............................................................. 9
1.6.4. Admissibilité de la statique .................................................................................. 9
1.6.5. Représentation de l'état de contrainte en un point par les cercles de Mohr ......... 9
1.7. Principe fondamental en MMC ................................................................................. 10
1.7.1. Equilibre ............................................................................................................. 10
1.7.2. Puissances virtuelles et travaux virtuels ............................................................. 11
1.8. Solutions en MMC .................................................................................................... 12
2. Comportement des matériaux – résistance des matériaux ............................................... 12
2.1. Classification du comportement mécanique des matériaux ...................................... 12
2.1.1. Préambule ........................................................................................................... 12
2.1.2. Définitions .......................................................................................................... 13
2.1.3. Modèles rhéologiques ........................................................................................ 16
2.1.4. Plasticité dans les métaux : micromécanismes de base et "comportement
émergent" ................................................................................................................................. 16
2.2. Notion de thermomécanique ..................................................................................... 17
2.2.1. Coefficient de dilatation thermique .................................................................... 17
2.3. Elasticité linéaire – loi de Hooke .............................................................................. 17
2.3.1. Elasticité isotrope ............................................................................................... 17
2.3.2. Thermo-élasticité isotrope .................................................................................. 18
2.3.3. Elasticité anisotrope ........................................................................................... 18
2.3.4. Limite élastique et critères de plasticité ............................................................. 18
2.4. Rupture ...................................................................................................................... 19
2.4.1. Rupture sous chargement monotone .................................................................. 19
2.4.2. Rupture sous chargement cyclique, fatigue........................................................ 19

1
3. Exemples de réalisations .................................................................................................. 21
4. Modélisation des conditions aux limites .......................................................................... 25
4.1. Déplacements imposés .............................................................................................. 25
4.2. Pression, contraintes imposées, surface libre ............................................................ 25
4.3. Liaisons extérieures ................................................................................................... 26
4.4. Appuis élastiques....................................................................................................... 27
4.5. Liaisons intérieures ................................................................................................... 27
4.5.1. Cas de la construction métallique....................................................................... 27
4.5.2. Cas du béton armé et/ou précontraint ................................................................. 28
4.5.3. Cas de la construction en bois – lamellé collé ................................................... 29
4.5.4. Cas de la construction mixte .............................................................................. 30
4.6. Remarques ................................................................................................................. 30
5. Spécificités de la théorie des poutres ............................................................................... 30
5.1. Poutre – barre ............................................................................................................ 30
5.1.1. Définition ........................................................................................................... 30
5.1.2. Equilibre des poutres et barres ........................................................................... 30
5.2. Hypothèse de Navier-Bernoulli ................................................................................. 31
5.2.1. Cinématique ....................................................................................................... 31
5.2.2. Déformation d'une section .................................................................................. 31
5.2.3. Admissibilité ...................................................................................................... 32
5.3. Principe de St Venant, torseur généralisé .................................................................. 33
5.3.1. Equilibre d'une section ....................................................................................... 33
5.3.2. Relation entre l'effort tranchant et le moment fléchissant .................................. 34
5.3.3. Répartition des contraintes dans une section...................................................... 34
5.3.3.1. Contraintes normales ................................................................................... 34
5.3.3.2. Moment quadratique et théorème d'Huygens .............................................. 35
5.3.3.3. Effort normal et Moment fléchissant .......................................................... 35
5.3.3.4. Contraintes tangentielles ............................................................................. 36
5.3.3.5. Diagrammes de sollicitations ...................................................................... 36
5.3.4. Déformation d'une poutre ................................................................................... 36
5.3.5. Remarque ........................................................................................................... 36
5.4. Energie de déformation ............................................................................................. 36
5.4.1. Energie et élasticité ............................................................................................ 36
5.4.2. Energie de déformation liée à l’effort normal .................................................... 37
5.4.3. Energie de déformation liée à l’effort tranchant ................................................ 38
5.4.4. Energie de déformation liée au moment fléchissant .......................................... 38
5.4.5. Energie de déformation élastique dans une poutre............................................. 38
6. Résolution des problèmes de poutres droites -systèmes isostatiques plans ..................... 39
6.1. Equation de la déformée ............................................................................................ 39
6.2. Déplacement ou rotation en un point ........................................................................ 40
6.3. Théorème de Vereschiaguine, intégrales de Mohr .................................................... 41
7. Résolution des problèmes de poutres droites - systèmes hyperstatiques plans ................ 41
7.1. Degré d'hyperstaticité ................................................................................................ 41
7.2. Résolution par la méthode des forces ........................................................................ 42
7.3. Formule des trois moments, démonstration à partir de la méthode des forces ......... 42
8. Méthode des déplacements............................................................................................... 44
8.1. Introduction ............................................................................................................... 44
8.2. Solution de base ........................................................................................................ 45
8.2.1. Composante liée à l'effort normal ...................................................................... 45
8.2.2. Exemple pour un système triangulé ................................................................... 46
8.2.3. Composantes liées à la flexion ........................................................................... 48
8.2.4. Matrice de rigidité d'une poutre ......................................................................... 50

2
 8.2.5. Exemple pour un système fléchi ........................................................................ 50
8.2.6. Formule des trois moments à partir de la méthode de déplacements ................. 51
8.2.7. Exemple de mise en œuvre de la méthode dite des 3 moments ......................... 52
8.3. Formulaire ................................................................................................................. 52
9. Flambement et instabilités de forme ................................................................................ 53
9.1. Hypothèses de calcul ................................................................................................. 53
9.2. Cas de la compression simple – flambement d'Euler ................................................ 53
10. Exemple de résolution d’un problème de portique ........................................................ 54
10.1. Résolution par la méthode des forces ...................................................................... 54
10.2. Remarque concernant les déplacements des têtes de poteau dans des portiques .... 57
Annexe 1 : Un peu de géométrie pour un tétraèdre.............................................................. 58
Annexe 2 : Section réduite liée au cisaillement ................................................................... 59
A3. Démonstration inverse du TFS à partir du PTV… ........................................................ 60
Annexe 4 : En quoi la théorie des poutres n'est pas une solution exacte ?........................... 60
A4. 1 Gauchissement des sections ................................................................................... 60
A4. 2 Flexion circulaire .................................................................................................... 61
A4. 3 Loin des points d'application des charges… .......................................................... 61
Annexe 5 : Théorème des travaux virtuels, que serait la démonstration pour une poutre
droite ?.............................................................................................................................................. 62
A5. 1 Préambule ............................................................................................................... 62
A5. 2 Champs cinématiquement admissible et statiquement admissible ......................... 62
A5. 3 Travail des efforts intérieurs................................................................................... 63
A5. 4 Termes liés à l'effort normal ................................................................................... 63
A5. 5 Termes liés à l'effort tranchant ............................................................................... 64
A5. 6 Termes liés au moment fléchissant ........................................................................ 64
A5. 7 Total du travail des efforts intérieurs ..................................................................... 65
Annexe 6 : Relation de conversion entre rayon de courbure et contraintes internes
résiduelles......................................................................................................................................... 66
A6. 1 Introduction ............................................................................................................ 66
A6. 2 Enoncé du problème ............................................................................................... 66
A6. 3 Equations d’équilibre de la poutre en absence de force extérieure ........................ 67
A6. 4 Cas "classique" du module d'Young uniforme ....................................................... 69
A6. 5 Moyenne du module d'Young mesuré par essai de flexion .................................... 70
A6. 6 Contrainte résiduelle éventuellement additionnelle liée à l'enlèvement de lamelle
...................................................................................................................................................... 70
A6. 7 Références bibliographiques de cette annexe ......................................................... 70

Préambule
Ce cours résume certaines notions de mécanique des milieux (MMC) continus déjà
abordées de sorte que la continuité apparaisse entre la MMC et la théorie des poutres. Les notions
déjà vues seront revisitées très rapidement.
Il est plus commode en MMC d'utiliser des indices (x1, x2, x3) pour repérer les axes que des
lettres différentes (x, y, z). Toutefois, pour certaines parties de ce cours, les indices servent à
repérer des objets et dans ce cas nous aurons recours à cette dernière notation.

1. Introduction
1.1. Nature de différents types de problème en mécanique des milieux continus
Pour qu'une structure demeure pérenne, il faut éviter de l'endommager. On cherchera donc
à rester dans le domaine d'utilisation des matériaux où les phénomènes sont totalement réversibles.

3
En pratique, il ne faudra pas dépasser la limite élastique des matériaux. On parlera de
comportement linéaire élastique. Dans le cas contraire, on parlera de non-linéarité "matériau".
Pour qu’une structure soit peu déformable, il faut éviter les grandes déformations du
matériau ainsi que les grands déplacements de la structure. On parlera alors des hypothèses des
petites perturbations (H. P. P.). Dans le cas contraire, on parlera de non-linéarité "matériau" en
grandes déformations1 et/ou de non-linéarité "géométrique" en grands déplacements.
Les variations de température induisent aussi des déformations et variations
dimensionnelles des structures. Les dilatations thermiques différentielles peuvent générer des
contraintes considérables dans des structures et des matériaux. Elles peuvent être causées par la
non-uniformité du champ de température et/ou la non-uniformité des propriétés de dilatation
thermiques des matériaux. On parle alors de problème thermo-mécanique. Dans certains cas, il
faut de plus considérer les variations du comportement mécanique avec la température, ce qui
conduit aussi à une forme de non-linéarité de comportement.
Pour la plupart des réalisations, la structure est réutilisée et les grandes déformations et
grands déplacements sont indésirables. Le matériau est donc souvent sollicité dans le domaine
d'élasticité linéaire. Par ailleurs, les effets de la température peuvent souvent être négligés. Dans
les cas d’H. P. P. et d’élasticité linéaire, quelle que soit l'histoire de chargement de la structure
pour passer d'un état de sollicitation à un autre, on aboutit toujours au même et unique état. Par
ailleurs, il sera dans ces cas loisible de superposer les sollicitations mécaniques et leurs effets.

1.2. Quelques considérations pour le génie civil

Les constructions doivent assurer la sécurité des personnes et des biens. Les règlements de
construction considèrent plusieurs situations dont les principales sont les états limites de service
(ELS) et les états limites ultimes (ELU). Aux ELS il s'agit de ne jamais générer de déformation
permanente ce qui conduit à rester dans le domaine élastique des matériaux. Aux ELU on admet
que la structure peut avoir subi des dommages causés par des événements rares mais elle doit avoir
assuré la sécurité. Dans ce dernier cas, la sortie du matériau du domaine d'élasticité permet
d'atteindre les limites du matériau et de changer la distribution des efforts dans la structure.

Ruine de la structure par rupture du matériau en Ruine de la structure par rupture du matériau en
traction sous l'effet de la flexion. traction sous l'effet du cisaillement.

La ruine d'une structure peut-être liée aux limites de résistance des matériaux de
construction. Elle peut aussi être liée à des instabilités de forme telles que flambement, voilement,
déversement… Enfin, un manque de rigidité d'un ou plusieurs éléments de la structure peut
conduire à des nuisances ou des dégâts.

1
Sauf cas très particulier, un matériau subissant une grande déformation présente une réponse mécanique non
linéaire.

4
 Ruine de la structure par flambement, la Souplesse excessive d'une partie de la structure
déformation amplifie la sollicitation. entraînant des déformations incompatibles.

1.3. A propos de la théorie des poutres

Quels ingrédients faut-il pour construire la théorie des poutres classique ?
- définition d’une poutre,
entre autres une dimension "grande" devant les deux autres,
- hypothèse de Navier : les sections droites restent droites, les fibres parallèle à la fibre moyenne
demeurent parallèles à la fibre moyenne,
excellente approximation de la réalité, valable même en grands déplacements et au-
delà de l’élasticité linéaire,
- hypothèse de Saint Venant : on est "loin" d’une charge ponctuelle,
au voisinage des charges ponctuelles les sections se gauchissent donc la théorie donne
de moins bonne prédictions,
- les fibres sont en situation de traction ou compression uniaxiale, on se donne donc comme loi de
comportement contrainte = module d’Young * déformation,
le coefficient de Poisson n’intervient donc pas dans cette théorie,
- on doit assurer l’équilibre,
pour simplifier les calculs, on l’écrira généralement au niveau de la fibre moyenne
afin d’éviter un couplage entre effort normal et moment fléchissant,

Il suffit ensuite de faire un peu de géométrie, postuler l’équilibre et exprimer la loi de

comportement pour tout reconstruire. La lecture de ce cours – et bien d’autres - montre comment
s’y prendre.

Les sections droites restent droites.

5
 Point d’application d’une charge ponctuelle, les sections droites se gauchissent.

Au-delà de l’étude de la flexion, il conviendrait d’étudier la torsion. En général on évite ce

genre de sollicitation sur les éléments de structure de génie civil. En génie mécanique, les efforts et
puissances sont souvent transmis par des arbres soumis à de la torsion. C’est donc dans ce
domaine de formation qu’on trouvera naturellement le plus de cours sur la torsion. Nonobstant, il
faudra aux futurs ingénieurs du génie civil découvrir, au-delà de ce cours, les notions de base
concernant la torsion. En effet, il est indispensable de poursuivre l’apprentissage par l’étude de la
théorie des plaques. Celle-ci fait intervenir les notions de flexion et de torsion...

1.4. Principales notations

Les vecteurs et tenseurs sont en gras, les tenseurs d'ordre deux sont soulignés deux fois.
Les opérateurs grad, div, etc, sont en italique. x1, x2 et x3 représentent les trois directions de
l'espace dans un repère orthonormé direct. Id représente le tenseur identité tel que Idii = 1 et
Idij = 0 si i  j.
notation unité désignation notation unité désignation
surface de 
3
 m domaine spatial  m²
ij accélération = 2u/t2
3
m/s²  Kg/m densité
fi N/m3 force volumique
ij ij
2
N/m contrainte (vraie) de - ou déformation (petite) de
ou Pa Cauchy m/m Green
2
T N/m vecteur contrainte F N force
P N/m charge répartie selon x1 q N/m charge répartie selon x2
N, Nij N effort normal A m² aire d'une section
V, Vij N effort tranchant Sr m² section réduite
M, Mij Nm moment fléchissant IGZ m4 Moment quadratique
T °K température d'une section (ou inertie)
W J énergie ou travail y(x1) m flèche en x1
dans la notation (x1, x2, x3) dans la notation (x, y, z)
ui m déplacement d'axe i ui m déplacement de i selon x
vi m/s vitesse d'axe i = ui/t vi m déplacement de i selon y

Un fond de page grisé indique qu'il s'agit de notions déjà vues, en général en cours de
mécanique des milieux continus et/ou non spécifiques au génie civil. Ces notions ne sont
généralement pas détaillées en cours.

1.5. Déformations
1.5.1. Admissibilité de la cinématique
Un champ de déplacements est dit cinématiquement admissible s'il respecte les conditions
aux limites, est dérivable et ne génère pas d'interpénétration de matière ni de vide.

6
 1.5.2. Allongement - raccourcissement
Il existe deux façons de considérer les mouvements de matière. La description eulérienne,
principalement utilisée en mécanique des fluides, considère les points de l'espace et la vitesse
matérielle en ces points. La description lagrangienne, principalement utilisée en mécanique des
solides, considère la vitesse ou le déplacement des points matériels dans l'espace. Sous les H. P. P.,
on assimile la position actuelle dans l'espace des points matériels à leur position initiale et les
descriptions se rejoignent. En conséquence les rotations d'axes x1, x2 et x3 attachés à la matière
sont considérées comme négligeables. L'expression d'une déformation longitudinale est donc
11 = u1/x1. Les translations et rotations solides n'induisent pas de déformation. Il existe
différentes façons de représenter les déformations (L/L0, L/L0 + ½ (L/L0)², ln(L/L0), …)
mais sous les H. P. P. toutes se rejoignent.

espace dessin, allongement d'une fibre

 u1   u1 
u0 =  u2  uP =  u2 
 u3 0  u3 P

1.5.3. Distorsion
Sous les H. P. P., l'expression d'une distorsion est 12 = ½ (u1/x2 + u2/x1) = ½ 12.

espace dessin, élimination de la rotation solide

1.5.4. Tenseur des déformations

Sous les H. P. P., l'expression du tenseur des déformations est :
ij = ½ (ui/xj + uj/xi)   = ½ (grad u + gradT u)
Le tenseur des déformations est symétrique, i. e. ij = ji. Il existe un repère principal dans lequel
les termes non diagonaux s'annulent. Dans ce repère, on note conventionnellement les
déformations principales 1  2  3. Pour trouver le repère principal il faut résoudre l'équation :
 n =  n  (   Id) n = 
Pour un tenseur symétrique, dét(   Id) = 0 admet 3 valeurs propres réelles qui sont les
déformations principales 1, 2 et 3. En effet :
 11 12 13   n1   n1   11   12 13   n1   0 
    . n  =  n          . n  = 0
 12 22 23   2   2  12 22 23
  2  
 13 23 33   n3   n3   13 23 33     n3   0 
pour que ce système admettent une autre solution que (0, 0, 0) il faut que son déterminant soit nul.
Dans le repère principal, on peut visualiser l'ellipsoïde des dilatations :

7
 1 = 2 = 3 1 > 2 > 3
 1     1  1  1     0  0  1     0  0
     . 0 =  0      . 1 =  1      . 0 =  0
 2
   1   2
   2   2
   3 
   3   0  0    3   0  0    3   1  1
Dans la direction principale 1, L / L0 = (L0+L) / L0 = 1 + 1
Dans la direction principale 2, e / e0 = (e0+e) / e0 = 1 + 2
Dans la direction principale 3, h / h0 = (h0+h) / h0 = 1 + 3

1.5.5. Equations de compatibilité

Les déformations étant liées par dérivations aux déplacements, elles ne peuvent pas évoluer
indépendamment. Les déformations doivent donc respecter 6 équations, obtenues par permutation
circulaire à partir de :
²ij/xkxl  ²ik/xjxl  ²jl/xixk  ²lk/xixj

1.6. Contraintes
1.6.1. Vecteur contrainte
Sous les H. P. P., on n'actualise pas les surfaces sur lesquelles s'exercent les forces car leurs
variations sont considérées négligeables. Si une force dF est appliquée sur un élément de surface
ds, le vecteur contrainte au point M est défini par :
V = ds0
lim F / ds = dslim F / ds0
00

où ds0 désigne l'élément de surface dans l'état initial et ds dans l'état actuel

espace dessin, vecteur contrainte en O

1.6.2. Tenseur des contraintes

Le tenseur des contraintes permet de connaître le vecteur contrainte sur une facette orientée
par n. Il est défini par :

8
 .n=T

espace dessin, équilibre du tétraèdre, tenseur des contraintes

Bilan des forces selon x1 :
2
11 ½bc  21 ½ab  31 ½ac + T1 ½ a12 b22 + a12 c32 + c32 b22 = 0 d'où
2
T1 = (11 bc  21 ab  31 ac) / a12 b22 + a12 c32 + c32 b22 = 11 n1  21 n2  31 n3
A l'instar des déformations, il existe un repère principal dans lequel les contraintes de
cisaillement s'annulent. Le repère principal est le même pour les contraintes que pour les
déformations. Par convention 1  2  3. La contrainte normale vaut n = n . T = n .  . n. La
contrainte tangentielle est la complémentaire à T dans le plan de la facette considérée, i. e.
n t = T  (T . n) n.

1.6.3. Réciprocité des contraintes tangentielles

Sauf cas très particulier de couple magnétique, l'équilibre des moments pour élément de
volume impose que ij = ji d'où la réciprocité des contraintes tangentielles.

espace calcul, pas de couple à l'échelle microscopique

1.6.4. Admissibilité de la statique

Un champ de contraintes est dit statiquement admissible s'il respecte les conditions aux
limites en effort et l'équilibre en tout point.

1.6.5. Représentation de l'état de contrainte en un point par les cercles de Mohr

Dans un plan (xi, xj), la normale est définie par (cos , sin ). Dans ce plan, il peut être
commode d'exprimer le vecteur contrainte dans le repère principal à partir de la contrainte

9
moyenne m = ½ (i + j) et  = ½ (i  j). La composante du vecteur contrainte dans ce plan
vaut (cos  i, sin  j). La composante normale n vaut (cos²  i + sin²  j) et la composante
tangentielle n vaut (sin  cos  i + cos  sin  j).
n = cos²  (m + ) + sin²  (m  ) = m +  (cos 2
n = m (cos  sin   cos  sin  ) +  (cos  sin   cos  sin  ) =  (sin 2
Il apparaît l'équation d'un cercle de centre (m, 0) et de rayon  dans le plan (n, n). La contrainte
tangentielle maximale correspond donc à  = 45° et vaut ½ (i  j). On peut procéder à des
permutations circulaires entre i, j et k pour couvrir les 3 plan définis par les trois directions
principales. On constate que le cisaillement maximal vaut max = ½ (1  3).

espace dessin, facettes perpendiculaires diamétralement opposées

1.7. Principe fondamental en MMC

1.7.1. Equilibre
L'équilibre peut être dit "dynamique" ou "statique" selon que, respectivement, l'on
considère les effets inertiels ou non. Les équations écrites en "dynamique" sont toujours valables.
On en tire celles de la "statique" en négligeant les termes inertiels.
Principe fondamental de l'équilibre
div  + f    = 0
Cette équation traduit l’équilibre des forces. L’équilibre des moments est lié à la symétrie du
tenseur des contraintes (voir la section réciprocité des contraintes tangentielles). En statique on
néglige le terme d'accélération pour obtenir : div  + f = 0
En absence de force volumique (pesanteur, …) on obtient en statique : div  = 0

espace dessin, id. avec les forces de volume

10
On parlera de puissances en dynamique et de travaux en statique.

1.7.2. Puissances virtuelles et travaux virtuels

Postulant le principe l'équilibre en tout point matériel qui est exprimé par :
div  + f    = 0
On peut multiplier 0 par un vecteur vitesse quelconque v* et intégrer sur un volume, on obtient
toujours 0 :
 {(div  + f   ) . v*} d = 0
En développant, on obtient :
 {(11/x1 + 12/x2 + 13/x3 + f1   2u1/t2) v1*
+ (21/x1 + 22/x2 + 23/x3 + f2   2u2/t2) v2*
+ (31/x1 + 32/x2 + 33/x3 + f3   2u3/t2) v3*} dx dx dx = 0

On peut réécrire :
 {11/x1 v1* + 21/x1 v2* + 31/x1 v3* + (f1   2u1/t2) v1*
+ (12/x2 v1* + 22/x2 v2* + 32/x2 v3* + (f2   2u2/t2) v2*
+ (13/x3 v1* + 23/x3 v2* + 33/x3 v3* + (f3   2u3/t2) v3*} dx dx dx = 0

Si v* est dérivable, en intégrant par partie on obtient :

 {(f1   2u1/t2) v1* + (f2   2u2/t2) v2* + (f3   2u3/t2) v3* +
[11 v1*]/x1  11 v1*/x1 + [21 v2*]/x1  21 v2*/x1 + [31 v3*]/x1  31 v3*/x1 +
[12 v1*]/x2  12 v1*/x2 + [22 v2*]/x2  22 v2*/x2 + [32 v3*]/x2  32 v3*/x2 +
[13 v1*]/x3  13 v1*/x3 + [23 v2*]/x3  23 v2*/x3 + [33 v3*]/x3  33 v3*/x3}
dx dx dx =
 {(f   ) . v*}  {11 v1*/x1  21 v2*/x1  31 v3*/x1  12 v1*/x2 
22 v2*/x2  32 v3*/x2 + 13 v1*/x3  23 v2*/x3  33 v3*/x3} dx dx dx
+  (11 + 12 + 13) v1*x1 dx dx +
 (21 + 22 + 23) v2*x2 dx dx +  (31 + 32 + 33) v3*x3 dx dx= 0

Si  est symétrique (toujours sauf cas très particulier de couplage mécano-magnétique) :

 {(f   ) . v*} dx dx dx   { : grad v*} dx dx dx
+  (11 + 12 + 13) v1*x1 dx dx +
 (21 + 22 + 23) v2*x2 dx dx +  (31 + 32 + 33) v3*x3 dx dx= 0

En remarquant que, n, la normale sortante du milieu, est telle que ni = dxj dxk / d (i, j et k tous
différents, voir l'annexe 1) :
 {(f   ) . v*} dx dx dx   { : grad v*} dx dx dx
+  (11 + 12 + 13) v1*x1 n1 d +
 (21 + 22 + 23) v2*x2 n d +  (31 + 32 + 33) v3*x3 n d =

11
  {(f   ) . v*} d   { : grad v*} d
+  V1 v1* ds +  V2 v2* ds +  V3 v3*d =

Théorème des puissances virtuelles :

 {(f   ) . v*} d   { : grad v*} d +  {V . v*} d = 0

En d'autres termes, la puissance des efforts intérieurs est égale à la puissance des efforts extérieurs.

Dans le cas statique, on remplace le vecteur v* vitesse par le vecteur déplacement u* et le terme
  est donc considéré comme négligeable. On aboutit alors au

Théorème des travaux virtuels :

 {f . u*} d   { : grad u*} d +  {V . u*} ds = 0

En d'autres termes, le travail des efforts intérieurs est égal au travail des efforts extérieurs.

NB. On a postulé ici le principe de l'équilibre et démontré le théorème des puissances virtuelles.
Si l'on postule le principe des puissances virtuelles, on démontre alors le théorème de l'équilibre.
Ces deux notions sont duales.

1.8. Solutions en MMC

Un champ partout statiquement et cinématiquement admissible qui respecte partout la loi de
comportement est solution exacte d'un problème de M. M. C. . Dans le cadre des h. p. p. – pas de
non-linéarité matériau ni géométrique -, il solution unique et exacte du problème.

2. Comportement des matériaux – résistance des matériaux

2.1. Classification du comportement mécanique des matériaux
2.1.1. Préambule
Le comportement mécanique d’un matériau à l’état solide peut présenter différents types de
composantes. La plus connue est l’élasticité. Elle se manifeste pour la plupart des matériaux en
petites déformations – typiquement moins de 1 ‰ pour les matériaux raides tels que les aciers,
céramiques et moins de 1 % pour les matériaux plus souples tels que le bois et les matières
plastiques2. Dès que l’on soumet un matériau à de plus grandes déformations, son comportement
mécanique peut présenter de la plasticité et/ou de l’endommagement (inélasticité). L’éventuelle
viscosité du matériau peut en revanche se manifester dès les petites déformations. Cette
composante du comportement est responsable du fluage, de la relaxation et de l’élasticité
recouvrée (anélasticité). Dès que l’on sort du cadre de l’élasticité linéaire on parle de " non-
linéarité du comportement mécanique ". On peut ainsi traiter de viscoélasticité, viscoplasticité,
viscoélastoplasticité, plasticité avec endommagement, élasticité non linéaire, etc.
La composante visqueuse du comportement est responsable de la sensibilité de la réponse
mécanique à la vitesse de la sollicitation3 ou au temps d’application de la sollicitation mécanique.
2
dans le langage usuel matières plastiques désigne les hauts polymères. En toute rigueur il conviendrait de
considérer que "matières plastiques" désignent les matériaux déformables présentant un comportement plastique et les
métaux entreraient dans cette appellation.
3
La sollicitation mécanique peut être une déformation imposée ou une contrainte imposée, voire une
combinaison des deux.

12
Le temps et la température sont liés, souvent de façon complexe. D’une façon générale, on
constate expérimentalement que plus le matériau solide est proche de sa température de
ramollissement ou fusion, plus les effets de sa viscosité sont importants. Ainsi, à température
ambiante un bitume ou un haut polymère seront bien plus sensibles à la vitesse de sollicitation
qu’un acier ou une céramique. Pour les aciers, on admet en général que les effets de la viscosité à
température ambiante sont du second ordre. L'eau liée dans les bétons et le bois est aussi
responsable d'une relativement faible viscosité.

2.1.2. Définitions
Elasticité : comportement réversible sans dissipation d’énergie, i. e. l’énergie fournie pour
déformer le matériau est totalement restituée au déchargement. On peut considérer une relation
contrainte-déformation linéaire aux petites déformations.
contrainte

décharge
charge
déformation
Comportement élastique linéaire – représenté schématiquement par un ressort.

Une conséquence importante d'un tel comportement mécanique est le principe de

superposition valide sous les H. P. P. en élasticité linéaire. En effet, les contraintes étant
proportionnelles aux déformations, la superposition de sollicitations conduit à la superposition de
leurs effets.
Pour les caoutchoucs et élastomères le comportement élastique peut être non linéaire aux
grandes déformations – typiquement au-delà de 20 % de déformation.
contrainte

décharge
charge

déformation
Comportement élastique non linéaire.

Plasticité : comportement irréversible avec dissipation d’énergie mécanique ne dépendant pas de la

vitesse de sollicitation ou de sa durée d’application. Il y a une déformation résiduelle lorsque la
sollicitation mécanique est supprimée.

charge
contrainte

décharge
déformation
Comportement plastique parfait – représenté schématiquement par un patin (avec frottement).

13
 contrainte
charge

décharge

déformation
Comportement élastique plastique parfait.

Viscosité : comportement de type fluide avec dissipation d’énergie mécanique. L’application

d’une contrainte induit une vitesse de déformation et réciproquement. En d’autres termes, si la
déformation imposée ne varie plus au cours du temps, toute contrainte tend à disparaître, avec
déformation irréversible.
contrainte

vitesse de déformation
Comportement visqueux linéaire représenté schématiquement par un amortisseur.

Viscoélasticité : comportement réversible avec dissipation d’énergie mécanique dépendant de la vitesse de

sollicitation ou de sa durée d’application. Il n’y a pas de déformation résiduelle lorsque la sollicitation
mécanique est supprimée.
contrainte

charge

décharge
déformation
Comportement viscoélastique.

Viscoplasticité : comportement partiellement irréversible avec dissipation d’énergie mécanique dépendant

de la vitesse de sollicitation ou de sa durée d’application. Il y a une déformation résiduelle lorsque la
sollicitation mécanique est supprimée. Cette déformation résiduelle dépend de l’histoire du chargement.

14
 350
AU4G trempé à 200 °C acier inox. 316 à 20°C
300 220

250
200

 (MPa)
 (MPa)

200

150
180
100 t=0 = 200 MPa
50 t=0 = 230 MPa
160 t=0 = 250 MPa
0
1 10 100 3 4 5
10 10 10
t (s) t (s)
Exemples de résultats d'essais de relaxation sur métaux.

Endommagement : la sollicitation mécanique engendre une variation des modules d’élasticité, voire de
l’anisotropie4. C’est typiquement une des composantes du comportement mécanique du béton hydraulique
en compression lorsqu'on dépasse la contrainte réglementaire.
contrainte

charge

recharge
décharge
déformation
Comportement endommageable.

Effets de la température : Les propriétés mécaniques évoluent avec la température. Selon les mécanismes
sous jacents impliqués et interagissant, les corrélations entre les évolutions peuvent être plus ou moins
complexes.

600 30
500 200
500
400
400 20
Fe, Re (MPa)
Fe, Re (MPa)

Ar (%)
E (GPa)

300
300
150
200 200 10

100 100

0 100
0 0
0 100 200 300 400 500 600
0 100 200 300 400 500 600
température (°C) température (°C)
Variations des caractéristiques de deux aciers de type acier doux avec la température.

4
Anisotropie : opposé à l’isotropie, l’anisotropie se manifeste par une réponse mécanique du matériau qui
dépend de la direction de sollicitation.

15
 2.1.3. Modèles rhéologiques
Ce sont des modèles schématiques simples pour la compréhension et l’analyse pédagogique du
comportement viscoélastique.

F = k x  F / S = E (L / L)


F=kv  F/S= (L / L)
t

Modèle rhéologique : élément de Maxwell pour le comportement viscoplastique.

Modèle rhéologique : élément de Kelvin-Voigt pour le comportement viscoélastique.

2.1.4. Plasticité dans les métaux : micromécanismes de base et "comportement

émergent"
Nombreux sont les phénomènes physiques qui nécessitent d'imbriquer les lois et
mécanismes à plusieurs échelles pour comprendre le comportement macroscopique. On parle de
comportement émergent. En effet, par exemple, la connaissance précise des interactions force-
distance entre les atomes à l'échelle nanométrique ne permet pas de prédire le comportement
plastique ni la rupture des métaux. En effet, des mouvements coopératifs à une échelle supérieure
rendent possible l'accommodation de la déformation imposée à bien moindre énergie potentielle
que des mouvements de blocs ou l'éloignement des atomes. Il s'agit en l'occurrence des
mouvements de défauts appelés dislocations.

16
 Schéma de principe montrant le déplacement d'une dislocation coin dans un cristal. Une
dislocation préexistante traverse une région du cristal pour laisser une déformation plastique
résiduelle (déformation permanente).

2.2. Notion de thermomécanique

2.2.1. Coefficient de dilatation thermique
Le coefficient de dilatation thermique pour les solides est aussi appelé coefficient de
dilatation linéique. Il relie les variations dimensionnelles aux variations de température. Supposons
un matériau dont les propriétés sont isotropes et un écart de température V suffisamment faible,
on peut alors écrire :
L = V L V = V L (Va – V0)
où L est la variation de longueur, V le coefficient de dilatation linéique, L la longueur initiale (à
V0), V0 la température initiale et Va la température actuelle. La longueur à la température Va est
L + L. On entend par écart de température suffisamment faible, un écart de température tel que
V puisse être considéré comme constant (pas de changement de phase…) et que L << L.

Matériau Module d’Young Coef. de Densité Coef. De Conductivité Température de

macroscopique Poisson (kg / m3) dilatation thermique fusion
(GPa) linéaire (W m1 °C1) (°C)
(°C1)
Al 68 0.33 2700 2.5 105 210 660
Cu 110 0.34 8960 1.65 105 385 1083
Fe 200 0.29 7870 1.2 105 76 1535
Ti 116 0.34 4500 0.9 105 17 1660
Zn 90 0.25 7100 3.1 105 110 420
Pb 14 0.42 11340 2.9 105 33 327
Sn 42 0.33 5765 2.4 105 63 232
Béton 20 à 40 0.18 2500 0.7 à 1.3 105 1.75  1500
hydraulique
Bois Isotrope tranverse Isot. tr. 400 à 1000  0.4 105 -
// fibres  10  0.45 % humid.! 0.12 à 0.23
Polymères <3 0.3 à 0.5 900 à 1600 5 à 20 105 0.1 à 0.4 100 à 250
Quelques caractéristiques physiques de matériaux à 20°C.

2.3. Elasticité linéaire – loi de Hooke

2.3.1. Elasticité isotrope
Lorsque la réponse à une sollicitation mécanique est indépendante de la direction de
sollicitation, on dit que le comportement est isotrope. Dans le cas de l'élasticité linéaire isotrope,
on peut relier les contraintes aux déformations à l'aide de deux paramètres indépendants du
matériau, couple de paramètres à choisir selon la commodité parmi :

17
 - module d'Young E et coefficient de Poisson  ( < 0.5),
- coefficients de Lamé,  et µ,
- modules de compressibilité K et de cisaillement G (G = µ).
Les relations suivantes relient ces paramètres :
E E
G=µ= K=
2 (1 + ) 3 (1  2)
 µ (3 + 2µ)
= E=
2 ( + µ) ( + µ)
E G
=  = 0.5  (2K + 2G/3)
(1 + ) (1  2)
Lorsque le matériau est beaucoup souple en cisaillement que sous contrainte hydrostatique alors on
dit qu'il se déforme à volume quasi-constant et   0.5. La loi de Hooke est telle que :
ij = 2µ ij + ij  (11 + 22 + 33)   = 2µ  +  Tr() Id
ij = (1 + )/E ij – ij /E (11 + 22 + 33)   = (1 + )/E  – /E Tr() Id

2.3.2. Thermo-élasticité isotrope

La déformation liée à la sollicitation mécanique s’additionne à celle d’origine thermique. Si
LT est le coefficient linéique de dilatation thermique à la température V, alors il vient :
méca. total

 =  LT V ij = (1 + )/E ij – ij /E (11 + 22 + 33)

ij ij
ij = 2µ ijméca + ij  (11méca + 22méca + 33méca) =
ij = 2µ (ijtotal  ijthermique) + ij  ((11total  11thermique) + (22total  22thermique) + (33total  33thermique))

ij = 2µ (ijtotal  LT T ij) + ij  (11total + 22total + 33total  3LT V)

NB. Le coefficient volumique de dilatation thermique VT vaut 3LT.

2.3.3. Elasticité anisotrope

Lorsque la réponse à une sollicitation mécanique est dépendante de la direction de
sollicitation, on dit que le comportement est anisotrope. Dans le cas du bois, il s'agit d'un
comportement orthotrope – 3 directions à considérer - et il y a neuf paramètres indépendants pour
définir l'élasticité. On peut éventuellement simplifier en moyennant le comportement
perpendiculairement aux fibres – 2 directions à considérer – ce qui revient à considérer que le
comportement est isotrope transverse et il ne reste que cinq paramètres indépendants.

2.3.4. Limite élastique et critères de plasticité

La limite élastique est souvent mesurée en traction uniaxiale. Ceci correspond, par
exemple, au cas particulier de l’état de contrainte dans une poutre fléchie en absence d’effort
tranchant significatif. Dans le cas général, il faut considérer que le matériau est dans un état de
contrainte triaxial et être capable de définir sa limite élastique. Pour les métaux, par exemple, les
critères de plasticité sont indépendants de la pression hydrostatique, i. e. seuls les cisaillements
interviennent dans le critère5. Chaque fois que l'on peut, si cela suffit convenablement, on utilise
commodément une variable scalaire issue du tenseur des contraintes pour définir la limite
5
Cette quasi-insensibilité à la pression hydrostatique n’est généralement pas le cas des sols ou des polymères
dont les critères sont dits piézo-dépendants.

18
élastique. Les critères les plus classiquement utilisés sont celui de Tresca, associé à la plus grande
contrainte de cisaillement, et celui de von Mises, associé à l’énergie de cisaillement. Dans le
repère principal des contraintes ces critères s’expriment :
Critère de Tresca : (1 – 3) / 2 < fe  élasticité
2
Critère de von Mises : 1/2{(1 – 2)2 + (2 – 3)2 + (1 – 3)2} < fe  élasticité
où fe désigne la limite élastique en traction et 1  2  3 désignent les contraintes dans le repère
principal.
Critère de Mohr- Coulomb : n  µn + c  élasticité
où µ = tan  désigne un frottement interne et c une cohésion. Ce critère est parfois utilisé pour les
sols à la limite des formalismes de plasticité et rupture, voir TD. Ce critère fait intervenir la
pression.

Exemple. Considérons une poutre fléchie par un moment fléchissant M3G et soumise à un effort
tranchant V(x1) et à un effort normal N(x1). En utilisant les notations usuelles du cours de
résistance des matériaux, si la théorie des poutres s’applique, on peut écrire que :
N(x1) M3G(x1) x2
= A  I3G
et que :
V(x1) S(x2, v)
 = b(x ) I
2 3G
à une altitude x2 dans la poutre à une abscisse x1. On calcule ensuite les contraintes principales qui
2 2
valent 1 = 1/2{ + 2 + 2}, 2 = 0 et 3 = 1/2{  2 + 2} dans ce cas particulier où les
contraintes autres que  et  sont nulles. Enfin, on calcule la combinaison de ces contraintes qui
correspond au critère convenable pour le métal utilisé (von Mises le plus souvent pour les aciers)
que l’on compare à la donnée matériau fe, mesurée en laboratoire. Cet exemple permet de
comprendre les règles de calcul aux ELU (couplages entre MELU, VELU et NELU) en construction
métallique.

2.4. Rupture
2.4.1. Rupture sous chargement monotone
On ne peut mesurer le comportement ultime d'un matériau qu'en allant jusqu'à la rupture.
Lorsque, par exemple, on mesure l'allongement à rupture d'une éprouvette de traction, les résultats
d'essai peuvent s'avérer peut reproductibles, d'autant plus que les matériaux sont fragiles. En effet,
dans les solides, une fissure finit par apparaître en s'amorçant sur un défaut soit intrinsèque à la
matière (impureté, cavité, …), soit extrinsèque (usinage, rayure, …). Le comportement ultime de
la matière en sommet de fissure n'est donc mesurable qu'au travers de la réponse d'une structure
fissurée. Pour tenter de s'affranchir des défauts préexistants, on peut en générer un de taille
contrôlée en procédant à une entaille dans l'échantillon. Ceci permet de mesurer la ténacité du
matériau. Toutefois, l'étude de la rupture sous chargement monotone demeure une affaire de
spécialistes.

2.4.2. Rupture sous chargement cyclique, fatigue

On désigne par fatigue l’endommagement des matériaux sous sollicitations mécaniques
cycliques, alternées ou simplement variables. En laboratoire, on caractérise les matériaux à la
fatigue presque toujours avec des sollicitations cycliques6. Lorsque le matériau est sollicité au-delà

6
On remarque que, la plupart du temps, comme dans le cas des séismes par exemple, ce ne sont pas les
données matériaux qui manquent mais le spectre de la sollicitation qui est méconnu.

19
de sa limite élastique, le nombre de cycles à rupture est faible et l’on parle de fatigue
oligocyclique. Certains matériaux admettent une limite d’endurance, c’est à dire qu’il existe une
sollicitation en dessous de laquelle la rupture n’est jamais observée. La rupture, lorsqu’elle
intervient, s’amorce sur des défauts intrinsèques (microfissure, crique…) ou extrinsèques (rayure,
poinçonnement…) qui croissent jusqu’à la ruine de la structure. On représente souvent la tenue en
fatigue à l’aide de courbe de Wöhler.

400
acier XC10

350
 (MPa)

300

limite
250 d'endurance

200
1 2 3 4 5 6
10 10 10 10 10 10
nombre de cycles à rupture
Courbe de Wöhler illustrant des résultats d'essais de fatigue.

20
3. Exemples de réalisations

21
22
23
24
4. Modélisation des conditions aux limites
4.1. Déplacements imposés

en 3D : (a) – encastrement, (b) – appui glissant, (c) – appui glissant

poutres et structures : (a) – encastrement, (b) – appui simple, (c) – appui glissant

Encastrement : rotation et déplacements nuls, 3 composantes R1, R2 et C3 de réaction d'appui à

déterminer.
Appui simple : déplacements nuls, 2 composantes R1 et R2 de la réaction d'appui à déterminer.
Appui glissant : 1 composante de déplacement nulle, 1 composante R2 de la réaction d'appui à
déterminer.

4.2. Pression, contraintes imposées, surface libre

Lorsque rien n'est spécifié, la surface est libre de contrainte, i. e. les composantes normales
 et tangentielles  du vecteur contrainte sont nulles. Lorsque seule la composante normale est non
nulle, il s'agit d'une pression.

en 3D : (d) – densité surfacique de forces.

poutres et structures : (a) – charge ponctuelle, (b) – charge uniformément répartie

25
 Représentation d'une charge de neige
espace dessin, représentation du poids propre, de la neige et du vent

4.3. Liaisons extérieures

Les trois principaux type de liaison extérieure pour les problèmes plans sont
- l’appui glissant qui empêche un déplacement – selon y sur la figure suivante - ,
- l’appui fixe qui empêche les deux composantes de déplacements dans le plan,
- l’encastrement glissant qui empêche les deux composantes de déplacements dans le
plan et la rotation – selon z sur la figure suivante - ,

De gauche à droite, appui glissant, appui fixe et encastrement.

Réalisation pratique d’un appui glissant pour un pont métallique de chemin de fer.

Les liaisons peuvent être idéales comme par exemple celles du Centre Pompidou de Paris
ou plus approximatives et bien plus économiques comme dans la plupart des cas pratiques.

26
 Appuis fixes ou rotules en construction métallique réalisations idéales et réalisation économique
courante.

4.4. Appuis élastiques

En toute rigueur, aucune liaison n'est parfaite. Toutefois, on peut souvent – cas isostatique
par exemple - procéder en deux étapes si nécessaire, c'est à dire dissocier la déformation des
liaisons du calcul des déformations dans la structure.

poutres et structures : (a) – appuis élastiques, (b) – encastrement élastique, (c) – massif élastique

4.5. Liaisons intérieures

Les liaisons intérieures reprennent les mêmes principes que les liaisons extérieures, c'est à
dire transmission ou non des deux composantes d'une force et/ou d'un couple. Elles sont souvent
approximatives d'un point de vue constructif. La modélisation simplifie la réalité.

4.5.1. Cas de la construction métallique

L’effort tranchant est transmis essentiellement par l’âme des profilés en I. Le moment
fléchissant est transmis essentiellement par les semelles. La rigidité en torsion est relativement
faible et on évite en général de solliciter en torsion.

27
 Constructions métalliques.
espace dessin, liaison âme et liaison semelles

4.5.2. Cas du béton armé et/ou précontraint

Le dimensionnement des éléments en béton armé et/ou précontraint nécessite un
approfondissement supplémentaire à ce cours. En effet, il s’agit d’un matériau composite dont la
plus grande longueur interne est proche d’une ou des dimensions des sections.

Ferraillage. Béton armé préfabriqué.

espace dessin, continuité ou non des armatures

28
 4.5.3. Cas de la construction en bois – lamellé collé
L’effort tranchant est transmis majoritairement par le centre des sections. Le moment
fléchissant est transmis majoritairement par les extrémités. La rigidité en torsion est relativement
faible et on évite en général de solliciter en torsion.

Lamellé collé
espace dessin, transmission des contraintes par connecteurs

29
 4.5.4. Cas de la construction mixte

Ponts à haubans Construction mixte acier-béton.

4.6. Remarques
Le choix des liaisons est en général fait par l’ingénieur des structures en bureau d’étude, en
respect des volontés de l’architecte. Ce choix repose sur plusieurs critères dont essentiellement le
la sécurité et le prix. Par exemple, on se prémunira au mieux des dilatations thermiques
différentielles, d’un éventuel tassement différentiel des appuis, du report des charges en cas de
ruine d’un des éléments de structure. Des dispositions parasismiques sont prévues dans les
réglementations. L’utilisation optimale de la matière conduira aussi le choix des liaisons. Toutefois
il faut de plus compter dans le prix de la structure les coûts des assemblages et des phases de
construction.

5. Spécificités de la théorie des poutres

5.1. Poutre – barre
5.1.1. Définition
Poutre : objet géométrique généré par une section perpendiculaire à une fibre moyenne, les
dimensions de la section étant petites devant la longueur de la fibre. Le rayon de courbure doit être
grand relativement à la longueur de la fibre. La poutre est dite droite si la fibre moyenne est
rectiligne.
Barre : poutre droite non fléchie.

Pour la suite, on se limitera aux poutres droites et aux structures bidimensionnelles.

5.1.2. Equilibre des poutres et barres

Moment d'une force appliquée en P dans le plan (x1, x2) par rapport à O : M3O = OP  FP
Il suffit d'appliquer le V. P. V. avec un mouvement de solide (translation, rotation) à une structure
pour se persuader que l'équilibre sera assuré si :

30
  Fext = 0
 MO(Fext) = 0, pour un point O quelconque.
Dans un plan, cela fournit 2 + 1 = 3 équations, dans l'espace 3 + 3 = 6 équations. Ces équations
découlent du P. V. V. donc du principe div  = 0.

Démonstration : supposons les forces inertielles négligeables. Soit v un champ de vitesses virtuel
uniforme quelconque. Dans un mouvement de solide rigide, la puissance des efforts extérieurs est
nulle, i. e. Pext =  Fi . v dt = 0. Comme v est uniforme, on peut écrire Pext = v dt .  Fi = 0. Si
Pext = 0 quel que soit v non nul alors  Fi = 0. Soit  une vitesse de rotation virtuelle quelconque
autour d'un point O quelconque. La vitesse d'un point Pi vaut (  OPi). De même
Pext =  Fi . (  OPi) dt = 0 en absence de déformation virtuelle. Le produit mixte est invariant
par permutation circulaire donc Pext =   . (OPi  Fi) dt =  dt .  (OPi  Fi). Si Pext = 0 quel
que soit  non nul alors  (OPi  Fi) = 0.

Si les équations d'équilibre permettent de tout (réactions d'appui, déplacements, contraintes et

déformations) déterminer partout dans la structure, alors elle est dite isostatique. Si ces équations
d'équilibre global ne suffisent pas, elle est dite hyperstatique.

5.2. Hypothèse de Navier-Bernoulli

5.2.1. Cinématique
Les sections droites restent droites. Cette hypothèse fournit une approximation de qualité
d'autant plus excellente que la hauteur de la poutre est petite par rapport à sa longueur. On peut
aisément vérifier expérimentalement ou par modélisation numérique que cette hypothèse est très
proche de la réalité, même lorsqu'on dépasse modérément la limite de linéarité du comportement.

Cinématiques possibles si les section droites restent droites.

espace dessin, allongement – raccourcissement des fibres

5.2.2. Déformation d'une section

Le rayon de courbure R est tel que dx1 = R sin d  R d i. e. 1 / R = d/dx1. Le déplacement en
(x1 + dx1, x2) vaut donc (x2 sin d, R  R cos d soit, pour des grandes courbures
(x2 d, R ( 1  (1  d(x2 d, dx1 d(x2 d, 0 au premier ordre. L'allongement ou

31
raccourcissement d'une fibre située en x2 vaut donc 11 = x2 d/dx1. Les déformations
longitudinales varient linéairement dans la section.
espace dessin et calculs, géométrie en flexion circulaire

Distribution des contraintes normales en flexion pure. La fibre neutre (contrainte normale nulle)
est confondue avec la fibre moyenne.

Exemple de distribution des contraintes normales en flexion composée. La fibre neutre n'est plus
confondue avec la fibre moyenne.

N(x1) M3G(x1) x2
(x1, x2) = A(x )  I (x )
1 3G 1

Où N est positif en traction. En notation (x, y, z) avec A et I3G constants et avec la notation
simplifiée de la théorie des poutres pour les contraintes, on réécrit :
N(x) M(x) y
(x) = A  I
Gz

La contrainte maximale en valeur absolue se situera soit en ymin, noté v', soit en ymax noté v, selon
les signes de N et M.
N M v' N Mv
max = sup{ A  I  , A  I }
Gz Gz

5.2.3. Admissibilité
L'hypothèse des sections droites est cinématiquement admissible. La loi de comportement
du matériau est respectée. Si les contraintes étaient statiquement admissibles, alors sous les
H. P. P. en élasticité linéaire, on obtiendrait ainsi l'unique et exacte solution du problème, loin des
points d'application des charges. Mais ce n'est pas le cas, la solution proposée n'est pas

32
statiquement admissible. Pour s'en convaincre, il faut étudier de plus près le cisaillement. On se
reportera aux annexes pour plus de détails.
espace dessin

Gauchissement sous l'effet du cisaillement, admissibilité statique et réciprocité des contraintes

tangentielles
espace schéma, surfaces sup. et inf. libres de contrainte

5.3. Principe de St Venant, torseur généralisé

5.3.1. Equilibre d'une section
Concernant l’équilibre d’une structure soit on postule le principe l’équilibre des forces et
des moments directement sur la structure, soit on postule le principe des puissances virtuelles
(P. P. V.) et on démontre le théorème de l’équilibre des forces et des moments sur la structure.
Pour cette dernière démonstration, il suffit d’imposer des champs « rigidifiants » de translations et
rotations7.
Pour des raisons évidentes, la théorie des poutres ne peut donner de résultats satisfaisants
au voisinage des points d'application des forces. Par contre, loin des points d'application des
forces, l'action sur une section de poutre peut se résumer à un torseur résultant (effort Normal,
effort tranchant, Moment fléchissant). Pour une structure plane, ce torseur pourra être calculé en
isolant la partie droite - ou gauche selon la convention arbitraire que l'on se fixe – à partir de la
section dans laquelle on désire connaître la répartition des contraintes. En utilisant la cinématique
simplifiée des sections qui restent droites, on peut, donc en utilisant les conditions d'équilibre et la
loi de comportement retrouver les équations de la théorie des poutres. On rappelle que les H. P. P.
sont faites et que le comportement du matériau est élastique linéaire. Pour plus de commodité, on
se situera au niveau de la fibre moyenne de la section de la poutre située en G d'ordonnée x2G.
Si l'on considère une section quelconque, l'action étant égale et opposée à la réaction en vertu de
l'équilibre, on peut écrire :

7
Notons que l’on a choisi de démontrer le théorème des puissances virtuelles (T. P. V.) à partir du principe de
l’équilibre des forces et des moments sur un élément de volume.

33
  Fi à gauche +  Fi à droite = 0
 GPi  (Fi à gauche) +  GPi  (Fi à droite) = 0
En isolant la partie gauche de la structure, on peut écrire l'équilibre ainsi :
 Fi à gauche + N(x1) x1 + V(x1) x2 = 0
 GPi  (Fi à gauche) + M3G(x1) = 0
où N représente la force résultante des contraintes normales, V celle des contraintes tangentielles
et M le couple en G résultant des contraintes normales. Défini ainsi, le torseur résultant (N, V, M)
correspond donc au torseur résultant des forces extérieures à droite qui doit équilibrer celui de
gauche impliqué dans les équations ci-avant.
espace dessin, sous-structure isolée et efforts intérieurs reportés

5.3.2. Relation entre l'effort tranchant et le moment fléchissant

M3G(x1) / x1 = V(x1)
espace dessin et calculs, section en équilibre

5.3.3. Répartition des contraintes dans une section

5.3.3.1. Contraintes normales
On considérera que 13 = 23 = 0 et aussi que 22 = 0 (en rapport avec l'hyp. de SV Venant).
Les intégrales des contraintes réparties dans la section doivent correspondre à ces résultantes.

u(x1) / x1 = N(x1) / E(x1) A(x1)

(x1) / x1 = M3G(x1) / E(x1) IG3(x1)

(x1, x2) = x2 M3G(x1) / IG3(x1) + N(x1) / A

espace dessin et calculs, superposition des contraintes et déformations

34
5.3.3.2. Moment quadratique et théorème d'Huygens
vsup.
IG3 = A x22 x22 b(x2) dx2
dx2 dx3 = 
vinf.
vsup.
où G est la position de la fibre moyenne telle que 
x2 b(x2) dx2 = 0. Ainsi défini, toute rotation
vinf.
différentielle /x1 autour de G ne génère d'effort normal. Si l'origine O du repère n'est pas en
fibre moyenne, on peut calculer ainsi le moment quadratique de A par rapport à O :
IO3 = IG3 + ║OP║² A
vsup.
b(x2) dx2.
où A est l'aire de la section A, i. e. A = 
vinf.
vsup.
Démonstration : x22 b(x2) dx2
IO3 = 
vinf.
vsup. vsup.
IO3 = 
(x2  xG2) b(x2) dx2 + 
2
(2 x2 xG2  xG2 ) b(x2) dx2
2

vinf. vinf.
vsup. vsup.
IO3 = IG3 + xG2 { 2 (x2  xG2) b(x2) dx2  xG2 b(x2) dx2}
vinf. vinf.
IO3 = IG3 + xG22 A
espace dessin et calculs, intégrations sur une section

5.3.3.3. Effort normal et Moment fléchissant

En béton armé par exemple, il peut parfois être commode de calculer le torseur résultant en
un point O, choisi comme origine des x2, qui n'est pas situé en la fibre moyenne. En intégrant les
contraintes dans la section, il apparaît que :
NP(x1) = NG(x1)
MP(x1) = MG(x1)  xG2 NG(x1)
espace dessin et calculs, exemple pour une section de béton armé

35
5.3.3.4. Contraintes tangentielles
(x1, x2) = V(x1) S(x1, x2) / IG3(x1) b(x1, x2)
espace dessin et calculs, équilibre ∂M/∂x1 et fibres supérieures entre x2 et vsup.
5.3.3.5. Diagrammes de sollicitations
Poutre sur deux appuis avec une force F. –x2 appliquée en son milieu :
Equilibre : RA1 = 0, RA2 + RB2  F = 0, RB2 l  F l/2= 0
NAB(x1) = 0
VAC(x1) = F / 2, pour x1  [0, l/2] VCB(x1) = F / 2, pour x1  [l/2, l]
MAC(x1) = F x1 / 2, pour x1  [0, l/2] MCB(x1) = F (l  x1) / 2, pour x1  [l/2, l]
espace dessin et calculs, diagrammes N, V & M

5.3.4. Déformation d'une poutre

Pour calculer la déformée d’une structure, il faut intégrer les déformations le long des
poutres et identifier les constantes d’intégrations à l’aide des conditions de liaison aux extrémités
de chaque poutre. En pratique, il n’est pas nécessaire de prendre en compte les déformations
d’effort normal et tranchant dans les éléments fléchis.
Déformation liée à l’effort normal : u1 / x1 = 11(x1) = N(x1) / E(x1) A(x1) noté (x1)
Déformation liée à l’effort tranchant : u2 / x1 = 12(x1) = V(x1) / µ(x1) Sr(x1) noté (x1)
Déformation liée au moment fléchissant : y"(x1) =  / x1 = 1 / R(x1)
E(x1) IG3(x1) y"(x1) = M3G(x1)

5.3.5. Remarque
La théorie des poutres permet de réduire d’une dimension les problèmes à résoudre. Plus
exactement, les variables x1 et x2 sont dissociées et interviennent successivement. En effet, on
commence par résoudre des problèmes dont la variable est x1 le long de la fibre moyenne. Ensuite,
pour connaître la distribution des contraintes dans une section située en x1 donné, la variable est x2.
Ceci fournit des résultats très satisfaisants « loin » - typiquement on écarte la section de (x1  h) à
(x1  h) où h est la hauteur de la section de poutre - des points d’application des forces. Dans ces
sections, les problèmes de répartitions des contraintes sont 2-D, voire 3-D. Ces problèmes
« d’assemblage » sont documentés par ailleurs et traités ensuite.

5.4. Energie de déformation

5.4.1. Energie et élasticité
L'énergie de déformation dans le cas de l'élasticité linéaire est analogue à celle d'un ressort.
Dans le cas du ressort la force F est proportionnelle au déplacement relatif x des extrémités d'où
F = k x, k étant la constante de rigidité. L'énergie élastique (énergie libre) stockée dans le ressort

36
vaut donc W = ∫ F dx = ∫ k x dx = ½ k x² = ½ F x . En considérant que tous les efforts sont
appliqués en même temps avec la même progression - étant donné que l'histoire du chargement
n'influe pas sur l'état actuel dans le cas des H. P. P. en élasticité linéaire, il est loisible de choisir le
cas le plus simple pour cette démonstration - l'énergie de déformation d'une structure peut donc se
mettre sous la forme :
k k
 d1   f1 
 k  k
n n 3
W = ½   d2  .  f2  = ½   dik fik
k = 1  d3   f3 
k k
k=1 i=1

Où d1 et d2 représentent des déplacements, d3 une rotation, f1 et f2 représentent des forces, f3 un

couple.
n 3
Facteurs d'influence ijkl : dik =    ijkl fjl
l=1 j=1
n n 3 3
W = ½       ijkl fjl fik
k=1 l=1 i=1 j=1
théorème de Maxwell-Betti : ijkl = jikl
n 3
W / fik =    ijkl fjl = dik
l=1 i=1
théorème de Castigliano : W / fik = dik
espace dessin et calculs, exemple simple sur une poutre isostatique

5.4.2. Energie de déformation liée à l’effort normal

V = 0, M = 0, dWx1→x1 + dx1 = N(x1) du*(x1)

N(x1) N²(x1)
En élasticité du* est proportionnel à N, du*(x1) = E A dx1 donc dWx1→x1 + dx1 = ½ E A dx1
Par E A, il faut bien sûr entendre E(x1) A(x1)
t = ta.
 
 t
(t) dt N(x1)
car dWx1→x1 + dx1 =  (t) N(x1) dx1
 EA
t=0
(t) définit une progression continue et dérivable quelconque de chargement dans le temps telle
que (t0) = 0 et (ta) = 1.
t = ta.
t = ta. ta
  
et car  (t) (t) dt =   (t) '(t) dt = ½ │ = ½[(ta)  (t0)] = ½
 t
t=0 t0
t=0

espace dessin et calculs, déformations et contraintes d'effort normal

37
 5.4.3. Energie de déformation liée à l’effort tranchant

N = 0, M = 0, dWx1→x1 + dx1 = V(x1) dv*(x1)

V(x1) V²(x1)
En élasticité dv* est proportionnel à V, dv*(x1) = G S dx1 donc dWx1→x1 + dx1 = ½ G S dx1
r r
Par G Sr, il faut bien sûr entendre G(x1) Sr(x1)

espace dessin et calculs, déformations et contraintes d'effort tranchant

5.4.4. Energie de déformation liée au moment fléchissant

M²
N = 0, V = 0, dWx1→x1 + dx1 = ½E I dx1
G3
N = 0, V = 0, dWx1→x1 + dx1 = M(x1) d*(x1)
M(x1) M²(x1)
En élasticité d* est proportionnel à M, d*(x1) = E I dx1 donc dWx1→x1 + dx1 = ½ E I dx1
GZ GZ
Par E IGZ, il faut bien sûr entendre E(x1) IGZ(x1)

espace dessin et calculs, déformations et contraintes de flexion

5.4.5. Energie de déformation élastique dans une poutre

dWx1→x1 + dx1 = ½ { EN²A + µV²Sr + EM²IG3 }dx1

l
[N(x1)]² [V(x1)]² [M(x1)]²
W0→l = ½  {
 E(x1) A(x1) + µ(x1) Sr(x1) + E(x1) IG3(x1) dx1 }
0
[V(x1)]² [M(x1)]² [N(x1)]² [M(x1)]²
en général << et <<
µ(x1) Sr(x1) E(x1) IG3(x1) E(x1) A(x1) E(x1) IG3(x1)
l
[M(x1)]²
W0→l  ½   E(x1) IG3(x1) dx1 dans les structures fléchies
0
l
[N(x1)]²
W0→l = ½  
E(x1) A(x1) dx1 dans les structures triangulées-rotulées
0
espace dessin et calculs, pas de travail de N dans  ni , …

38
6. Résolution des problèmes de poutres droites -systèmes isostatiques
plans
6.1. Equation de la déformée
Poutre sur deux appuis avec une force F. –x2 appliquée en son milieu, E IG3 est constant.
Conditions aux limites (conditions d'appuis) : y(0) = 0 et y(l) = 0
Pour x1  [0, l/2] : E IG3 yAC"(x1) = F x1 / 2, d'où :
E IG3 yAC'(x1) = F x12/4 + a
et E IG3 yAC(x1) = F x13/12 + a x1 + b
Pour x1  [l/2, l] : E IG3 yCB"(x1) = F (l  x1) / 2
E IG3 yCB'(x1) = F (l  x1)2/4 + c
et E IG3 yCB(x1) = F (l  x1)3/12  c (l  x1) + d
Pour ce problème, il y a 4 constantes d'intégrations inconnues a, b, c, d à identifier. Conditions de
raccordement : yAC'(l/2) = yCB'(l/2) et yAC(l/2) = yCB(l/2). On en déduit :
b = 0, d = 0, c = a + Fl2/8 puis a = Fl2/16 donc c = Fl2/16. Le résultat est donc :
Pour x1  [0, l/2] :
yAC'(x1) = (F x12/4  Fl2/16 ) / E IG3 et yAC(x1) = (F x13/12  Fl2 x1/16) / E IG3
Pour x1  [l/2, l] :
yCB'(x1) = (F (l  x1)2/4  Fl2/16 ) / E IG3 et yCB(x1) = (F (l  x1)3/6  Fl2 (l  x1)/16) / E IG3
La flèche maximale correspond à y'(x1) = 0 donc pour x1 = l/2 et vaut :
ymax = y(l/2) = Fl3/ 48 E IG3

AN : bois résineux de classe C24, en flexion fmax = 24 MPa, E = 11 GPa. ymax < l / 200.
Portée de 2.5 m. Section de chevron 80 x 60 mm². IG3 = 0.06 x 0.083/12 = 2.56 106 m4
Condition de résistance :
F = 800 N, Mmax = 1250 N.m, max = 19.53 MPa.
max < fmax  Mmax v / IG3 < fmax  Fl v / 4 IG3 < fmax  F < 4 IG3 fmax / l v
F < 983 N
Condition de flèche :

39
ymax = 9,25 mm
ymax = Fl3/ 48 E IG3< l / 200  F < 48 E IG3 / 200 l²
F < 1081 N
_________________________________________________________________
V
S Venant ? - Charge F = 800 N appliquée sur 60 x 60 mm².
22  800 / 36 104 m2 = 0.22 MPa << max
déformation liée à V : u2   l/2  u2  r l/2 où V = F/2 et Sr = 5A/6. Considérons µ de
l'ordre 1 GPa, on obtient u2  1.25 104 et on remarque que u2 << y(l/2)
espace dessin et calculs, déformations en parallélogramme de V négligeables, contrainte 22 sous
la charge qui tend vers 0 quand on s'éloigne de F

6.2. Déplacement ou rotation en un point

Si l'on a besoin que d'une valeur en un point et non de l'équation complète de la déformée,
on peut se tourner vers les méthodes "énergétiques" pour résoudre efficacement le problème. Le
théorème de Castigliano peut être utilisé. Toutefois, le T. T. V. permettra de résoudre plus de
problèmes sans que les calculs ne soient plus fastidieux. Lorsqu'on cherche un déplacement d (ou
une rotation) en un point P, il suffit d'appliquer une force virtuelle F* (ou un couple) au point P et
de faire travailler les efforts virtuels dans les déformations réelles. Le travail des efforts extérieurs
doit être égal au travail des efforts intérieurs. En pratique, on applique l'effort virtuel sur la même
structure mais pas forcément avec les mêmes liaisons, l'essentiel étant que seule la force (ou le
couple) travaille du côté "extérieur" et que le calcul soit le plus simple possible. Pour obtenir le
déplacement d dans la direction de F*, on écrira donc :
N*(x1) N(x1) V*(x1) T(x1) M*(x1) M(x1)
F* dcherché = 
{ E(x1) A(x1) + µ(x1) Sr(x1) + E(x1) IG3(x1) }dx1
Structure

espace dessin et calculs, exemple simple de la poutre isostatique traitée ci-avant

40
 6.3. Théorème de Vereschiaguine, intégrales de Mohr
Dans la plupart des cas l’effort virtuel varie linéairement et EI ne varie pas sur le domaine
d’intégration. On peut alors facilement montre que :

L
M(x1)M*(x1) (M) M*(x1G)

 { EI dx }
1=
EI
0
où (M) est la surface de M(x1) et x1G le centre de gravité de (M).
En effet, comme M*(x1) varie linéairement, on peut le mettre sous la forme a x1 + b. Il vient :
L L L L
M(x1)M*(x1)  

 {
EI
} dx1 = 
EI
(a x1 + b) M(x1 ) dx1 = {a 
EI  1
x M(x 1) dx1 + b 
 M(x1) dx1 }
0 0 0 0
L
 x1 M(x1) dx1
 L
 0
= {a x1G (M) + b (M)} puisque x1G  où (M) =   M(x1) dx1
EI L

 M(x1) dx1
0
0

7. Résolution des problèmes de poutres droites - systèmes

hyperstatiques plans
7.1. Degré d'hyperstaticité
Une structure en équilibre est dite isostatique si les équations d’équilibre suffisent pour
tout déterminer. Si elles ne suffisent pas, elle est dite hyperstatique.
Si une structure ne peut être équilibrée – pas assez de liaisons compte tenu des degrés de
liberté - alors il s’agit d’un système hypostatique ou mécanisme.
Le degré d’hyperstaticité est indépendant du chargement.
Lors de la conception d’une structure, le choix du degré d’hyperstaticité – des liaisons ou
des degrés de liberté – dépend des considérations de sécurité, de pérennité (tassements ou
dilatations différentiels…) et bien sûr économiques.

Pour calculer le degré d'hyperstaticité, il suffit de faire le bilan du nombre d'équations et du

nombre d'inconnues après avoir fait l'équilibre des nœuds puis des poutres. Dans le plan, il y a 3
équations d'équilibre par nœud ainsi que par poutre. Chaque extrémité de poutre connectée à un
nœud possède a priori 3 inconnues (N, V, M). Il faut ajouter les liaisons avec l'extérieur,
i. e. réactions d'appui. Les charges appliquées sont des données. Le degré d'hyperstaticité D peut
donc être calculé à l'aide de :
n noeuds
D° =  li  3 nb
i=1

où li est le résultat du bilan inconnues – équations au nœud i et nb le nombre de barres. Le degré

d'hyperstaticité est indépendant du chargement. En toute rigueur, il faut faire le calcul sur toutes
les sous structures possibles pour être certain que la structure ne présente pas la somme d'un
mécanisme local et d'une sous structure hyperstatique.

41
 dessins et calculs, équilibre des nœuds

7.2. Résolution par la méthode des forces

1 - On rend la structure isostatique en supprimant des liaisons, on obtient une structure
isostatique associée.
2 - Le système de base est la structure isostatique associée sollicitée par les efforts donnés.
3 - Les systèmes suivants sont obtenus avec la structure isostatique associée sollicitée par
un des efforts correspondant à une liaison supprimée. On procède pour tous les efforts de liaison
supprimés.
4 - On écrit que la superposition de tous ces systèmes avec le système de base reconstitue
le système de départ du point de vue des conditions de liaison.

espace dessin et calculs, exemples simples d'une poutre continue sur 3 appuis puis 4 appuis

NB. Le théorème de Ménabréa, conséquence directe du théorème de Castigliano, peut aussi être
utilisé. Celui-ci ne fait apparaître indirectement la superposition mais la superposition reste
néanmoins à la base de la résolution.

7.3. Formule des trois moments, démonstration à partir de la méthode des

forces
On considère une poutre droite continue sur 3 appuis telle que le produit du module
d’Young E et du moment quadratique IG3 soit constant, i. e. « EI = constante ». Les deux travées
peuvent être de longueurs différentes. Le chargement des travées est quelconque. Il s’agit de
déterminer le moment fléchissant MB sur l’appui central. Le système est hyperstatique de degré 1.
On choisit d’utiliser la méthode des forces et le T. T. V. pour résoudre.

42
 Poutre continue sur appuis.

On considère la poutre AB de gauche de longueur LAB. La rotation en B dans ce système

isostatique associé au système original résulte de la superposition de l’action du chargement entre
appuis, de l’action de MA et de l’action de MB. On calcule tout d’abord la rotation en B, g iso, sous
l’action du chargement extérieur. En négligeant les déformations d’effort tranchant et d’effort
normal devant celles du moments fléchissant, le T. T. V. permet d’écrire :

LAB
M(x1)M*(x1)
C* 
g iso
{
= 
E IG3
} dx1
0
où M(x1) est associé au chargement entre les appuis et M*(x1) au couple virtuel C* appliqué en B.
On calcule ensuite la rotation en B, g MA, liée à l’action de MA :
1 1 MA
C* g MA = 3C* 2 E I LAB
G3
Puis la rotation en B, g MB
, liée à l’action de MB :
2 1 MB
C* g MB = 3C* 2 E I LAB
G3

Poutre continue sur appuis, superposition et couple virtuel.

On procède de façon analogue pour la poutre BC et il vient :

LBC
M(x1)M*(x1)
C* d iso = 
 { E IG3 }
dx1
0
1 1 MC
C* d MC = 3C* 2 E I LBC
G3
2 1 M B
C* d MB = 3C* 2 E I LBC
G3

43
Il s’agit de reconstituer la structure de départ, pour cela il faut écrire la continuité en rotation de la
poutre au passage de l’appui B. Cela se traduit par :
g iso + g MA + g MB = d iso + d MC + d MB
En remplaçant par les expressions obtenues précédemment, on obtient finalement la formule des
trois moments ou formule de Clapeyron :
LAB MA + 2 (LAB + LBC) MB + LBC MC = d iso  g iso)
Bien sûr, il faut connaître d iso  g iso) pour résoudre le problème. Dans le cas d'une charge
uniformément répartie q portée par Z, on sait que d iso = qLBC3/24EI et g iso = qLAB3/24EI. On
trouvera plus loin un exemple de mise en œuvre de la méthode.

8. Méthode des déplacements

8.1. Introduction
La méthode des déplacements est basée sur la solution analytique obtenue pour une poutre
dont les extrémités sont soumises à des déplacements et rotations données8. Elle constitue une
méthode systématique de résolution qui se prête bien à la programmation informatique.
Malheureusement, sans outil informatique, elle conduit souvent à des calculs fastidieux.
Nonobstant, elle présente un grand intérêt pour l'apprentissage et peut constituer une introduction à
la méthode des éléments finis.
Sans rien perdre conceptuellement, nous nous limiterons à des structures planes dans un
plan (X, Y). Les déplacements seront notés Ui selon X, Vi selon Y et la rotation dans ce plan i
selon Z. L'indice i désigne le nœud numéro i. (Ui, Vi, i) sont appelés les degrés de liberté du
nœud9. Il n'y a donc pas de torsion autour de l'axe x de la génératrice dans les poutres. Nous
noterons en minuscules (ui, vi) les déplacements exprimés dans le repère local (x, y) de la poutre et
en majuscules ceux exprimés dans le repère global (X, Y), généralement celui de la structure.
Deux choix sont possibles pour les efforts aux extrémités des poutres, soit on considère les efforts
intérieurs (N, V, M) – ce sera le cas ici -, soit on considère les forces et couples appliqués.

Poutre orientée entre les nœuds i et j.

Les données du problème sont :

- les nœuds définis par leurs coordonnées spatiales et les types de liaisons (rotule en
extrémité de poutre, encastrement sur appui…),
- les poutres chacune définie par deux nœuds et ses caractéristiques géométriques et
matériau,
- les conditions aux limites définies par les charges imposées et les déplacements
imposés.
En utilisant la solution de l'équilibre et la déformée de la poutre calculée une fois pour toutes dans
ce cours, le principe de résolution est :
- on procède à l'écriture des équations d'équilibre de tous les nœuds de la structure,

8
Elle est aussi appelée méthode des relaxations car on peut considérer que les poutres sont chargées pendant
que les nœuds sont bloqués en déplacements et rotation puis débloqués pour minimiser l'énergie potentielle de la
structure et ainsi satisfaire aux conditions d'équilibre.
9
Dans le plan, il y a 3 degrés de liberté par nœud. Dans l'espace, il y a 6 degrés de liberté par nœud, 3
déplacements et trois rotations.

44
 - on met en forme le système linéaire d'équation, les efforts donnés sont remplacés par
leur valeur (0 par exemple pour une extrémité libre non chargée) et les efforts inconnus
sont exprimés à partir des déplacements et rotations,
- on résout le système linéaire d'équations,
- on trace les diagrammes, on calcule éventuellement des déplacements entre des nœuds.

Rappel des conventions de signes pour les efforts intérieurs, traction, cisaillement et flexion.

8.2. Solution de base

8.2.1. Composante liée à l'effort normal
Il s'agit de calculer une fois pour toutes la solution d'un problème de poutre chargée sans
connaître a priori les conditions aux limites aux extrémités. En pratique, on se donne un
déplacement ui selon x et un effort normal Nij – pour l'instant inconnus - à l'extrémité de la poutre
du côté du nœud i et on intègre les efforts et les déplacements jusqu'à l'autre extrémité j. Il faut
bien sûr respecter l'équilibre en tout point et le comportement du matériau. Les sections droites
restent droites puisque l'on utilise la théorie des poutres. On obtient ainsi le déplacement uj et
l'effort normal Nji à l'autre extrémité aux constantes de départ ui et Nij près, c'est à dire que l'on
obtient le déplacement relatif et l'effort normal relatif.
Un chargement astucieux qui permet de couvrir une large gamme de types de chargement
consiste à ne charger uniformément qu'entre l'abscisse x1 et l'abscisse x2. Donc :
x.
N(x) = Nij  
P(s) ds et P(s) = p si s  [x1, x2], P(s) = 0 si s  [x1, x2]
0
où est orienté selon x. L'effort normal dans la poutre de longueur L prend donc les valeurs
suivantes :
0  x  x1 x1  x  x2 x2  L
N(x) = Nij N(x) = Nij  p (x  x1) N(x) = Nij  p (x2  x1) = Nji

Chargement normal entre les nœuds i et j.

La poutre se déforme en allongement et/ou raccourcissement selon :

x. x.

u(x) = ui +  N(s)
N(s)
E(s)A(s)ds = ui +  EA ds
 
0 0
car on considère le cas où E(s)A(s) est constant le long de la poutre10.
0  x  x1 x1  x  x2 x2  L

10
Si cela n'est pas le cas, on pourra obtenir une approximation en procédant à une sous-structuration en créant
des nœuds intermédiaires entre i et j puis en considérant que EA est constant par morceau.

45
 Nij x. Nij p (x2  x1)²
u(x) = ui + EA x
Nij  P(s) u(x) = ui + EA x  2 EA
u(x) = ui +  ds
 EA p (x2  x1)(x  x1)
0  EA
Nij p (x  x1)²
= ui + EA x  2 EA

On peut maintenant calculer le déplacement du nœud j :

Nij p (x2  x1)² p (x2  x1)(L  x1) Nij p (2L  x2  x1)(x2  x1)
uj = u(L) = ui + EA L  2 EA  = u i+ L 
EA EA 2 EA
On rappelle que Nij est indépendant de vi, vj, i et j . Exprimons maintenant Nij en fonction de ui
et uj :
EA p (2L  x2  x1)(x2  x1)
Nij = (uj  ui) L  2L
p (2L  x2  x1)(x2  x1) p (x2²  x1²)
Et notons Nij0 = puis Nji0 =
2L 2L
EA p (x2²  x1²)
car Nji = Nij  p (x2  x1) = (uj  ui) 
L 2L

Diagramme d'effort normal pour un chargement entre x1 et x2. Trois possibilités sont représentées
chacune pour une valeur possible de Nij.

 Nij   EA/L EA/L   ui   Nij 

0
On peut mettre ces équations sous forme matricielle :  N  =  .  u  +  N 0
ji  EA/L EA/L  j ji

Enfin notons que Nij  Nji = p (x2  x1) = Nij0  Nji0 ce qui traduit que la différence d'effort normal
entre les extrémités est due au chargement extérieur et que celui-ci est reporté aux nœuds à l'aide
des termes Nij0 et Nji0.
Pour les systèmes dits triangulés, i. e. constitué de barres rotulées aux deux extrémités ou
assimilées comme telles, ces dernières équations suffisent pour résoudre le problème.

8.2.2. Exemple pour un système triangulé

Considérons deux barres assemblées par des rotules constituant la structure de la figure
suivante. On considérera que le produit de l'aire par le module d'élasticité, i. e. EA, est constant et
le même pour les deux barres.

46
 Structure triangulée à deux barres.

Pour cette structure isostatique, les trois équations d'équilibre peuvent suffire pour tout
déterminer. Rappelons que l’équilibre d’une barre – droite ou non - avec une rotule à chacune de
ses extrémités impose que la force transmise soit colinéaire à la direction définie par les deux
rotules. Nommons F1X et F1Y les composantes de réactions d'appui selon X et Y au nœud 1 et de
façon analogue F2X et F2Y pour le nœud 2. L'équilibre statique impose F1X + F2X + F = 0 et
F1Y + F2Y = 0 et F L + F2Y L = 0 d'où F2Y = F et F1Y = F et F2X = 0 et F1X = F. Il est aisé d'en
déduire ensuite les efforts dans les barres puis les déformations des barres.
Vérifions ceci par la méthode des déplacements. Les conditions aux limites imposent que
seul le déplacement du nœud 3 est non nul. Les 6 inconnues du problème sont F1X, F1Y, F2X, F2Y,
U3 et V3. Ecrivons les 6 équations d'équilibre des actions sur les nœuds :
1 1
Nœud 1 : F1X + N13 = 0 et F1Y + N13 = 0
2 2
Nœud 2 : F2X = 0 et F2Y + N23 = 0
1 1
Nœud 3 : F  N31 = 0 et  N31  N32 = 0
2 2
Il est commode d'utiliser les expressions suivantes pour les changements de repère. On rappelle les
caractères gras désignent des vecteurs.
2
 = (Xj – Xi)² + (Yj – Yi)²
Xj – Xi Yj – Yi Yj – Yi Xj – Xi
x= X+ Y y= X+ Y
   
Xj – Xi Yj – Yi Yj – Yi Xj – Xi
X= x y Y= x+ y
   
Il faut maintenant faire apparaître les déplacements. Les déplacements dans le repère local de la
1 1
barre 1-3 s'expriment à l'aide de leurs correspondants dans le repère global : u3 = U3 + V3 et
2 2
1 1
v3 =  U3 + V3. Etant donné que U1 = 0, V1 = 0, U2 = 0 et V2 = 0, il vient :
2 2
EA EA
Nœud 1 : F1X + (U3 + V3) = 0 et F1Y + (U3 + V3) = 0
2 2L 2 2L
EA
Nœud 2 : F2X = 0 et F2Y + L V3 = 0

47
 EA EA EA
Nœud 3 : F  (U3 + V3) = 0 et  L V3  (U3 + V3) = 0
2 2L 2 2L
Il suffit maintenant de résoudre ce système de 6 équations à 6 inconnues. On trouve finalement
FL FL
V3 =  EA, U3  EA (1 + 2 2), F1X = F, F1Y = F, F2X = 0 et F2Y = F.

8.2.3. Composantes liées à la flexion

Si l'on considère les déplacements selon l'axe y vi et vj ainsi que les rotations i et j, il
apparaît que tous déplacements relatifs des nœuds (vi  vj  0 et i = j, par exemple) induit de la
flexion dans la poutre. On ne peut donc pas dissocier ces quatre degrés de liberté du point de vue
des équations de flexion. Pour plus de simplicité, nous négligerons les déformations d'effort
tranchant devant celles du moment fléchissant, i. e. V/GSr <<M/EI. Comme précédemment, nous
considérons une charge q orientée selon y et uniformément répartie entre l'abscisse x1 et l'abscisse
x2, i. e. Q(s) = q si s  [x1, x2] et Q(s) = 0 si s  [x1, x2]. L'équilibre global de la poutre est traduit
par :
L. x2.
Vji  Vij + 
Q(s) ds = 0  Vji  Vij + q 
ds = 0  Vji  Vij = q (x1 – x2)
0 x1
L. x2.
s q(s) ds + Vji L = 0  Mji  Mij + q 
Mji  Mij +  s ds + Vji L = 0
0 x1
 Mji  Mij = q (x1 – x2)(x1 + x2)/2  Vji L
respectivement pour les forces projetées sur y et pour les moments selon Z calculés ici par rapport
au nœud i. L'effort tranchant varie donc selon :
0  x  x1 x1  x  x2 x2  L
V(x) = Vij V(x) = Vij  q (x  x1) V(x) = Vij  q (x2  x1) = Vji
Et le moment fléchissant selon :
0  x  x1 x1  x  x2 x2  L
M(x) = Mij  Vij x x. x2.
M(x) = Mij  Vij x + (x  s) Q(s) ds M(x) = Mij  Vij x + q 
(x  s) ds
0 x1
x. = Mij  Vij x
(x  s) ds
= Mij  Vij x + q  + q [(x – x1)²  (x – x2)²] / 2
x1
= Mij  Vij x + q (x – x1)² / 2

Diagrammes d'effort tranchant (g) et de moment fléchissant(d) pour un chargement entre x1 et x2,
Vij et Mij étant donnés.

On peut maintenant calculer la déformée en intégrant les rotations depuis le nœud i. On fait
le calcul pour le cas où EI est constant entre les nœuds i et j.
0  x  x1 x1  x  x2 x2  L
EI y1"(x) = Mij  Vij x EI y2"(x) = Mij  Vij x EI y3"(x) = Mij  Vij x
+ q (x – x1)2 / 2 + q [(x – x1)2  (x – x2)2] / 2

48
EI y1'(x) = Mij x  Vij x2/2 + C1
EI y2'(x) = Mij x  Vij x2/2 + C3 EI y3'(x) = Mij x  Vij x2/2 + C5
EI y1(x) = Mij x /2  Vij x /6
2 3
+ q (x – x1)3 / 6 + q [(x – x1)3  (x – x2)3] / 6
+ C1 x + C2
EI y2(x) = Mij x2/2  Vij x3/6 EI y3(x) = Mij x2/2  Vij x3/6
+ C3 x + C4 + C5 x + C 6
+ q (x – x1)4 / 24 + q [(x – x1)4  (x – x2)4] / 24

Les conditions de raccordement – continuité des déplacements et rotations - des tronçons 1, 2 et 3

se formalisent par y'1(x1) = y'2(x1), y'2(x2) = y'3(x2), y1(x1) = y2(x1), y2(x2) = y3(x2). Ceci conduit à
C1 = C3 = C5 et C2 = C4 = C6. Remarquons que EI y1'(0) = C1 donc C1/EI = i. De même
EI y1(0) = C2 donc C2/EI = vi. Faisons maintenant apparaître les rotations et déplacements
différentiels :
1 L2 q
j  i = y'3(L)  y'1(0) = {Mij L  Vij + [(L – x1)3  (L – x2)3]}
EI 2 6
1 L2 L3 q
vj  vi = y3(L)  y1(0) = EI {Mij 2  Vij 6 + 24 [(L – x1)4  (L – x2)4]} + i L

où le terme i joue le rôle d'une rotation solide dans cette dernière équation. En effet, dans le cas
d'une rotation solide, i. e. j  i et M(x) = 0, V(x) = 0 et q(x) = 0, la différence de déplacement des
nœuds i et j vaut vj  vi = tan(i) L. Or tan(i)  i puisque dans les hypothèses des petites
perturbations i << 1, d'où vj  vi = i L . Il faut ensuite éliminer le terme Vij dans les équations
précédentes pour pouvoir exprimer Mij à l'aide des degrés de liberté en flexion i, j,vi et vj. Il
vient :
L L2 q qL
EI (vj  vi)  3 (j  i)} = Mij 6 +24 [(L – x1)4  (L – x2)4]}  18 [(L – x1)3  (L – x2)3]}+ i L

–6 EI 6 EI 4 EI 2 EI
puis Mij = L2 vi + L2 vj  L i  L j + Mij0

q
où Mij0 = 12 L2 [(L + 3x1)(L – x1)3  (L + 3x2)(L – x2)3]

En éliminant de même Mij, il vient :

–12 EI 12 EI 6 EI 6 EI
L3 vi + L3 vj – L2 i – L2 j + Vij
0
Vij =

q
où Vij0 = 2 L3 [(L + x1)(L – x1)3  (L + x2)(L – x2)3]

En utilisant les deux équations d'équilibre global de la poutre établie en début de chapitre on
trouve en complément :
6 EI 6 EI 2 EI 4 EI
Mji = L2 vi – L2 vj  L i  L j + Mji0

q
où Mji0 = 12 L2 [(4L – 3x2) x23  (4L – 3x1) x13]

–12 EI 12 EI 6 EI 6 EI
L3 vi + L3 vj – L2 i – L2 j + Vji
0
Vji =

q
où Vji0 = 2 L3 [(2L – x1) x13  (2L – x2) x23].

49
On peut mettre sous forme matricielle ces équations pour obtenir :

 EI/L EI/L2 EI/L3 EI/L2

3
MVij
  vi Vij0
   Mij0 
 EI/L2 EI/L EI/L2 EI/L  i
 Vji   EI/L3 EI/L2   vj   Vji 
ij
= . +
EI/L2 EI/L3
0

 Mji   EI/L2 EI/L EI/L2 EI/L   j   Mji 

0

8.2.4. Matrice de rigidité d'une poutre

Les signes des termes de la matrice de rigidité suivante peuvent bien sûr changer selon la
convention de signes adoptée pour les efforts intérieurs mais aussi avec le choix de travailler avec
les efforts intérieurs aux extrémités de la poutre ou les forces et couples appliqués aux extrémités
de la poutre ou encore aux nœuds. Avec nos choix et conventions, on obtient :
EA/L
 
Nij0
 
0 0 +EA/L 0 0
   
Nij ui
0 EI/L EI/L
3 2
0 EI/L EI/L
3 2
vi Vij0
       
Vij
Mij 0 EI/L2 EI/L 0 EI/L2 EI/L i Mij0
= . +
 Nji
  EA/L 0 0 +EA/L 0 0
 
uj
  Nji0

 Vji
Mji  0
0
EI/L3 EI/L2
+EI/L2 +EI/L
0
0
EI/L3 EI/L2
EI/L2 +EI/L  
vj
j   Vji0
Mji0 
8.2.5. Exemple pour un système fléchi
Considérons une poutre fléchie par exemple sous l'action du vent en toiture représentée sur
la figure suivante. On considérera que le produit du moment quadratique par le module d'élasticité,
i. e. EI, est constant le long de la poutre. Il s'agit d'un cas où il n'y aura pas de changement de
repère à faire.

Poutre fléchie par exemple sous l'action du vent en toiture.

Les conditions aux limites sont telles que U1 = 0, V1 = 0, V2 = 0, M1 = 0 et M2 = 0. Les

inconnues non triviales de ce problème sont les réactions d'appui, F1X, F1Y, F2y, le déplacement U2
et les rotations 1 et 2. Il faut d'abord calculer les charges reportées sur les nœuds comme si ceux-
ci étaient bloqués en déplacement et rotation. Cela donne : N120 = 0, N210 = , V120 = qL / 2,
V210 = qL / 2, M120 = qL² / 12, et M210 = qL² / 12. L'équilibre des nœuds impose que :
Nœud 1 : F1X + N12 = 0 , F1Y + V12 = 0 et M12 = 0
Nœud 2 : N21 = 0 , F2Y  V21 = 0 et M21 = 0
Faisons apparaître les déplacements :
EA
Nœud 1 : F1X + L U2 = 0

6 EI 6 EI
F1Y  L2 1  L2 2 + qL / 2 = 0

4 EI 2 EI
 L 1  L 2 + qL2 / 12 = 0

50
 EA
Nœud 2 : L U2 = 0
6 EI 6 EI
F2Y  L2 1 + L2 2 + qL / 2 = 0

2 EI 4 EI
L 1 + L 2 + qL / 12 = 0
2

Il vient U2 = 0, F1X = 0 puis 1  2 = 0, 1  2 = qL3/12EI donc 1 = qL3/24EI et 2 = qL3/24EI

et enfin F1Y = qL/2 et F2Y = qL/2.
Conformément aux équations précédentes, on peut ensuite calculer, par exemple, le
moment fléchissant entre les nœuds. M(x) = Mij  Vij x + q (x – x1)² / 2, i. e.
M(x) = qL/2 x + q x² / 2. Pour calculer y(x) et y'(x), on se reporte également aux équations
précédentes dans lesquelles C1 = i/EI et C2 = vi/EI sont désormais connus.

8.2.6. Formule des trois moments à partir de la méthode de déplacements

Poutre continue sur appuis.

A partir de la matrice de rigidité de la poutre on peut écrire, en supposant EI constant :

4 EI 2 EI
MA =  A   + MAB0
LAB LAB B
4 EI 2 EI
MB =  L B  L A + MBA0
AB AB

4 EI 2 EI
MB =  L B  L A + MBC0
BC BC

4 EI 2 EI
MC =  L C  L B + MCB0
BC BC

où LAB et LBC sont respectivement les longueurs des poutres A-B et B-C. En éliminant A et C on
obtient :
6 EI
MA + 2 MB = L B + MAB0 + 2 MBA0
AB
6 EI
MC + 2 MB =  L B + MCB0 + 2 MBC0
BC
En éliminant B il vient, à condition que EI soit le même dans les poutres de gauche et droite :

LAB MA + 2 (LAB + LBC) MB + LBC MC = LAB (MAB0 + 2 MBA0) + LBC (MCB0 + 2 MBC0)

On peut déjà utiliser cette formule pour résoudre relativement facilement certains
problèmes de poutres continues sur appui si l’on ne dispose que du formulaire de la méthode des
déplacements. Pour faire apparaître la formule des trois moments classique, on doit considérer que
l'on connaît les rotations sur appuis des poutres A-B et B-C celles-ci étant séparées en condition
isostatique. Ces rotations sont désignées par d iso et g iso, d et g signifiant que l'on est à droite ou à

51
gauche de l'appui. Puisque les déplacements sur appui vA, vB et vC sont annulés (par construction),
remarquons que la matrice de rigidité nous permet de calculer aisément ces rotations qui peuvent
se mettre sous la forme :
g iso = LAB (MAB0 + 2 MBA0)
d iso = LBC (MCB0 + 2 MBC0)
On obtient finalement la formule des trois moments ou formule de Clapeyron :
LAB MA + 2 (LAB + LBC) MB + LBC MC = d iso  g iso)
Dans le cas d'une charge uniformément répartie q portée par Y, on utilise d iso = qLBC3/24EI et
g iso = qLAB3/24EI.

8.2.7. Exemple de mise en œuvre de la méthode dite des 3 moments

Poutre continue sur 3 appuis uniformément chargée.

Supposons une poutre continue sur 3 appuis uniformément chargée par une charge q portée par
Y. Aux extrémités non chargées par des couples extérieurs A et C on peut affirmer que
MA = MC = 0. A l’aide des équations précédentes, on peut écrire que :
2 (LAB + LBC) MB = qLBC3/4  qLAB3/4)
d’où MB = q(LBC3  LAB3)/[8(LAB + LBC)]. Si LAB = LBC = L, on retrouve bien que MB = qL2/8.
Ensuite, sachant que le moment fléchissant en B, désigné par MB, peut être calculé à partir des
efforts à gauche de B, il vient MB = FAY L  qL2/2 on déduit que la réaction d’appui selon Y en A
vaut FAY = qL/2  qL/8 = 3qL/8. De même, en calculant le moment fléchissant en B à partir des
efforts à droite de B il vient FcY = 3qL/8. Enfin, la réaction d’appui selon Y en B doit assurer
l’équilibre statique des efforts selon Y de sorte que FAY + FBY + FCY  2qL = 0 donc FBY = 5qL/4.

8.3. Formulaire

Action des nœuds sur la poutre.

Action de la poutre sur les nœuds.

52
 Vij0 = q [(L + x1)(L  x1)3  (L + x2)(L  x2)3] / (2L3)
Vji0 = q [(x13 (2L  x1)  x23 (2L  x2)] / (2L3)
Mij0 = q [(L + 3x1)(L  x1)3  (L + 3x2)(L  x2)3] / (12L2)
Mji0 = q [x23 (4L  3x2)  x13 (4L  3x1)] / (12L2)
Vij0 = qL / 2
Vji0 = qL / 2
Mij0 = qL2 / 12
Mji0 = qL2 / 12

Vij0 = F [(L + 2a) (L  a)2] / (L3)

Vji0 = F [(2a  3L) a2] / (L3)
Mij0 = F [(L  a)2 a] / (L2)
Mji0 = F [(L  a) a2] / (L2)

Nij0 = p (x2  x1)(2L  x2  x1) / (2L)

Nji0 = p (x2²  x1²) / (2L)

Si x1  x2  a et F = p(x2  x1)
Nij0 = F (L  a) / L
Nji0 = F a / L
Report de charge sur les nœuds provenant de chargement entre les nœuds.

9. Flambement et instabilités de forme

9.1. Hypothèses de calcul
Les H. P. P. ne sont plus totalement vérifiées car les déplacements sont grands, les
déformations restent petites et en élasticité linéaire. Il faut travailler dans la configuration
déformée pour connaître les sollicitations extérieures. Comme il faut connaître les sollicitations
extérieures pour calculer la déformée et réciproquement, on aboutit inévitablement à une équation
différentielle.

9.2. Cas de la compression simple – flambement d'Euler

Considérons un poteau simplement comprimé tel que E(x1) IG3(x1) et N(x1) sont constants.
y"(x1) = M3G(x1) / E IG3
M3G(x1) = N y(x1) N < 0 en compression
D'où y"(x1) = N y(x1) / E IG3
La famille des solutions de ce problème est de la forme y(x 1) = a cos  x1 + b sin  x1.
Considérons un poteau rotulé en tête et pied de longueur L, i. e. y(0) = 0 et y(L) = 0. On en déduit
que b = 0 puis a cos  L = 0 donc que  = n  / L (n entier), enfin que ² = N / E IG3. Il vient
donc n² ² / L² = N / E IG3. Concernant la constante a, elle peut prendre une valeur quelconque ce
qui indique une situation instable, la déformée pouvant s'amplifier à effort normal constant. Par
ailleurs, la valeur minimale de N sera obtenue pour n = 1. On en déduit qu'il ne faut pas dépasser
une force de compression critique évaluée à :
Nc = ² E IG3 / Lf²
où Lf est appelée longueur de flambement et vaut L dans le cas de deux liaisons avec rotules. La
longueur de flambement prend d'autres valeurs pour d'autres conditions de liaison : L/2 pour deux
2
encastrements et L/ 2 pour un côté avec rotule et l’autre encastré.

53
 espace dessin et calculs, longueurs de flambement avec encastrement(s)

Pratiquement, le flambement - si l'on n'a pas réussi à l'éviter – en compression simple est
déclenché par des imperfections géométriques et/ou matériau puis amplifié par la déformation. Les
sections finissent par tourner significativement et certaines des équations utilisées pour démontrer
l'existence de l'instabilité deviennent approximatives. Le flambement – ou plutôt non-flambement -
d'un élément de structure peut être étudié relativement simplement mais le flambement d'ensemble
d'une structure nécessite un ordinateur. Les mots clés pour une éventuelle recherche sont :
flambage, flambement, cloquage, voilement d’âme, déversement et instabilité de forme.
Lorsque que la sollicitation combine flexion et compression, il convient d’étudier avec soin
l’instabilité de forme moment fléchissant et effort normal admissibles sont fortement liés.
Heureusement, les règlements de construction prévoient que de nombreux cas pour lesquels
l’élancement est relativement faible puissent être examinés de façon forfaitaire.

10. Exemple de résolution d’un problème de portique

10.1. Résolution par la méthode des forces
Considérons le portique de la figure suivante de gauche. Les poteaux sont encastrés en A et
E, leur rigidité en flexion est définie par EIA. Les poutres sont encastrées sur les poteaux et
rotulées entre elles, leur rigidité en flexion est définie par EIC. Une charge q de type poids propre
est considérée en toiture. Il s’agit de montrer que si l’on modifie les rigidités relatives des poutres
et poteaux, i. e. le rapport IA / IC, on change des valeurs de réaction d’appui, contrairement à une
structure isostatique.
Le calcul du degré d’hyperstaticité donne d = 3 + 3 + 2 + 3 + 3  3x4 = 2. En effet si, par
exemple, on supprime les blocages en rotation des pieds de poteaux, i. e. on remplace les
encastrements par des appuis simples, alors on obtient un portique à 3 rotules que nous
connaissons isostatique.
La symétrie de la structure et du chargement permet d’affirmer que le déplacement du
point C ne peut être que selon Y. La cinématique de la figure suivante de droite pour un demi-
portique permet donc de résoudre le problème plus simplement en réduisant le nombre
d’inconnues de 2 à 1. Rappelons que le degré d’hyperstaticité est indépendant du chargement et,
même avec un chargement symétrique, la structure demeure hyperstatique de degré deux. Du point
de vue des équations, poser par exemple Uc = 0 par raison de symétrie comme équation
supplémentaire permet de ne plus avoir à trouver qu’une seule équation manquante qui, s’ajoutant
aux trois équations d’équilibre permettra de résoudre le problème11.

11
. On aurait aussi pu affirmer que les réactions d’appui verticales valent chacune 5qL par raison de symétrie
et d’équilibre.

54
Structure symétrique symétriquement chargée. Simplification autorisée par la symétrie.
Pour résoudre, on peut par exemple choisir la structure isostatique des figures suivantes,
l’ « inconnue hyperstatique » étant la réaction d’appui glissant au point C. Il s’agit d’exprimer que
la superposition des efforts du système 0 (S0) et de FXC fois les efforts système 1 (S1) fournit le
système à résoudre. Pour cela, la superposition des déformations doit respecter les conditions du
système à résoudre. Donc, l’addition du déplacement de C selon X dans S0, noté uc0, et du
déplacement de C selon X dans FXC.S1, noté uc1, doit donner 0.
FXC doit être tel que uc0 + uc1 = 0
Rappelons que l’utilisation de la superposition nécessite de respecter les hypothèses des petits déplacements et du
comportement élastique linéaire de la matière.

Système 0 (S0) structure isostatique associée Système 1 (S1), structure isostatique associée
chargée par les données de la descente de chargée par 1 N à multiplier éventuellement
charge. par FXC.
Calcul de uC0 :
On fait travailler les efforts de S1 dans les déformations de S0. Pour plus de lisibilité, la
charge unitaire virtuelle est notée 1*. En vertu du principe des travaux virtuels, en négligeant les
déformations d’effort normal et tranchant devant celles du moment fléchissant, on écrit donc que

1* uC0 =  m1*(x) m0(x) / EI(x) dx =


toute la Structure


 m1*(x) m0(x) / EIA dx +  m1*(x) m0(x) / EIC dx
poteau poutre

55
Le calcul de l’intégrale donne 1* uC0 = H/EIA (-10qL²) 1/2 ({-3L-H-3L}1*) + 5L/EIC 1/3 (-10qL²)
3/4 (-3L1*) d’où l’on déduit que uC0 = 5qL2H(6L+H)/EIA + 75qL4/2EIC porté par X.

Calcul de uC1 :
On fait travailler les efforts de S1 dans les déformations de FXC S0. Comme précédemment,
x

on écrit donc que

1* uC1 =  C
 m1*(x) FX m1(x) / EI(x) dx
x

toute la Structure
Le calcul de l’intégrale donne 1* uC1 = H/3EIA 1/6 2(9L²+{3L+H}²+3L{3L+H}) FXC 1* + 5L/EIC
1/3 (-3FXCL) (-3L1*) d’où l’on déduit que uC1 = H/3EIA ({27L²+H²+9LH}FXC) + 15L3/EIC FXC
porté par X. FXC doit être tel que uc0 + uc1 = 0 donc
5L2H(6L+H)/IA + 75L4/2IC
FXC = -q H(27L²+H²+9LH)/3I + 15L3/I
A C

Si l’on désigne par  le rapport IC/IA on peut faire mieux apparaître la dépendance de cette
réaction d’appui à rigidité relatives poutres/poteaux.
5L2H(6L+H) + 75L4/2
FXC = -q
H(27L²+H²+9LH)/3 + 15L3/
Considérant toujours le demi-système hyperstatique, l’équilibre impose que la composante selon X
de la réaction d’appui en A vaut -FXC et le moment d’encastrement vaut 10qL² + (3L+H) FXC porté
par Z. Dans ce cas, seule la composante selon Y en A qui vaut 5qL est indépendante de IC/IA.

Equilibre et diagrammes du système 0. Equilibre et diagrammes du système 1.

Convention avec T = -V pour l'effort tranchant. Convention avec T = -V pour l'effort tranchant.

56
 10.2. Remarque concernant les déplacements des têtes de poteau dans des
portiques
En général les déformations liées à l’effort normal sont négligées devant celles de flexion,
liées au moment fléchissant. Bien sûr, dans les systèmes triangulés (constitués de « bielles ») pour
lesquels seul l’effort normal est non nul, il convient de considérer les déformations de cet effort
pour calculer les déformations. Si l’on néglige les déformations d’effort normal, une tête de
poteau peut-elle monter ou descendre sous l’effet de la rotation du poteau ? Considérons un cas
simple où le poteau est rotulé en pied (si le poteau est encastré, la rotation globale sera moindre).

Rotation du poteau et cinématique en négligeant les déformations d’effort normal.

Pour la figure précédente cela revient à calculer le terme v sachant que la tête de poteau se
déplace sur un cercle de rayon L et de centre A. v est la projection sur la verticale de la différence
entre la longueur de l’hypoténuse du triangle rectangle (ABB’) et la longueur L. v est donc
inférieur à cette différence. Les flèches sont typiquement limitées réglementairement à L/200 et la
structure devra être dimensionnée pour respecter, entre autres, ce critère. Il vient donc :
V ≤ √(L2 + (L/200)2) - L
Ce qui donne v ≤ 0.0000125 L. La cinématique de la déformation permet donc d’assimiler l’état
initial et l’état déformé. Négliger les déformations d’effort normal revient donc dans ce type de cas
à considérer que le déplacement vertical de B est nul.
Rappelons tout de même que l’étude des instabilités de flambement ne permet pas de faire
l’assimilation entre la géométrie à l’état déformé et celle de l’état initial …

57
Annexe 1 : Un peu de géométrie pour un tétraèdre
On définit un triangle passant par les points a, b et c distincts. On se place hors des cas
particuliers où un de ces points correspond à 0 qui ont des solutions évidentes. Les points sont
définis par Oa = (a1, 0, 0), Ob = (0, b2, 0), Oc = (0, 0, c3). La question posée est : quelle est la
surface s de ce triangle ?

espace dessin, tétraèdre défini par a, b et c

__
Pour cela, il faut situer le point d sur la droite ab de sorte que a, c et d ainsi que b, c et d
définissent deux triangles rectangles. On peut donc calculer s ainsi :
s = ½ (║ad║½ ║cd║½ + ║db║½ ║cd║½)
Etant donne que ad est colinéaire à db on peut réécrire :
s = ½ ║ab║½ ║cd║½
__
d  ab définie par x1/a1 + x2/b2 = 1. De plus cd doit être perpendiculaire à ab. Puisque
cd = cO + Od = (d1, d2, c3) et ab = aO + Ob = (a1, b2, ), on écrit donc :
(d1, d2, c3) . (a1, b2, ) = 0  d1 = d2 b2 / a1
en reportant dans l'équation de la droite on obtient d2 b2 / a1² + d2 / b2 = 1, d'où :
d1 = a1 b2² / (a1² + b2²) et d2 = a1² b2 / (a1² + b2²)
On peut maintenant calculer ║cd║ :
║cd║2 = (a14 b22 + a12 b24) / (a12 + b22)2 + c32 = (a12 b22) / (a12 + b22) + c32
2
Puis enfin s=½ a12 b22 + a12 c32 + c32 b22
Ensuite on cherche la normale unitaire On = (n1, n2, n3) à cette surface. Le vecteur On devra être
perpendiculaire à ab, ac et par conséquence bc. On peut donc écrire que a1 n1 + b2 n2 = 0 et que
a1 n1 + c3 n3 = 0 d'où n2 = n1 a1 / b2 et n3 = n1 a1 / c3. Pour que On soit unitaire il faut que
n12 + n22 + n32 = 1 d'où :
n12 (1 + a12 / b22 + a12 / c32) = 1  n12 (a12 b22 + a12 c32 + c32 b22) = c32 b22  n12 = b22 c32 / 4 s2
Et on obtient finalement n1 = b2 c3 / 2 s et de même n2 = a1 c3 / 2 s et n3 = a1 b2 / 2 s. Donc :
2 s n1 = b2 c3 2 s n2 = a1 c3 2 s n3 = a1 b2

58
Annexe 2 : Section réduite liée au cisaillement
La loi de Hooke de l’élasticité linéaire exprime que la contrainte de cisaillement  est
proportionnelle à la distorsion  :

=µ
ou encore pour la notation de la MMC ij = µ ij = 2 µ ij = 2 ij car ij  2 ij

L’équilibre (ceci revient à écrire que div  = 0 pour calculer l'évolution de(x2) /x2 en fonction
de (x1) /x1, cette dernière évolution étant liée à l'évolution du moment fléchissant selon x1 donc
à l'effort tranchant) de la section impose qu'en x1 :

T S(x2) T S(x2)
(x2) = I b(x ) donc  = µ I b(x )
G3 2 G3 2

vsup.
On rappelle que S(x2) = 
b()  d
x2
1
La densité d’énergie élastique de cisaillement vaut 2   Un élément de volume pour un tronçon de
poutre de longueur dl vaut d = b(x2) dx2 dl. En intégrant sur la section on obtient :

vsup.
= 2(x2) (x2) b(x2) dx2
1
dW dl

vinf.

en remplaçant  et  par leurs expressions on fonction de V, S(x2), µ, IG3 et b(x2) on obtient :

vsup.

= 
1 T² S²(x2)
dW  2µ I ² b(x ) dx2 dl
 G3 2
vinf.

vsup.
dW 1 T² S²(x2)
 dl = 2µ IG3²  dx
 b(x2) 2
vinf.
où apparaît un terme caractérisant la géométrie de la section. On introduit le terme "section
vsup.
T 1 1 S²(x2)
réduite" de sorte que équivalent = µ S d’où :
Sr IG3² 
= dx2
r  b(x2)
vinf.
Par exemple pour une section rectangulaire de dimensions b x h, on obtient :
h/2
1 144 6 5bh 5A
= 2 6 
 {b [(h/2)²  x 2²] / 2}² / b dx2 =  Sr = 6 = 6
Sr b h 5bh
h/2
A l’aide de ces notations, on peut donc réécrire :
dW T²
=
dl 2 µ Sr

59
A3. Démonstration inverse du TFS à partir du PTV…
…ou du théorème fondamental de la dynamique (T.F.D.) à partir du principe des
puissances virtuelles (P.P.V.). Si l’on avait postulé le principe des travaux virtuels (P.T.V.), on
aurait procédé à la démonstration du théorème fondamental de la statique. Pour ce faire, il suffit de
considérer un mouvement de solide rigide quelconque tel que le travail des efforts intérieurs est
nul. Le travail des efforts extérieurs doit donc être nul lui aussi. On démontre donc par translations
que la somme des forces est nulle et par rotations que la somme des moments par rapport à tout
point est nulle. Il est plus facile de faire cette démonstration inverse mais généralement l'équilibre
est plus intuitif que les puissances virtuelles.

Annexe 4 : En quoi la théorie des poutres n'est pas une solution

exacte ?
A4. 1 Gauchissement des sections
Etudions plus en détail l'admissibilité de la solution proposée par la théorie des poutres. On
peut dire qu'elle est partout cinématiquement admissible et qu'elle respecte partout la loi de
comportement. Est-elle cinématiquement admissible. En annexe 2, pour connaître la distribution
V S(x2)
de la contrainte de cisaillement, nous avons établi que (x2) = en écrivant que l'évolution
IG3 b(x2)
de (x2) satisfait à l'équilibre. Ceci revient à écrire que div  = 0. Ainsi, il est montré que la
contrainte de cisaillement est nulle en fibre supérieure (x2 = v) et en fibre inférieure (x2 = -v'). Cela
est-il compatible avec la cinématique proposée, i. e."les sections droites restent droites". Comme le
montre la figure suivante, il est clair que, les faces supérieures et inférieures de la poutre étant
libres de contrainte de cisaillement, du fait de la réciprocité des contraintes tangentielles, on ne
peut pas équilibrer une contrainte de cisaillement en fibre supérieure et inférieure avec l'hypothèse
des sections qui restent droites. Les sections se gauchissent donc obligatoirement en présence
d'effort tranchant. Ce gauchissement reste négligeable dans la plupart des cas.

Incompatibilité des conditions aux limites en cisaillement.

On peut tout de même remarquer que seule une évolution du gauchissement compliquerait
le problème d'évaluation des contraintes normales. En effet, si le gauchissement est constant, il est
toujours aisé de calculer les déformations (allongement ou raccourcissement) des fibres
longitudinales dans la section et la distribution des contraintes normales est linéaire malgré le
gauchissement si celui-ci est constant.

60
 Gauchissement des sections sans variation du gauchissement.

A4. 2 Flexion circulaire

Supposons que seul le moment fléchissant soit non nul, i. e. "flexion circulaire", comme
par exemple en section centrale d'une poutre en flexion "4 points". Comme le montre la figure
suivante, la cinématique ne conserve pas les angles droits des coins. Par exemple, pour le coin
supérieur droit, la distorsion vaut /2. Nous avons établi que x1 = y" = M/EI. Pour un
segment de longueur dx1 nous avons donc d = x1 dx1 = M dx1/EI. La distorsion vaut donc
 = M dx1/2EI qui est un terme du second ordre. La flexion seule n'engendre donc pas de
gauchissement considérable. Dans le cas de flexion circulaire, la théorie des poutres fournit donc
un outil précis.

La flexion ne génère pas de distorsion au premier ordre.

A4. 3 Loin des points d'application des charges…

C'est ici que la théorie des poutres connaît son point faible. Elle ne prédit pas de contraintes
normale selon x2 et ne peut donc pas assurer l'équilibre. Toutefois, le calcul de la déformée
intégrant les déformations partout, ces erreurs locales sont difficilement mesurables lorsqu'on
s'intéresse aux déplacements (intégration de y"). En outre, si l'on se reporte à l'exemple du calcul
de déformée, on constate que la contrainte 22 est généralement petite devant la contrainte 11
normale prédite par la théorie des poutres. Néanmoins, il faut procéder à des vérifications
particulières dans les zones de charge concentrée.

Charge concentrée et contrainte normale 22.

61
Annexe 5 : Théorème des travaux virtuels, que serait la
démonstration pour une poutre droite ?
A5. 1 Préambule
Cette partie n'est là que pour montrer que cette démonstration n'est pas plus simple avec la
théorie des poutres - bien au contraire - mais tout de même possible. Il faut donc la considérer
comme un exercice de style. Nonobstant cette démonstration est bien peu élégante. Postuler le
P. P. V. est généralement la solution la plus simple pour les démonstrations.

A5. 2 Champs cinématiquement admissible et statiquement admissible

Considérons une déformation de poutre constituée d'une composante associée à l'effort
normal, N = u1 / x1, une à l'effort tranchant, V = u2 / x1, et une à la flexion, y"(x1) =  / x1.
Cette déformation est cinématiquement admissible, elle respecte les conditions d'appui, ne génère
pas de brisure ni d'interpénétration de matière. Considérons par ailleurs une sollicitation de cette
même poutre constituée de charges réparties (p(x1), q(x1), r(x1)) et de torseurs résultants
(N*, V*, M*) aux extrémités de sorte que cette poutre soit en équilibre.

Attention à la notation :V pour l'effort tranchant et v pour sa déformation.

espace dessin et calculs, déformation cinématiquement admissible, sollicitation statiquement

admissible.

62
Pour un tronçon de poutre tel que dx1 est petit, l'équilibre impose que :
N(x1) + N(x1 + dx1) + p.dx1 = N(x1) + N(x1) + dN(x1) + p.dx1 = 0
V(x1) + V(x1 + dx1) + q.dx1 = V(x1) + V(x1) + dV(x1) + q.dx1 = 0
M(x1) + M(x1 + dx1) + r.dx1 + V(x1).dx1  q(x1).dx1²/2 
M(x1) + M(x1) + dM(x1) + r.dx1 + V(x1).dx1 = 0
Donc :
dN(x1) = p.dx1  N(x1) /x1 = p
dV(x1) = q.dx1  V(x1) /x1 = q
dM(x1) = {V(x1) + r}.dx1  M(x1) /x1 = {V(x1) + r}

A5. 3 Travail des efforts intérieurs

Le travail des efforts intérieurs (N*, V*, M*) dans les déformations (N, V, y") vaut :
L
Wint. = 
{N*(x1).N(x1) + T*(x1).T(x1) + M*(x1).y"(x1)} dx1
0

A5. 4 Termes liés à l'effort normal

Calculons le terme Wint. N lié à l'effort normal :
L
Wint. N = 
N*(x1).N(x1) dx1
0
En remarquant que (N*u)’ = N*.N  N*'.u puisque u’ = N et en intégrant par partie, on obtient :
L L L
Wint. N = 
 [(N*u)’  N*'.u] dx 1 = N*.u│  
p*.u dx1
0 0 0
L
Wint. N = N*(L).u(L)  N*(0).u(0)  
p*.u dx1
0
Lorsque pF* représente une charge concentrée (ou plusieurs, il suffira de considérer une
superposition), c'est à dire si :
x1F  
pF*(x1) = 0 si x1  [x1F   x1F   et 
pF*(x1) dx1 = F1
*

x1F  
où x1F désigne le barycentre de la charge concentrée et F1* une force ponctuelle d'axe x1. Si  tend
vers 0 alors on peut considérer que u1(x1) est quasiment constant dans l'intervalle
[x1F   x1F  et vaut u1(x1F). On peut donc dire dans ce cas que :
L

pF*.u1 dx1 = F1*.u1(x1F)
0

63
La démonstration rigoureuse fait intervenir la fonction de Dirac (x1F) de sorte que l'on écrit
pF*(x1) = F1.(x1F), puis :
L L

F1.(x1F) u1(x1) dx1 = F1 
(x1F) u1(x1) dx1 ≡ F1*.u1(x1F)
0 0
Toutefois, il faut se rappeler que, dans le cadre de l'hypothèse de Saint Venant, la théorie des
poutres ne s'applique que loin des points d'application des charges ponctuelles. Pratiquement, la
modélisation mathématiquement rigoureuse n'a donc pas beaucoup plus de sens que
l'approximation d'une charge ponctuelle à une charge répartie sur une courte distance.

A5. 5 Termes liés à l'effort tranchant

Calculons de même le terme Wint. V lié à l'effort tranchant. Notons v la déformation d’effort
tranchant d’axe x2.
L
V*(x1).T (x1) dx1
Wint. V = 
0
En remarquant que (V*v)’ = V*.V  V*'.v puisque v’ = V et en intégrant par partie, on obtient :
L L L
[(V*v)’  V*'.v] dx1 = V*.v│  
Wint. V =  q*.v dx1
0 0 0
L
q*.v dx1
Wint. V = V*(L).v(L)  V*(0).v(0)  
0
Lorsque qF* représente une charge concentrée (ou plusieurs, il suffira de considérer une
superposition), c'est à dire si :
x1F  
qF*(x1) = 0 si x1  [x1F   x1F   et 
qF*(x1) dx1 = F2
*

x1F  
où x1F désigne le barycentre de la charge concentrée et F2* une force ponctuelle d'axe x2. Si  tend
vers 0 alors on peut considérer que v(x1) est quasiment constant dans l'intervalle
[x1F   x1F  et vaut v(x1F). On peut donc dire dans ce cas que :
L

qF*.v dx1 = F2*.v(x1F)
0
La démonstration rigoureuse fait intervenir la fonction de Dirac (x1F) de sorte que l'on écrit
qF*(x1) = F2.(x1F), puis :
L L

Fy.(x1F) v(x1) dx1 = F2 
(x1F) v(x1) dx1 ≡ F2*.v(x1F)
0 0

A5. 6 Termes liés au moment fléchissant

Calculons de même le terme Wint. M lié au moment fléchissant. Notons y la déformation de
flexion d’axe x2.

64
 L
Wint. M = 
M*(x1).y"(x1) dx1
0
En remarquant que (M*y')’ = M*.y"  M*'.y' et en intégrant par partie, on obtient :
L L L
Wint. M = 
[(M*y")’  M*'.y'] dx1 = M*.y'│  
(V* + r*).y' dx1
0 0 0
L L
Wint. M = M*(L).y'(L)  M*(0).y'(0)  
r*.y' dx1  
V*.y' dx1
0 0
Lorsque rF* représente un couple ponctuel (ou plusieurs, il suffira de considérer une
superposition), c'est à dire si rF*(x1) = C3.(x1F) où x1F désigne le point d'application du couple et
C3* le couple d'axe 3, alors il vient :
L L

C3.(x1F) y'(x1) dx1 = C3 
(x1F) y'(x1) dx1 ≡ C3*.y'(x1F)
0 0
où y'(x1F) est la rotation d'axe 3 du point x1F aussi notée (x1F). Concernant le terme d'effort
tranchant, on remarque que (V*.y)' = V*.y' + V*'.y d'où
L L L L

 V*.y' dx 1 = 
 [(V*.y)' + q*y] dx 1 = V*.y│  
q*.y dx1
0 0 0 0
L
= V*(L).y(L)  V*(0).y(0) + 
q*.y dx1
0
où l'on peut à nouveau considérer dans le cas d'une charge ponctuelle F2* (ou plusieurs charges
ponctuelles F2i* associées aux abscisses x1Fi, en superposant) que :
L

q*.y dx1 = F2*.y(x1F)
0

A5. 7 Total du travail des efforts intérieurs

Notons u2 le déplacement d’axe 2, i. e. u2 = y + v. Il vient :

Wint. = N*(L).u(L)  N*(0).u(0) + V*(L).v(L)  V*(0).v(0)

L L L
+ M*(L).(L)  M*(0).(0)   q*.u2 dx1 + 
p*.u1 dx1 +  r*. dx1
0 0 0
n m p
 F1i*.u1(x1Fi)  F2j*.u2(x1Fj)  C3k*.(x1Fk)
i=1 j=1 k=1
où l'on reconnaît dans le second terme de cette égalité le travail des efforts extérieurs. On peut
donc affirmer que, si les champs sont admissibles, le travail des efforts intérieurs est égal au travail
des efforts extérieurs. Ceci est aussi vrai dans le cas de poutres courbes, de plaques etc. mais, par
souci de simplicité, l'on se reportera à une démonstration plus générale (voir ci-avant) plutôt que

65
de compliquer la démonstration avec une courbure ou d'autres sources de complexité. On
remarque que l’on n’a pas fait ici d’hypothèses de petits déplacements ni d’élasticité linéaire.

Annexe 6 : Relation de conversion entre rayon de courbure et

contraintes internes résiduelles
A6. 1 Introduction
On montre ici comment la théorie des poutres peut servir en mécanique des matériaux.
Lors d'un refroidissement ou d'une trempe, la matière se fige de façon non uniforme et des
contraintes résiduelles apparaissent. Elles sont recherchées, par exemple dans le cas du verre, pour
mettre la "peau" en compression afin de rendre plus difficile un démarrage de fissuration en
surface. En outre, les vitesses de refroidissement diffèrent entre le cœur et la peau. La
microstructure du matériau peut dépendre de la vitesse de refroidissement et donc influer sur ses
propriétés mécaniques. Le cas des polymères semi-cristallins est un exemple de ce type de
situation. La cristallisation est influencée par la vitesse de refroidissement. A cœur, là où la vitesse
est plus lente, les sphérolites sont plus gros qu'en peau. La proportion de phase amorphe rapportée
à la phase cristalline peut aussi varier. Ainsi, les propriétés mécaniques peuvent varier dans
l'épaisseur, en même temps que les contraintes résiduelles.
On peut tenter de mesurer ces propriétés mécaniques locales par microindentation (Alchikh
et al., 2010; Alchikh, 2009). On peut aussi obtenir une information à l'aide de la théorie des
poutres.

A6. 2 Enoncé du problème

Les contraintes internes résiduelles dans une plaque de polymère, par exemple, peuvent
être mesurées par la méthode d’enlèvement de couches ou lamelles. La formule de Stoney (Stoney,
1909) permet de calculer ces contraintes résiduelles en mesurant le rayon de courbure de la plaque
après l’enlèvement successif de lamelles. Dans cette formule, le module d’Young est supposé
constant dans la hauteur de la plaque (axe y). Cependant dans notre développement original
suivant, le module d’Young peut varier en fonction de la hauteur y.
La plaque du polymère de hauteur initiale h0 est supposée multicouche et composée de
n lamelles, le module d’Young E  y  peut varier continûment dans la hauteur de la plaque mais
pas dans l'épaisseur (axe z). La contrainte résiduelle pour chaque fibre  R  y  est donc constante
selon z dans le plan de la fibre.

Figure A. 6. 1. : schéma de calcul : à gauche : la plaque est divisée en plusieurs lamelles. A droite
une distribution supposée de contraintes internes résiduelles, pour les fibres comprimées et
pour les fibres tendues.

Les déformations générées par les contraintes résiduelles dans les différentes fibres sont
supposées élastiques. D’après l’hypothèse de Navier-Bernoulli de la théorie des poutres « Les
sections planes et normales à la fibre moyenne avant la déformation restent planes et normales à la

66
fibre moyenne après la déformation » (Laroze, 1988), le déplacement d’une section droite est un
déplacement de corps solide. Les sections droites restent droites (ne gauchissent pas) et se
translatent et tournent simplement les unes par rapport aux autres. Nous écrivons l’équation du
comportement de déformation élastique initiale à faible amplitude pour chacune des fibres. La
« loi de Hooke » en théorie des poutres se résume à :
 R  y   E  y  0  y 

avec une déformation initiale d’une fibre 0  y   L l l0y y   L  lL0  y  , où L désigne la

0 
longueur de la plaque en absence de force extérieure et l0  y  la longueur de la fibre "libre de
contrainte", c'est à dire non déformée par les contraintes résiduelles.

A6. 3 Equations d’équilibre de la poutre en absence de force extérieure

L’effort normal N (égale à la somme des contraintes par rapport à la section) est nul :

 R  y  dy  ép.  0 0 E  y   0  y  dy  0
h0 h
N0  ép. 
0

où ép. désigne l'épaisseur des lamelles dans la direction z. La somme du moment fléchissant M
par rapport à la section est nulle aussi :

 R  y   y  ym  h0  dy  ép.  0 0 E  y   0  y   y  ym  h0  dy  0
h0 h
M 0  ép. 
0

où la fibre moyenne ym est définie comme :

0 y E  y  dy
h

 y  ym  E  y  dy  0  ym  h  
h
0
0 E  y  dy
h

Si maintenant nous enlevons une lamelle, la plaque retrouve une situation d'équilibre statique par
une déformation de flexion et d'allongement ou rétrécissement. Nous pouvons écrire alors :

  y   0  y   a y  b   inf. 0inf. h sup. 0 sup.

    

  y  désigne la déformation à y après l’enlèvement de la lamelle. L'indice 0 correspond à l'état

avant enlèvement de matière.
La courbure a est définie par :

 inf.  0inf.   sup.  0 sup.

a 1  1 
R R0 h
R0 et R : Rayon de courbure de la plaque mesuré avant et après l’enlèvement de lamelle.
0inf. et  inf. : Déformation de la fibre inférieure avant et après l’enlèvement de lamelle.
0 sup.et  sup. : Déformation de la fibre supérieure avant et après l’enlèvement de lamelle.
b : Déformation de la fibre inférieure (en y = 0) liée à l'enlèvement de lamelle.
Les équations d’équilibre de la plaque deviennent :
L’équation de l'effort normal par unité d'épaisseur :

67
 0 E  y  0  y   a y  b dy  0
h

0 E  y   0  y  dy  h E  y   0  y  dy  0 E  y   a y  b  a ym (h)  a ym (h)  dy  0
h0 h0 h

0  h 0  R  y  dy  0 E  y  a  y  ym (h)  dy  0 E  y   a ym (h)  b  dy  0
h h h

h  R  y  dy  0   a y
h0
 m ( h) b  E h  0

   
E h   0  R  y  dy
h
  ym (h)  1  1   b  Équation A. 6. 1
 R R0   h

0 E  y  dy
h
E h est la moyenne d’ E sur l'épaisseur h : E h


0 = N(R) + - N(R)
0 = M(R) + -M(R)
Figure A. 6. 2. : Schéma d'un principe de raisonnement possible : la variation de courbure est
induite par le moment lié aux contraintes résiduelles dans la lamelle [h, h0] sur la nouvelle poutre
obtenue en enlevant de la matière.

L’équation du moment fléchissant par unité d'épaisseur :

0 E  y  0  y   a y  b  y  ym  h   dy  0
h

E  y   0  y   y  ym  h   dy   0 E  y   0  y   y  ym  h   dy
h0 h
0 h

  E  y   a y  b   y  ym  h   dy  0
h
0

Le premier terme de cette équation est nul. En effet :

E  y   0  y   y  ym  h   dy 0 E  y  0  y   y  ym h   ym h0   ym h   dy
h0 h0
0  0

  E  y   0  y   y  ym  h0   dy   ym  h0   ym  h    E  y   0  y  dy
h0 h0
0 0

 M 0   ym  h0   ym  h   N0  0  0  0
L’équation du moment devient donc :

68
 h
   R  y   y  ym  h   dy  b  E  y   y  ym  h   dy
h0
h
0


a  E  y  y y  ym  h  dy  0 
h
0


   R  y   y  ym  h   dy  b  0  a  E  y  y y  ym  h  dy  0 
h0 h
h 0


   R  y   y  ym  h   dy  a  E  y  y y  ym  h  dy  0 
h0 h
h 0

0 E  y  y  y  ym  h  dy
h
Posons (h)  de sorte que l'on identifie ici le correspondant du
classique produit (moment quadratique x module d'Young) de la théorie des poutres à élasticité
uniforme.

a (h)    R  y   y  ym  h  dy
h0
Équation A. 6. 2
h

La valeur moyenne <E>h, entre 0 et h, du module E peut être mesurée par un essai de
traction (mesure globale). En enlevant une couche de matière ou lamelle de hauteur (h0 – h) finie,
nous n'aurons accès qu'à une moyenne de la contrainte résiduelle dans cette lamelle. Par souci de
simplicité, nous considérons que  R  y  est constante - et implicitement E(y) aussi - sur
 h  h0 
l’épaisseur de la lamelle enlevée h, h0  , i. e.  R  y    R   pour y  h, h0  .
 2 

 1 1   
 h h  
L’équation A. 6. 1. devient  ym (h)  
 
R R 
 b

E h R 

0  h h
2  0
 
  0   
1 1    h0  h     h0  h  
 R R  (h)   R  2   h0  h   2  ym  h  
L’équation A. 6. 2. devient a (h)   
 0    
Ceci valide dans le cas d'un module variable avec y le schéma de pensée classique de la
Figure A. 6. 1. Pour continuer le calcul, nous devons donc connaître E  y  . Avec les nouvelles
technologies de micro et nanoindentation nous pourrons bientôt mesurer le module E après
l’enlèvement de chaque lamelle. L’évolution d’ E avec l’épaisseur donnée par ces mesures de
nanoindentation plus la valeur moyenne d’ E mesurée par un essai de traction permettent de
remonter à la fonction E  y  et de calculer la contrainte résiduelle par les précédentes équations
A. 6. 1. et A. 6. 2.

A6. 4 Cas "classique" du module d'Young uniforme

Dans le cas où E  y  est constant par rapport y alors on retrouve bien la formule de Stoney :
h2
y E  y  dy
h

ym  h  
0  2 h
E

E  y  dy
h
Eh 2
0

69
   h  h 
  h0  h  0 2  aE h 12
L'équation A. 6. 2. devient  R  0
h 3

 2 

  h0  h   a E h3 E h2  1 1 
 R       
 2  6  h0  h  h0 6  h0  h   R R0 

A6. 5 Moyenne du module d'Young mesuré par essai de flexion

Si l'on cherche à mesure le module moyen par flexion, <E>(I, h), on retrouve la formule de
h
y(y  ym(h)) dy. La moyenne a donc un sens différent de celle
(h) et <E>(I, h) = (h)/Ih où Ih = 
0
h
1
mesurée en traction uniaxiale notée <E>h qui vaut 
h
E(y) dy.
0

A6. 6 Contrainte résiduelle éventuellement additionnelle liée à l'enlèvement de

lamelle
En pratique on peut enlever des couches par fraisage. Le fraisage induit un traitement
thermique en surface. En effet, l'usinage ou la coupe perturbe l'état de contraintes résiduelles sur
une faible épaisseur en surface. Si cela n'est pas négligeable, il convient alors de considérer qu'une
couche surfacique est perturbée par le fraisage avec ses propres contrainte résiduelle et module.
Cependant, nous avons fait l’hypothèse que l’enlèvement d’une lamelle ne perturbe pas la
distribution des contraintes résiduelles dans le reste de la plaque. Introduisons cette perturbation P u
en faisant l'hypothèse qu'elle est constante, concentrée en y = h (donc Pu s'exprime en N/m et R en
N/m²), et vient se superposer à R. L'équation A. 6. 2. devient :
1 1    h0  h     h0  h  
a (h)     (h)   R    h0  h    ym  h    Pu  h  ym  h  
R R0   2   2 
Si cette perturbation est considérable alors il conviendra de corriger les rayons de courbure.
S'il est possible d'effectuer un traitement thermique permettant de relaxer les contraintes internes
alors Pu peut-être estimé. Afin de vérifier que cette perturbation est bien considérable, on peut, par
exemple, proposer d'enlever des lamelles après traitement thermique. Dans un premier temps, si
des lamelles sont enlevées en fibres supérieure et inférieure de façons non symétriques alors une
perturbation Pu considérable entraînera des courbures de la plaque et il conviendra de se soucier de
prendre en considération Pu. Dans le cas contraire, il semblera inutile de se préoccuper de cette
perturbation. Si Pu est considérable, elle sera mesurée par enlèvements successifs de lamelles à
l'aide des équations précédentes. Bien sûr, il conviendra de vérifier que le traitement thermique ne
modifie pas de façon conséquente le module d'élasticité local du matériau.

A6. 7 Références bibliographiques de cette annexe

M. Alchikh, C. Fond, Y. Frere et H. Pelletier, "Mechanochemical degradation of poly (vinyl
fluoride) by Sodium Hydroxide measured by Microindentation", Journal of Materials
Science, (2010), 45, p. 2311-2316.
M. Alchikh, C. Fond et Y. Frere, "Discontinuous crack growth in poly (vinyl fluoride) by
mechanochemical ageing in Sodium Hydroxide", Polymer Degradation and Stability,
(2010), 95, p. 440-444.
Mohamad Alchikh, "Effet du vieillissement chimique dans la soude sur le comportement
mécanique du PVDF", thèse de l'Université de Strasbourg, soutenue le 26 mars 2009.

70
Laroze, S., 1988, Mécanique des structures - Tome 2 : Théorie des poutres, 2 éd.,
Eyrolles/Masson.
Stoney, G. G., 1909, The tension of metallic films deposited by electrolysis: Proc. Roy. Soc.,
v. A82, p. 172.

71