REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

ECOLE NATIONALE POLYTECHNIQUE DEPARTEMENT DU GENIE CIVIL

5ime anne

Mini projet de PONT

LA METHODE DE GUYON-MASSONNET-BARE

Ralis par :

TERMOUL Tarek

2009/2010

PONTS Introduction

5ime anne 2009/2010

La mthode de GUYON-MASSONNET-BARE repose sur la thorie des plaques orthotropes. Elle fut dveloppe par Guyon [GUY.46] dans le cas d'une dalle orthotrope rigidit torsionnelle ngligeable. Massonnet en 1950 gnralisa les relations trouves par Guyon en introduisant l'effet de la torsion dans les calculs [MAS.50]. En 1966, Massonnet et Bare publirent un recueil de ces mthodes illustr par un nombre d'exemples [MAS.66]. Cette mthode est simple et utilise pour le calcul des ponts poutres multiples en tenant compte de la rsistance la torsion du pont, elle donne des rsultats satisfaisants par rapport aux autres mthodes, car la rigidit transversale du tablier n'est pas ngligeable.

Principe de la mthode et hypothses

Cette mthode vise dterminer les efforts transitant dans un grillage de poutres soumis un chargement quelconque, ponctuel ou rparti. Le systme dalle-poutre discret est remplac par un systme uniforme compos d'une dalle anisotrope ou orthotrope ayant des caractristiques constantes suivant chacun de ses axes transversal et longitudinal. Ce passage d'une rpartition discrte de la rigidit, une rpartition continue, est l'hypothse principale sur laquelle repose cette mthode. La deuxime hypothse consiste admettre que le coefficient de Poisson du matriau constitutif est nul. Cette hypothse est plus ou moins contestable, mais, dans la mesure o le but de la mthode est de dterminer la rpartition des efforts dans les diffrentes parties de la structure et o les variations de ces efforts ne sont pas trs importantes, l'erreur qui en rsulte peut tre considre comme ngligeable. Le rseau de poutres est assimil une dalle orthotrope possdant deux bords libres (selon Ox) et deux bords simplement appuys (selon Oy). La mthode s'appuie sur la rsolution approche de l'quation diffrentielle d'un grillage simple constitu, dans le sens y, de m poutres (longerons) espaces les unes des autres de b0, et dans le sens x, de n entretoises espaces de l0; les poutres sont, dans les deux sens, faites du mme matriau, rsistantes la torsion et assembles aux nuds rigidement l'une l'autre.

ENP. Dpartement G.CIVIL

PONTS

5ime anne 2009/2010

Dsignons par BP = EIP la rigidit flexionnelle des poutres BE = EIE la rigidit flexionnelle des entretoises et par CP = GIPt et CE = GIEt les rigidits la torsion des poutres et entretoises avec G= E/2(1+) ; coefficient de poisson. Si nous remplaons cette construction par une construction quivalente o les rigidits sont rparties continument sur la longueur et la largeur, les rigidits flexionnelles seront par unit de longueur : de mme, les rigidits torsionnelles par unit de longueur de la construction quivalente seront
et = = et = =

Myx

Mxy

Pour un dplacement vertical (la dforme de la dalle) de la construction selon la surface W = w(x,y) dont les courbures ont pour valeurs flchissants par unit de longueur

et

, il nat des moments

=
l'expression

et =

La gomtrie montre que la torsion de la surface de dplacement W est donne par et dans la construction natront donc les moments de torsion unitaires

L'effort tranchant vertical Qy se calcul suivant la formule = remplaant les valeurs de Mx et Mxy :

et =

, en

De la condition d'quilibre de l'lment dans le sens vertical on trouve aprs remplacement de Mx ,My , Mxy , Myx par leurs valeurs, l'quation diffrentielle d'un grillage simple dont les rigidits sont rparties continument

+ +

p( x,y ) est le chargement de la dalle. Cette quation est quivalente celle d'une plaque orthotrope.

= ,

ENP. Dpartement G.CIVIL

PONTS

5ime anne 2009/2010

Les deux paramtres fondamentaux (paramtres sans dimension)

paramtre d'entretoisement : a t dduit par Guyon en calculant les grillages sans tenir compte de la torsion :

= .

Avec : b : demi largeur de la dalle l : porte de la trave Le paramtre dtermine la souplesse de l'entretoisement, plus grand est , plus souple est l'entretoisement.

paramtre de torsion :
Dans la pratique, le coefficient de rigidit torsionnelle (P+E) est toujours compris entre les valeurs correspondant aux deux cas particuliers, celui de la dalle et celui du grillage simple. Si l'on pose + = 2 . l'effet de torsion est caractris par le paramtre dont la valeur, pour couvrir le domaine entier entre les deux cas particuliers prcits, variera de 0 1. = avec 0 1 . . Dans le cas gnral, nous pouvons calculer, par cette formule, le paramtre de torsion mais l'valuation des rigidits la torsion P et E tant, ordinairement, trs difficile, il faut introduire souvent des hypothses simplificatrices pour obtenir une valeur approche de . Le paramtre de torsion prend une expression particulirement simple dans le cas d'une construction mixte. En effet, nous pouvons admettre que la rigidit propre de torsion des poutrelles mtalliques est ngligeable et assimiler par consquence le pont une plaque dont les rigidits la flexion dans les deux sens p et E sont celle de la dalle isotrope en bton majores dans les rapports =

et =

avec =

et ed la hauteur de la dalle en bton). Tout calcul fait, on obtient pour un pont mixte :

(Ed est le module de Young

Par ces deux paramtres et , le comportement de la construction est compltement dfini. La rsolution analytique directe de l'quation diffrentielle de la plaque orthotrope conduit des calculs compliqus et peu pratiques mettre en uvre. La mthode de Massonnet permet de s'affranchir de cette difficult en utilisant une mthode approximative base sur les coefficients de rpartitions.

ENP. Dpartement G.CIVIL

PONTS Mthode des coefficients de rpartition

5ime anne 2009/2010

Deux hypothses servent de base la configuration de la mthode : la construction relle est remplace par une dalle orthotrope prsentant les mmes rigidits moyennes de flexion et de torsion et qui est au sens technique exactement soluble par le calcul diffrentiel. la rpartition transversale relle du chargement est remplace par celle qui nat sous une charge rpartie le long de l'axe X de la construction et d'excentricit e suivant la loi sinusodale de la forme = . sin o p1 =la valeur constante de chargement.

Les hypothses cites n'influencent que le calcul de la rpartition transversale de la charge, les autres calculs obissent aux rgles ordinaires de la stabilit des constructions. Le coefficient de rpartition transversale La construction prend une dforme en demi-onde de sinusode selon l'quation , = sin

Si la charge p(x), au lieu d'tre rpartie sur une droite, est rpartie uniformment sur la largeur 2b de la construction (tout en restant sinusodale dans le sens longitudinal), la construction prend une dforme en surface cylindrique d'quation =
2 1 2

x b

L /4

L /4

L /4

L /4

Le moment par unit de largeur pour un chargement de cette nature est

ENP. Dpartement G.CIVIL

PONTS

5ime anne 2009/2010

Dsignons le rapport du dplacement vertical w(x,y) d'un point de la construction sous l'effet d'une charge linaire p(x) celui w0(x) du mme point mais sous l'effet de la charge p0(x) uniformment rpartie sur la largeur du pont, par le coefficient de rpartition transversale K(y) 2 Ce coefficient K dpend du paramtre et , de l'excentricit relative e/b de la charge linaire et de l'ordonne relative y/b du point considr de la construction. =

Le calcul de K partir des relations complexes d'intgral est en gnral difficile c'est pourquoi on divise la construction dans le sens de la largeur en 8 bandes de mme largeur pour simplifier l'intgral. On trouve que le coefficient K dpend, entre autre, de la valeur du paramtre , pour viter de calculer sparment K pour chaque valeur de on utilise les formules d'interpolation dans lesquelles on emploie les coefficients K0 et K1 pour les valeurs extrmes = 0 et = 1 qui sont dfinis dans des tableaux pour des valeurs de donnes dans les abaques de Bare et Massonnet. Sur la base de calculs d'un grand nombre de cas, Massonnet a dtermin dans ces tableaux les valeurs des coefficients K0 et K1 pour compris entre 0.05 et 5.00 et selon les diffrentes excentricits de charges (e = b; 3b/4; b/2; b/4; 0) et pour les sections de la largeur de la dalle ( y = 0; y = b/4; y = b/2; y = 3b/4; y = b ) . Pour un calcul rigoureux de K , il est ncessaire d'appliquer les formules d'interpolation tablies par Sattler et qui dpendent aussi de la valeur de : Si 0 < 0.1 Si 0.1 < 1 >1 Si = + . . = + . = + . .

O = 1 ./. Dans le cas o calcule ne figure pas dans le tableau on doit faire une interpolation, et si la poutre en question se trouve entre deux sections dont les lignes d'influence sont connues on peut aussi faire une interpolation.

Calcul des efforts

Chaque type d'effort (moments, efforts tranchants, etc.) fait intervenir un coefficient de rpartition transversale des charges diffrent. Ce dernier est alors multipli par le moment moyen pour obtenir le moment flchissant existant dans une poutre dtermine. Le moment moyen correspond au moment de flexion de la poutre seule sous la charge et appuye ses extrmits.

ENP. Dpartement G.CIVIL

PONTS

5ime anne 2009/2010

On dplacera les charges de faon obtenir les plus grandes ordonnes et on retiendra pour le calcul des efforts l'excentricit qui donne les plus grandes valeurs des coefficients. Dans le cas des ponts poutres multiples la section d'tude sera impose par la position de la poutre, ce qui nous amne tracer les lignes d'influences pour les diffrentes excentricits de charge et on retiendra la section qui donne les plus grandes valeurs des coefficients.

Calcul des moments flchissant longitudinaux (dans les poutres)

Le rapport du moment flchissant rel Mx(x,y) au moment moyen M0(x) dans le sens transversal est gal au coefficient de rpartition transversale

Ky =

Le moment flchissant dans le sens X par unit de largeur pour une charge sinusodale est:

, =

Pour les charges , , , , ,

Le moment flchissant moyen la distance x est = (1) et M0 (x) reprsente le moment flchissant moyen produit par une charge linaire sinusodale unitaire, uniformment rpartie sur la largeur de la construction.

Ky =

Pour obtenir en un point quelconque de la construction le moment flchissant rel il suffit de calculer en ce point le moment moyen M0(x) et de le multiplier ensuite par le rapport:

Pour une charge linaire sinusodale l'expression du moment longitudinal est:

, = .

ENP. Dpartement G.CIVIL

PONTS =

5ime anne 2009/2010

Pour une charge sinusodale distribue sur la surface du tablier : , = . 0 t et reprsente respectivement les surfaces totale et partielle sous le chargement, du diagramme de K(y). Si une charge concentre unitaire applique un point gnral de coordonnes x=d, y=e se dplace sur la construction, l'quation de la surface d'influence des moments flchissant longitudinaux est

, =
Exemple

l'indique la figure, aux points e=-b, e=b/4, e=3b/4, on doit dterminer les flches et les moments flchissant aux sections transversales x=l/2 et x=l/4. p1 =1t/m, P =30000tm2/m, E =2845tm2/m.

Un pont est charg d'une charge linaire = comme

= 0,49 = . = = 0,72

La flche moyenne selon l'quation de la dforme est:

..30000

= 0,000427

(m)

Le moment flchissant longitudinal moyen est:

= 1,264 sin

(t.m)

Par interpolation linaire dans les tableaux de Massonnet pour =0,70 et =0,75 on trouve les valeurs des coefficients K0 et K1 pour =0,72. Dans la colonne y=b/4 on lit pour e=-b, b/4, 3b/4:

ENP. Dpartement G.CIVIL

PONTS

5ime anne 2009/2010

=0 =0 =0 =1 =1 =1

= 0,70 = 0,75 = 0,72 = 0,70 = 0,75 = 0,72

e = -b - 0,3589 - 0,4324 - 0,3883

e = b/4 1,7118 1,8140 1,7528

e = 3b/4 1,1934 1,1305 1,1683

0,4808 0,4351 0,4625

1,3307 1,3825 1,3514

1,1589 1,1584 1,1587

K (=0,49)

0,207

1,472

1,161

Les flches et les moments moyens sont: en = = 0,000427 m en =

= 0,000302 m

et = 1,264 t.m et = 0,895 t.m

Donc on obtient les flches et les moments longitudinaux en multipliant les flches et les moments moyens par le coefficient K calcul: e = -b x = l/2 w mm Mx tm w mm Mx tm 0,088 0,262 0,062 0,185 e = b/4 0,630 1,861 0,445 1;317 e = 3b/4 0,496 0,470 0,351 1,039

x = l/4

ENP. Dpartement G.CIVIL

PONTS

5ime anne 2009/2010

Calcul des moments flchissant transversaux (dans les entretoises)

Le moment flchissant par unit de largeur d'entretoise c'est--dire par unit de longueur dans le grillage est exprim, pour une charge sinusodale, par la relation:

, = . . .

Le coefficient (, , e, y) dpend des valeurs du paramtre d'entretoisement , du paramtre de torsion , de l'excentricit e de la charge et de l'ordonne y de la section faite dans l'entretoise dans laquelle on cherche l'effet.

Comme pour le coefficient K on dtermine pour quelconque par les formules interpolation Si 0 < 0.1 Si 0.1 < 1 Si >1 = + . . = + . = + . .

O = 1 ./. 0 correspondant = 0 1 correspondant = 1 (les valeurs de 0 et 1 sont donne dans les tableaux) Si la construction est soumise un systme de n charges linaires sinusodales: , , , , , ; le moment existant dans la section y de l'entretoise s'obtient par superposition des effets individuels

, = .

Pour une charge unitaire localise une distance d des appuis se dplace sur la construction, la surface d'influence des moments transversaux est

, =

ENP. Dpartement G.CIVIL

PONTS Calcul du moment de torsion dans le grillage:

On sait que : = et =

5ime anne 2009/2010

La diffrence de ces moments vaut : = + En introduisant les relations suivantes :

+ = 2

et

, =

Dans lquation prcdente, nous trouvons : = 2 . Pour l'emploi pratique, il est avantageux d'exprimer la diffrence des moments par la relation:

= 2

m est le coefficient sans dimension qui dpend des valeurs du paramtre d'entretoisement , du paramtre de torsion , de l'excentricit e de la charge et de l'ordonne y du point cherch. = + = + Massonnet a trouv, en calculant ces coefficients pour diffrentes valeurs de entre 0 et 1, que pour dterminer l'expression gnrale de on pouvait employer la relation simple = La valeur de 1 peut tre lue directement dans les tables et comme prcdemment.

Utilisation pratique de la mthode de Guyon-Massonnet:

On peut appliquer cette mthode dans les cas suivants : Un tablier de pont poutres sans entretoises intermdiaires cest alors hourdis qui joue le rle des entretoises. Dans ce cas, pour calculer la rigidit de torsion dune poutre comprenant une certain de largeur de hourdis, on doit : Dcomposer la section de la poutre en lments rectangulaires Faire la somme des inerties de torsion de ces lments. Attribuant la fraction de table de compression appartenant au hourdis un moment dinertie de torsion gale la moitie des moments dinertie calcule comme si ctait une poutre. Ossatures mixtes comportant un membre de poutrelles suprieures 2, dans ce cas, on doit admettre que linertie de torsion des poutrelles (et des entretoises ventuelles) est nulle. Les ponts biais, condition que le biais soit suprieure a 65.

ENP. Dpartement G.CIVIL

10

PONTS
Remarque:

5ime anne 2009/2010

Pour une dalle plaine doublement prcontrainte, hypothse de mellite du coefficient de poisson peut tre contestable. Pour le bton, ce coefficient est de lordre de 0,15. Rowe a montr que les coefficients K0, K1, et

ne sont gure affects, par contre, le

coefficient 1 est beaucoup plus sensible, en dfinitive, il a propos de calculer le coefficient

( )0.15

correspondant = 0.15 par la formule suivante :

( )0.15 = ( 1 )0.15 + ( ( 1 )0.15 0 ) (1)0.15 Est donn par des tables spcifiques.

Commentaires sur cette mthode

Comme mentionn prcdemment, la mthode de Massonnet est simple utiliser grce aux diffrentes tables permettant de calculer les coefficients de rpartitions transversales. Les raideurs des poutres et des entretoises sont tartines sur la longueur et la largeur de la dalle quivalente. L'excentrement des poutres et des entretoises par rapport la dalle est nglig. Cette mthode est aussi efficace que celle de Courbon dans le cas d'un calcul de prdimensionnement mais elle ncessite un apprentissage plus important. De plus, les entretoises ne sont pas considres comme infiniment rigides. Elle est particulirement performante dans le cas d'un trs grand nombre de poutres et d'entretoises puisque l'erreur commise par le tartinage des raideurs est alors plus faible. Cependant, avec les conceptions modernes des ponts mixtes, le nombre d'entretoises est souvent assez faible. La mthode de Massonnet est alors moins performante pour ces ponts. Cependant, elle permet de prvoir le comportement du pont sans entretoise, en considrant uniquement la raideur de la dalle dans la direction transversale et celle des poutres mixtes dans la direction longitudinale.

Conclusion
La mthode de Guyon-Massonnet-Bare permette de raliser un prdimensionnement des structures poutres et dalle. Elle est beaucoup t employe avant les progrs informatiques des annes 70-80 car elle offrait des moyens simples pour calculer les sollicitations et les flches dans ces structures. Elle est performante pour des structures comportant un grand nombre de poutres et d'entretoises. Elle permette de rsoudre les quations aux drivs partielles des plaques en passant par des dcompositions comme celles des sries de Fourier. Cependant, elles reposent sur des hypothses fortes et la gomtrie de l'ouvrage n'est jamais compltement respecte. Elle fut plus ou moins abandonne lors de la conception de logiciels bass sur la mthode des lments finis. Ces logiciels ont suivi le dveloppement des capacits des ordinateurs,

ENP. Dpartement G.CIVIL

11

PONTS

5ime anne 2009/2010

grant de plus en plus d'quation simultanment et de plus en plus vite. Dans les annes 7080, la capacit des ordinateurs tait encore trs faible et la mthode des lments finis ne pouvait s'appliquer qu' des lments linaires et peu nombreux. C'est pourquoi, les premiers calculs de ponts multi poutres entretoises utilisant la M.E.F ncessitaient la traduction gomtrique de l'ouvrage en grillage de poutres comportant la raideur de la dalle. Bibliographie [1]- CHARLES MASSONNET et RICHARD BARE Le calcul des grillages de poutres et dalle orthotropes selon la mthode Guyon-Massonnet-Bare , DUNOD Paris 1966. [2]- Yannick SIEFFERT l'entretoisement des ponts mixtes multipoutres ferroviaires INSA de Lyon 2004.

ENP. Dpartement G.CIVIL

12