Calcul des ponts biais à poutres multipes sans

entretoises par la méthode des éléments finis

Autor(en): Jirousék, Jaroslav

Objekttyp: Article

Zeitschrift: Bulletin technique de la Suisse romande

Band (Jahr): 98 (1972)

Heft 10

PDF erstellt am: 08.02.2016

Persistenter Link: http://dx.doi.org/10.5169/seals-71546

Nutzungsbedingungen
Mit dem Zugriff auf den vorliegenden Inhalt gelten die Nutzungsbedingungen als akzeptiert.
Die ETH-Bibliothek ist Anbieterin der digitalisierten Zeitschriften. Sie besitzt keine Urheberrechte an
den Inhalten der Zeitschriften. Die Rechte liegen in der Regel bei den Herausgebern.
Die angebotenen Dokumente stehen für nicht-kommerzielle Zwecke in Lehre und Forschung sowie für
die private Nutzung frei zur Verfügung. Einzelne Dateien oder Ausdrucke aus diesem Angebot können
zusammen mit diesen Nutzungshinweisen und unter deren Einhaltung weitergegeben werden.
Das Veröffentlichen von Bildern in Print- und Online-Publikationen ist nur mit vorheriger Genehmigung
der Rechteinhaber erlaubt. Die Speicherung von Teilen des elektronischen Angebots auf anderen
Servern bedarf ebenfalls des schriftlichen Einverständnisses der Rechteinhaber.

Haftungsausschluss
Alle Angaben erfolgen ohne Gewähr für Vollständigkeit oder Richtigkeit. Es wird keine Haftung
übernommen für Schäden durch die Verwendung von Informationen aus diesem Online-Angebot oder
durch das Fehlen von Informationen. Dies gilt auch für Inhalte Dritter, die über dieses Angebot
zugänglich sind.

Ein Dienst der ETH-Bibliothek

ETH Zürich, Rämistrasse 101, 8092 Zürich, Schweiz, www.library.ethz.ch

http://retro.seals.ch
BULLETIN TECHNIQUE DE LA SUISSE ROMANDE 98- année 13 mai 1972 N° 10

COMMUNICATION DE LA CHAIRE DE STATIQUE ET DE RÉSISTANCE DES MATÉRIAUX DE L'EPFL

Calcul des ponts biais à poutres multiples sans entretoises par

la méthode des éléments finis
par J. JIROUSEK, Dr se. techn., ing. SIA'

Dotée pendant très longtemps de moyens relativement modestes, la Chaire de statique et

de résistance des matériaux de VEcole polytechnique de Lausanne s'est consacrée essentielle¬
ment à sa tâche prioritaire, l'enseignement, laissant par la force des choses la recherche au
second plan. Bénéficiant maintenant de possibilités plus étendues, elle a eu le bonheur de pou¬
voir s'attacher M. J. JirouSek qui, après avoir enseigné dans plusieurs universités étrangères,
lui apporte le concours de ses connaissances et de sa compétence, qui sont grandes.
Nous sommes heureux de présenter ici cette publication ; si elle représente pour son auteur
la suite d'un travail commencé ailleurs, elle est pour notre chaire le témoignage d'une activité
que nous espérons promise à de nouveaux développements.
Professeur M.-H. Derron.

Introduction Les méthodes de calcul des ponts multiples à charnières

Le pont à poutres multiples sans entretoises est formé longitudinales sans frottement peuvent être réparties en
trois catégories principales. La première, la plus fréquente,
par un système de poutres préfabriquées placées l'une à est la théorie de la plaque articulée (voir par exemple [20],
côté de l'autre et liées le long de la portée (voir, par exemple,
[25] s) ; elle comprend toutes les méthodes qui ana¬
fig. 1). L'obliquité de l'ouvrage en plan peut être obtenue
lysent l'ouvrage approximativement comme un conti¬
par un décalage des extrémités des poutres. La liaison
entre les poutres est réalisée par le remplissage des loge- nuum sans rigidité transversale de flexion. La seconde
ments (shear key) aménagés le plus souvent dans la partie
catégorie comprend les méthodes dès déformations compa¬
tibles (voir par exemple [1], [4], [12], [14] ; elles ana¬
supérieure des joints (fig. 2). De plus, dans certains cas,
lysent le système des poutres en se basant sur les principes
l'ouvrage peut être aussi précontraint latéralement. Si la
continuité ainsi obtenue est suffisante, on peut calculer généraux de la méthode des forces. La troisième catégorie,
celle des méthodes d'équilibre des joints (voir [9], [10]
approximativement l'ouvrage comme une dalle ortho-
adopte au contraire les principes généraux de la méthode
trope. Par contre, si la précontrainte latérale est faible ou
des déformations.
inexistante, les déformations de l'ouvrage et le retrait peu¬
vent entraîner un décollement des clavettes de remplissage
des joints. La rigidité transversale, et par conséquent aussi
1
Chargé de cours, prem ier assis tan t à la chaire de statique et
de résistance des matériaux de l'EPFL.
la possibilité de transmission des moments transversaux 2 National
Cooperative Highway Research program Report 83
disparaissent. On peut alors considérer que les poutres (1970).
sont liées le long de la portée par des charnières longitudi¬ 8 Les chiffres
entre crochets renvoient à la bibliographie en
nales sans frottement. Cette hypothèse, très généralement fin de l'article.
admise en pratique, a été adoptée par exemple aux Etats-
Unis par le « Highway Research Bord » (voir le récent
Niveau théorique Remplissage
rapport n° 83 2). C'est également ce que nous ferons dans des charnières (Shear Key J
l'étude qui va suivre.

A-A

¦
G Zs
.igne aeppu
E
n
G
„igne a appui
D
Èf

BEBE G > centre de gravité de la section

y i—i^t S
¦ centre de cisaillement de la section
Fig. 2. — Géométrie de la section charnières fictives formées
Fig. 1. — Pont biais à poutres multiples. par le remplissage des joints.

169
A cause des complications mathématiques, les méthodes f) les appuis sont constitués de façon à ne pas entraver
connues se limitent presque en totalité à des ponts rectan¬ le gauchissement des sections d'appui.
gulaires. De plus, on suppose aussi presque toujours que
toutes les poutres ont la même section. D'après ce que nous D'autres hypothèses auxiliaires seront formulées au
savons, les travaux concernant les ponts biais se limitent à cours du développement de la méthode de calcul.
ceux (voir [9], [10] et [16]) qui ne traitent que le cas où les
charnières entre les poutres sont au niveau des centres de
cisaillement des sections. Ils ont, néanmoins, clairement 2. Efforts internes des poutres et énergie potentielle
démontré l'effet non négligeable du biais sur l'interaction de déformation
des poutres du pont.
Le but de la présente étude est de montrer les possibilités Considérons une poutre de l'ouvrage (voir par exemple
d'application de la méthode des éléments finis au calcul la figure 3, où G désigne le centre de gravité de la section
des ponts biais à poutres multiples. Vu l'utilisation crois¬ et S son centre de torsion). Dans le cas particulier où S
sante des ordinateurs dans la pratique, on peut envisager se confond avec G et où les charnières entre les poutres se
que l'élaboration d'un programme standard de calcul des trouvent au niveau des centres de torsion et si les sections
ponts biais (en général avec des poutres de sections ne se gauchissent pas, les charges verticales appliquées sur
inégales) permettrait de résoudre la majorité des pro¬ l'ouvrage n'engendrent qu'une flexion verticale et une tor¬
blèmes. sion des poutres autour de S. Alors, les points du plan
Afin d'obtenir une solution suffisamment simple, il a horizontal passant par les charnières ne se déplacent pas
fallu formuler certaines hypothèses simplificatrices. Celles- dans ce plan (« v 0) et les forces d'interaction d'une"
ci sont en l'occurrence toujours très bien satisfaites pour poutre sur l'autre, par l'intermédiaire des charnières, sont
les poutres-caissons, mais le sont en général beaucoup toutes verticales. Par contre, dans le cas général où ces
moins pour les poutres à section ouverte. La méthode de conditions simplificatrices ne sont pas réalisées, les points
calcul développée dans ce travail s'applique donc essen¬ du plan passant par les charnières subissent des déplace¬
tiellement aux ponts à poutres de section fermée. ments horizontaux et les forces d'interaction ont, en plus
La technique de calcul par éléments finis utilisée ici se des composantes verticales, des composantes horizontales
base sur la subdivision en plan de chaque poutre en une latérales et longitudinales (respectivement perpendicu¬
série de macro-éléments trapézoïdaux et sur les fonctions laires et parallèles aux axes des poutres). Ainsi, les poutres
de déplacement de ces> éléments. Ces fonctions, exprimées sont en général soumises non seulement à la flexion verti¬
par rapport aux paramètres nodaux des angles de l'élé¬ cale et à la torsion, mais aussi à la flexion horizontale et
ment, assurent aussi bien la compatibilité des déformations aux efforts normaux.
des éléments de la même poutre, que celle des poutres elles- La figure 3 montre les efforts internes de la section
mêmes aux joints. Partant du théorème du minimum de transversale d'une poutre dans ce cas général. On désigne
l'énergie potentielle, on développe la matrice de rigidité de par My, Mz les moments fléchissants, par T le moment de
l'élément liant les forces nodales aux déplacements nodaux. torsion, par Qy, Qz les efforts tranchants et par TV l'effort
Ceux-ci sont ensuite calculés à partir des conditions d'équi¬ normal. Les composantes des déplacements selon les direc¬
libre aux nœuds de l'assemblage. tions du système global d'axes x, y, z sont respectivement
m, v et w (fig. 3). De plus, y/ désigne la rotation autour de
l'axe de la poutre.

1re partie: Théorie

1. Formulation du problème et hypothèses de base
k^-Z
Considérons un pont biais à poutres multiples sans
entretoises (fig. 1). Les poutres peuvent être de section X,u
différente mais constante le long de la portée. L'ouvrage
n'est pas précontraint latéralement, ou cette précontrainte
est négligeable. Il s'agit d'étudier l'interaction des poutres y.v
sous l'effet des charges verticales et du tassement des
appuis, et de développer une méthode pratique de calcul. z,w

HypothiÉps de base :

a) les sections transversales des poutres possèdent un

axe vertical de symétrie ;
b) les liaisons entre les poutres sont équivalentes à des CjJ>
charnières longitudinales travaillant sans frottement ;
c) toutes les charnières liant les poutres entre elles se
trouvent dans le même plan horizontal (fig. 2) ;
d) l'obliquité en plan de l'ouvrage n'est pas excessive,
e
de sorte que les dimensions des sections biaises des
poutres, parallèlement aux lignes d'appui, sont
Y
B
petites vis-à-vis de la longueur des poutres ;
e) la section transversale des poutres est indéformable Fig. 3. — Composantes des déplacements («, v, w, y/) et efforts
dans son plan ; internes d'une poutre du pont.

170
Avec les conventions de la figure 3, les moments de horizontale des poutres. Cette flexion provoque le long des
flexion, le moment de torsion et l'effort normal s'expriment poutres des déplacements horizontaux vs vs (x), diffé¬
en fonction des déformations par rents d'une poutre à l'autre. Mais, puisqu'on suppose en
même temps que les sections transversales des poutres
dws restent indéformables dans leur plan (les éléments de l'ou¬
M„ El, (la)
vrage travaillent comme des poutres et non pas à la manière
d'une ossature plissée), on doit admettre que dans le plan
M, EIZ '
(lb) horizontal passant par les axes des charnières, les déplace¬
dx2
ments v de toutes les poutres doivent être égaux. Puisque
l'assemblage des poutres est toujours très rigide dans ce
GJ
dy
-EC dy (le) plan, il est logique d'admettre que les déplacements v au
~cbc~ d? niveau des charnières sont négligeables vis-à-vis des déplace¬
d«g ments verticaux w. Il en résulte (voir la figure 4) que le
N= EA dx
' (ld) déplacement horizontal v$ du centre de rotation de n'importe
quelle poutre est dû uniquement à sa distorsion y/ et s'ex¬
D'autre part, le bimoment dû aux contraintes normales prime simplement sous la forme
de torsion non uniforme (torsion fléchie) vaut (3)
vs ¦zsy/'

B= -EC dy
d?' (le)

Dans ces relations,

4r -JÜ
Iy> h moment d'inertie de la section par rapport à
l'axe principal parallèle respectivement äyetz
et passant par G [cm4] ;
r
r 11

p. 1z>
J moment de rigidité de torsion uniforme
(Saint-Venant) [cm4] ; J
C moment d'inertie sectorielle [cm6] ;
A aire de la section [cm2] ;
ws> vs — déplacements latéraux ve et v du centre de
torsion de la section ;
uq déplacement longitudinal « du centre de gra¬
vité de la section. =-Zs-4'
En négligeant l'effet des contraintes tangentielles dues Fig. 4. — Relation entre le déplacement latéral v„ et la rota¬
à la flexion, ainsi que celles dues à la torsion fléchie, l'éner¬ tion i// de la section.
gie potentielle des efforts internes (énergie potentielle de
déformation) de la poutre s'écrit :
2. La solution se simplifie considérablement si l'on
dhvs d2^ admet que les charnières n'empêchent pas les fibres adja¬
Ut M„ '
dx2
Mz
d*2
+ centes de deux poutres voisines de glisser longitudinalement
o l'une par rapport à l'autre. (Ce glissement correspond à la
différence de déformations dues à l'effet simultané de la
_dy
T (dy/ +B ~~dx* "(£ dx flexion et du gauchissement). Ceci implique que l'on peut
supprimer le dernier terme de l'équation (2), puisque les
efforts normaux des poutres seront dorénavant nuls. La diffé¬
&Vs rence des déplacements longitudinaux des fibres adjacentes,
EL, EL +
dx" dx2 due à la flexion, est très petite, et tant qu'il ne s'agit pas
de sections ouvertes, il en est de même pour l'effet du gau¬
chissement. La confrontation des résultats numériques
+ Gjfê)' + EC dyy EA
dueV
dx. (2) donnés plus loin (2e partie, problème 3) montre que la
dx
possibilité de glissement longitudinal ne se manifeste pra¬
tiquement pas et que cette hypothèse simplificatrice peut
di//
Dans cette expression, Tg GJ -f- est la partie du être adoptée avec confiance pour l'analyse des ouvrages à
dx partir de poutres à section pleine ou fermée auxquelles est
moment de torsion selon Saint-Venant, tandis que la tor¬ consacrée cette étude. Elle serait, en revanche, beaucoup
sion fléchie s'exprime en fonction du bimoment B. moins satisfaisante pour les systèmes de poutres à section
ouverte, dont le gauchissement est beaucoup plus impor¬
tant.
3. Hypothèses simplificatrices
Avant de développer notre méthode de calcul, nous for¬
mulerons encore deux hypothèses simplificatrices, tirées de 4. Subdivision de l'ouvrage en macro-éléments et
l'analyse détaillée des propriétés des ouvrages en question : introduction des fonctions de déplacements
1. On a vu au paragraphe précédent que, dans le cas Pour analyser l'ouvrage par la métìiode des éléments
général, les charges verticales engendrent aussi une flexion finis, on le remplace par un assemblage de macro-éléments

171
I
2 3
on trouve

dx (\-rf)^'i + nw2, (4a)

J? —_ (Wl- w2)- tg a [(1 -rf)wi + nw£ (4b)

dy b

Pour t] 1
/2 on en tire
n-2 n-1
*>a (w)* ç2 =3(wi
2
+ wa) » (5a)

Fig. 5. — Assemblage d'éléments finis. dw\

dxli 2
i,
2
~ (5b)

selon la figure 5. Chaque poutre est, en général, subdivisée

dans le plan en une série de M macro-éléments trapézoï¬
W MphnJ — t(wi - Wjs)--tga.(w[ + w£ (5c)

daux de longueur
et avec (3) on a encore

/'-»¦ vs
zs 1

-rD (wi— W2) + Zr zs tg a (wx + w^. (5d)

On utilisera deux systèmes d'axes dans le plan (voir
fig. 5 et 6) : le système dit « global » formé des axes x, y
orthogonaux, et le système dit « local » formé par les axes Soit maintenant pour i 1, 2 les fonctions de déplace¬
biais je, tj, liés à un élément particulier. Le plan des deux ment des charnières pour un élément
systèmes correspond à celui des axes des charnières liant w{ wt (x) Bw«, (6)
les poutres (voir la figure 2). où
D'après la figure 6, les relations entre les coordonnées
globales du système et les coordonnées locales d'un élé¬
B =[/!(*) Â(x) fa(x) /,(*)] (6a)

ment s'expriment par est la matrice de quatre fonctions

x x0 + x + btj.tg a, 1

y y0 + bn. fi(x)= -^{Az-lx*A + 2x*),

On admettra que le déplacement vertical au niveau des
charnières (z 0) s'exprime : h &) -i ÇxA*- 2 x*A + x3)
w (x, rj) (1 — rf) W! (x) + jjw^ (je), (4) (6b)
fa(x)= -zQx*A-2x*),
où h>i (x) sont les déplacements verticaux des char¬
(x), m>2

nières entre les nœuds Au Bx respectivement An, B2. Pour 1

exprimer les dérivées de w (x, rf) par rapport à x et y, on Â(x) A" (- x*A + Xs)
écrit:
et
w,y X, w
WA{
dw x,n dw x, n
dx x,y x,y <Pa,
3 (6c)
x, t] g
où 9 sont les jacobiens de transformation. En dési¬
gnant pour simplifier
lu
est le sous-vecteur deppatre déplacements nodaux généra¬
d lisés.
dx
(...)' Les fonctions fix) ont été choisies afin de satisfaire aux
conditions

A (0) =1 Â (0) | 0 /8 (0) =0 h (0) 0

/i(0) -0 0 /^(0) 1 /8(0) 0 f'4(0) 0

y0[
AGd) Â(A) 0 Â(A)=l A(A) 0
fi<A) 0 /aW) 0 f'3(A) 0 f'iiA)^!
n=0 En même temps, les paramètres du nœud (égaux respec¬
WX.ri) tivement au déplacement w et à la rotation q> de la char¬
nière en ce point) seront toujours communs à tous les
éléments autour du même nœud. Ainsi, les fonctions de
1-1 déplacement (6) assureront automatiquement, aussi bien
la continuité des déplacements ws, vs et des pentes y/, <p
X-0 X X=A entre les macro-éléments de la même poutre, que la compa¬
Fig. 6. Système d'axes « global » et système « local ». tibilité des déplacements w entre les poutres de l'ouvrage.

172
Notons toutefois (voir (5d)) qu'on ne peut pas assurer auto- Za Zr
dvg
matiquement la continuité des pentes —— entre les élé-
dx
ments dès que le pont n'est plus rectangulaire (tg a ^zé 0).
La conséquence de cette imperfection, qui d'ailleurs ne.
peut intervenir que si zs ^zé 0, est toutefois minime. Dans
la deuxième partie (problème 5), nous démontrerons que A,° 90° Z. "B
l'interaction des poutres n'est que très peu influencée par
la position zs des charnières et par conséquent aussi par A2
la flexion horizontale des poutres, de sorte qu'il n'est pas
nécessaire de rechercher une très grande précision.

É?B
5. Matrice de rigidité du macro-élément
Soient S le vecteur des forces nodales généralisées
(figure 7) et w le vecteur des déplacements nodaux asso¬
ciés d'un élément :

Wai
MA <Pax

wBl
Si MB <Pbx
(7a, b)
Sa" W2 WH Fig. 7. — Elément fini et ses efforts nodaux.
MA (PA2

MR En utilisant les équations (5), (6) et (7b)

On cherche la relation
S kw (8)

où k est la matrice de rigidité d'ordre 8x8.

Soit Ut l'énergie potentielle de déformation de l'élément
correspondant aux déformations données par les équations
(4), (5) et (6). Soit Ue l'énergie potentielle des efforts
externes concentrés aux angles du même élément et repré¬
sentés par le vecteur des forces nodales. La somme
-y,"

u=u(+ue B"
représente l'énergie potentielle totale de l'élément. Si l'élé¬
ment est en équilibre, on sait que, selon le théorème du
minimum de l'énergie potentielle, la variation d'énergie
1,1
-ÌB'-ItgaB» -B'--tgaB"
o
®-
2,
potentielle totale correspondant à une variation infinitési¬
male de ses déformations compatible avec les liaisons est
nulle
--B--tgaB -B--tgaB
ôU ô (U( + Ue 0

Pour une série de fonctions de déplacement approchées

V+ JtgaB*
b 2
-ÌB"+^tgaB"
o Z

(voir (4)-(6)) et caractérisées par un vecteur w de dépla¬

cements nodaux, la meilleure solution sera celle qui donne Cw, (H)
le minimum de U par rapport à ces déplacements consi¬
dérés comme variables. On doit, par conséquent, poser de sorte que
A A
du l- /Vm dx Ì wr Tc^D Cdiê-w
dvt dvi (Uf+Ue) 0, (9)

où la dérivation partielle s'opère par rapport au vecteur

défini par (7b). (l'indice T signifie que la matrice doit être transposée).
Pour exprimer l'énergie potentielle de déformation (voir L'énergie potentielle des efforts extérieurs étant
(2), où selon les hypothèses adoptées on néglige l'effet de
N), on introduira le vecteur U.- -wTS,
My EL, w", l'énergie potentielle totale vaut

m
Ts GJ ¥' De. (10)
A

Mz EL £/= Ut+ tf.-iw* fi iCdi-w-wTS.

B EC y," 0

173
Par dérivation (voir (9)), On cherche maintenant la relation

du P Kr, (15)
dv/
CrDCd*w-S 0,
où K est la matrice de rigidité de l'assemblage ; on l'ob¬
tiendra en assemblant les matrices de rigidité (k)} k des
d'où
éléments. A cet effet, on écrira pour un élément Q) quel¬
conque (voir la figure 5) la relation (8) sous la forme
S= / CTDCd3cw
[sf
1
kii k12 k13 k14" r<

Compte tenu de (8), la matrice de rigidité de l'élément Si+i k22 k2s k24 ri+i
est donc (S), (k), (w). Wm k33 k34 ,(16)
A
Symétrie
k= /CrDCd3c. (12) Sfc+i } *44_IJ r*+i
où les éléments ku, ki2 de la matrice de rigidité (k)^
La matrice de rigidité k obtenue selon (12) est donnée sont des sous-matrices d'ordre 2x2. Pour l'assemblage
dans les tableaux 1 et 2. Pour des raisons pratiques, on d'éléments, la relation (15) s'écrira par analogie :
l'exprime sous la forme
k k + Ak (13) Pi Kn Ki2- "Kln
où k (tableau 1) est la matrice fondamentale représentant P2 K22_ "Kg»
uniquement l'effet de la flexion verticale et de la torsion de I I

Saint-Venant et Ah. (tableau 2) est la matrice complémen¬ (17)

taire représentant la correction due à l'effet de la flexion
Symétrie
horizontale et de la torsion fléchie.
K„

6. Matrice de rigidité de l'assemblage où les éléments Ku, K12 de la matrice' de rigidité de

l'assemblage sont de nouveau des sous-matrices d'ordre
Une fois connue la matrice de rigidité d'un macro¬
élément, l'établissement de la matrice de rigidité de l'as¬
2x2. On obtiendra la matrice K en considérant à tour de
rôle la contribution Kj de chaque élément Q). Ainsi la
semblage n'est plus qu'une question de routine. Toutefois,
contribution de l'élément (J) de la figure 5 est par exemple
pour que même un lecteur moins au courant de la méthode
des éléments finis puisse suivre le texte sans difficulté, il
parait opportun de décrire brièvement le procédé. 1 2 i M-l k k+\ n
Avant d'établir la matrice de rigidité de l'assemblage 0"
d'éléments finis, il faut passer (fig. 5) de la notation locale
o
o-p--o 0 0 0 1

0 0 0 0 0 0 2
des nœuds (Ai, Blt A%, B^) à la notation globale (1,2
n) et numéroter également les éléments (©, @... (in)). En
remplaçant pour chaque élément la notation locale par la no¬
\\ 1

1
1

1
1

1
1

1
1

1
1

tation globale (c'est-à-dire en remplaçant par exemple pour

11 1 1 1 1 1

l'élément (J) de la figure 5 les déplacements locaux wA kn ki2 k13 kM 0

k22 k2S k24 0 /+
9a\\ wb, • ¦ ¦ P31 les déplacements wt, tpt, wi + 1
l'assemblage) on établit automatiquement la continuité
de
K, i\ 1
1

1
1 1

1
1

entre les éléments.

Définissons sur l'assemblage deux vecteurs
\l k83 ks4
1 1

0
1

k
Symétrie lt44 0 k-i
Pi 1 1 1

Pa
et (14a, b)
\ 11 1 1

0 n

et la matrice de rigidité totale est par conséquent

qui soient respectivement le vecteur des forces nodales
et le vecteur des déplacements nodaux. Leurs éléments
sont les sous-vecteurs K=^K,. (18)
/=!
Pi
Pi et (14c, d) Notons que cette façon de former la matrice K revient à
M,
satisfaire à chaque nœud de l'assemblage la condition
où Pi et Mi sont les forces nodales équivalentes. Elles statique
remplacent dans le calcili les véritables forces extérieures,
qui peuvent être par exemple réparties. Les nœuds d'ap¬
pui mis à part, on peut les obtenir au moyen des fonc¬
tions de déplacement (4) et (6) en exprimant l'égalité de
2 ^
où la somme concerne les quatre éléments autour du
leur travail virtuel et celui de la charge donnée. nœud /.

174
Tableau 1

Matrice de rigidité fondamentale — Effet de la flexion verticale et de la torsion pure

01+Û8 03+07 —a!—o3 Û5+Û7 ai—as 05—07+09 —01+03 05—07—09

a2—û4+2û8 —af,—ai —ag—as 05—07 —a9 —o2+2a8 —05+07+09 —fl4+a6+a8

ai+a3 —a&—a7 —ai+a3 —05+07—09 ai—a8 —05+07+09

a2+at+2a$ 05—07+09 04+06 + 08 —05+07—09 —a2+2o8

Ol+«8 05+07 —ai—as 05+07

02+04+208 —05—07 —06+08

Symétrie

01+03 —05—07

a2—04+2o8

GJ 1 G/ 1
G/
tg2„ a
al
^' 2A
+ ~^r
GJ
£/,
a4
26iz;tga a?
10 62 ET,
1
/ f77"
o2 .05 2^V1 + ^tg2a a8 —Â
2^
1
+ Ü!
£/„
tg2 a

6 G/ zt G/ 1 G/
6J £Ä tga
03 06 09
5 62J EL 30 62I/i,

Tableau 2
Matrice de rigidité complémentaire — Effet de la flexion horizontale et de la torsion fléchie

61+63 65—67+69 —61—63 65+67+69 -61+63 -65+67+69 61—63 —65—67+69

62+64—268 —65+67—69 64+65 —65—67+69 —62+64 65+67—69 64—66—268

61+63 —65—67—69 61—63 65—67—69 -61+63 65+67—69

62+64+268 -65+67+69 64—6g+26s 65—67—69 —62+64

E(Itz* + C)
61+63 65+67+69 —61—63 65—67+69

62+64+268 —65—67—69 64+65

Symétrie

61+63 —65+67—69

62+64—268

12
61
62J3
6' ÎZstga
h Wa 68 ^-8tg«
36
u
6s 2
a *6 tg2 a
2B tg
69
P2 X4

175
La matrice de rigidité donnée par (18) est singulière et Pour exprimer les déformations ws, va, y/ et leurs déri¬
comme telle ne possède pas de matrice inverse. Ceci est dû vées par rapport à x on utilisera les relations (5) et (6).
au fait qu'il est impossible de résoudre le système d'équa¬ Pour les forces nodales Z, compte tenu des équations (7)
tions (15) sans prescrire un nombre nnnimum de déplace¬ et (8), on se servira directement des lignes correspondantes
ments empêchant l'assemblage de se déplacer à la façon de de la matrice de rigidité k k + Ak donnée aux tableaux 1
solides rigides, puisque les déplacements wt et les rotations et 2. On exprimera ainsi le vecteur n en fonction des para¬
(Pi ne peuvent pas être uniquement déterminés par les condi¬ mètres nodaux aux angles de l'élément sous la forme
tions d'équilibre. Le système d'équations (15) doit, par
conséquent, être encore modifié conformément aux condi¬ Fw (21)
tions d'appui. Sur un nœud d'appui simple /, on peut poser
par exemple où F est la matrice de coefficients d'ordre 12x8 (voir
tableau 3) et w le vecteur de huit déplacements nodaux
w, wt,
défini par l'équation (7b).
où m>j est soit nul, soit égal au tassement d'appui donné. En En examinant l'équation (20), on se rend compte que le
cas d'encastrement, on aurait de manière analogue vecteur n ne contient pas l'effort tranchant Qy. Son absence
est due au fait que les hypothèses simplificatrices adoptées
au paragraphe 3 rendent impossible la détermination des
Wi
réactions horizontales entre les poutres au niveau des char¬
<Pi
nières de liaison. De façon générale, on a

Les valeurs ci-dessus peuvent être introduites dans le sys¬ dM2

tème (15), tout en éliminant les équations correspondantes
dx
devenues superflues, ce qui revient à réduire l'ordre de la
matrice K et des vecteurs P et r.
où »iz mz(x) représente une charge répartie, constituée
Pour éviter de devoir procéder à la réorganisation com¬
plète de la matrice K et des vecteurs P et r dans la mémoire par des couples horizontaux [tm/m] qui s'obtiennent
de l'ordinateur, la pratique courante consiste à multiplier
comme une des composantes de l'effet de la transmission
des interactions des poutres sur l'axe de torsion. Notons
sur la diagonale de K le coefficient correspondant à wt
que les valeurs de Qy sont très faibles et pratiquement
ou Çi par un très grand nombre (mettons 1012) tout en
sans importance.
remplaçant en (15) l'élément correspondant du vecteur P
En appliquant à tour de rôle la relation (21) aux éléments
par le même grand nombre multiplié par wt ou 7pi prescrit. d'une poutre du pont, on constate, en général, une légère
Cette opération a pour conséquence de remplacer dans (15)
oscillation des efforts obtenus respectivement à l'extrémité
l'équation superflue par une autre exprimant pratiquement droite de l'élément Q) et à l'extrémité gauche de l'élément
l'égalité wi vpj ou <pi ç?j. La résolution du système
(y+1) (et ceci même en l'absence des charges isolées qui
d'équations ainsi modifié
pourraient justifier une certaine discontinuité). Cette petite
K- (19)
imperfection est propre à la méthode habituelle des élé¬
ments finis, et pratiquement peu importante. Notons
est le vecteur r complet des déplacements nodaux, y com¬ qu'une meilleure approximation serait obtenue en prenant
les moyennes arithmétiques des valeurs de gauche et de
pris les déplacements et (ou) les rotations prescrits aux
droite.
appuis.

r-P =25 dkg

7. Calcul des efforts internes c

1 î
Soit n le vecteur des efforts internes aux extrémités
i
0 et x A de l'élément :
i
~x

-ElyWl
i
MyA f (0)
MzA -EIzv"e (0) i Y
BA -EC y/" (0)
15 cm 15 cm
GJ y/' (0) EIy GJ
*SA
TfA -EC y/" (0) Il 2 3 4
,-12,5 dkg
(\ 5 6 7 e
QzA ~Zai~Za» b) no 11 12 ,3 m IS 16 17
(20)
19 20 21 22 23 24 25 26
MyB -EIyw"a (A)
MZB
BB
-EIzv"B (A)
-EC y/" (A)
k-
i
29 30 31 32 33 34 35

mai 38 39 40 41 42 43 44
GJ y/' (A)
I

TgB
NÉ 47 48 49 SO SI 52 53
TfB -EC y/m (A)
QzB Zb1 + Zb3 8x3,75 cm

où T, et Tf sont respectivement les parties du moment de Fig. 8. — a) Modèle du pont rectangulaire de l'essai de
Kopecky [14].
torsion T dues à la torsion de Saint-Venant et à la torsion
b) Assemblage d'éléments finis (M 8) pour l'ana¬
fléchie. lyse numérique du même ouvrage.

176
Deuxième partie: Analyse numérique et applica¬
tions
Pour faciliter le calcul pratique des ponts biais à poutres n'entrerons pas dans les détails, puisque ce programme,
multiples sans entretoises, nous avons élaboré un pro¬ qui sera mis plus tard à la disposition des clients du centre
gramme standard en langage FORTRAN pour l'ordina¬ de calcul de l'EPFL, fera l'objet d'une publication ulté¬
teur (programme PBPMSE). Ce programme, basé sur la rieure.
théorie développée dans la première partie de ce travail, Dans le texte qui suit, on donne les résultats de quelques-
calcule pour le pont donné, avec les charges et les tasse¬ uns des nombreux problèmes résolus à l'aide du programme
ments d'appui choisis, les déplacements (w, <p) et les PBPMSE. Les trois premiers ont été choisis pour vérifier
efforts internes (My, Mz, B, Ts, Tf, Qz) dans les poutres la théorie, en confrontant nos résultats avec ceux de la
de l'ouvrage. Le nombre d'éléments par poutre M l: A littérature connue. Les autres exemples nous ont permis
peut être choisi en fonction de la précision voulue. Nous d'étudier quelques questions intéressant la pratique.

Tableau 3

Matrice F pour le calcul des efforts internes aux extrémités des éléments

Cl 2C2 —Cl C2 Cl 2c2 —Cl C2

—zs fois la ligne 3

-rfi+fls —2a\4-di di—d3 —dz+di 01+03 M2+04 —di—da 0*2+04

«i 2c2—«a —ex ez Cl 2e2+es —ci C2

h h -h h -h -h h -h
—2fl5+09 —205—09 —2an—OQ —205+09
—lai— 26a 2oi+263 —2aj.—2b3 201+263
+267-269 -267-269 +267-269

—Cl —ca Cl —2C2 —ci —c2 Cl -2c2

—z, fois la ligne 9

di+d3 02+04 —0*1— as 2rf2+04 -dì+03 —0*2+04 di—ds -202+04

—ei —e2 «1 —2«2—e$ —ci —e% Cl —2ea+es 10

comme la ligne 5 11

comme la ligne 6 12

Coefficients «1, 05, 09, voir tableau 1

Coefficients 63, 67, 6g, voir tableau 2

Eh,
ci-3EI» ca

6 Eh 2 Eh 6 Eh
rfi bA1
dz
bA
da
A» tga

3GJ GJ GJ
ci -2â-tga ea
-j tg a
C3=T
12 EC 6 EC
A-- bA» h ~bW

177
MOMENTS FLECHISSANTS My [kgcm] MOMENTS DE TORSIONI [kgcm] Problème 1 : Confrontation des calculs avec l'expé¬
Rjutre I rience
Pour la première vérification de la théorie, nous avons
choisi le modèle d'un pont rectangulaire à cinq poutres
Poutre I
(figure 8a) qui a fait l'objet d'une étude expérimentale de
J. Kopecky [14]. Le modèle était réalisé en résine polyester
durcie (extroplex) et les charnières étaient représentées par
des encoches minces, fraisées dans la plaque compacte du
matériau.
Pour les dimensions de la section transversale des poutres
et les qualités physiques du matériau (E, G), la réfé¬
rence [14] donne EIy GJ. Les charnières étant à mi-
hauteur des sections, on a zs 0 et les inerties Iz n'inter¬
viennent pas (voir le tableau 2). De plus, si l'on néglige
pour les sections pleines le petit effet de la torsion non
uniforme (C 0), la matrice de rigidité supplémen¬
taire Ak est nulle. Ainsi, pour le cas envisagé, la matrice
de rigidité k des éléments sera représentée uniquement par
la matrice fondamentale k du tableau 1.
Le calcul du modèle a été exécuté pour deux réseaux
différents d'éléments finis : celui de la figure 8b (M 8),
et le réseau deux fois plus fin (M =16). Excepté dans la
région voisine de la charge concentrée P (fig. 8a), les deux
calculs conduisent à des résultats peu différents. La
figure 9 montre que les résultats du second calcul (M =16)
concordent bien avec les valeurs expérimentales.
On remarquera, sur la figure 9, que sous la charge
concentrée P, le diagramme des moments fléchissants pré¬
sente une pointe très aiguë qui, probablement, influence
considérablement la répartition des moments entre les
Fig. 9. — Confrontation des résultats du calcul du pont de la
figure 8 avec les résultats expérimentaux obtenus sur le modèle poutres du pont. Ce fait devrait être pris en considération
(valeurs entre parenthèses et diagrammes en traitillé). dans le calcul d'un pont réel, où de vraies charges concen¬
trées n'existent pas. Remplacer une charge répartie par des
charges isolées pourrait conduire à des résultats assez
différents.

100 cm
a) KL 61 -18
94 Problème 2 Confrontation avec d'autres méthodes
de calcul

20 Iy 3,444. lO'V Comme seconde vérification de notre méthode, nous

GsS
65 85 avons choisi le pont biais de la figure lia, formé de sept
J 5.900 .10'2m'1 poutres-caissons préfabriquées de type KL 61-18 (fig. 10a),
soumis à une charge uniformément répartie p 1 t/m2.

70
)-.
98cm

b» 2x0,375"
U)
\tâk\ I
n
m—m
||
-

i
M
"itt* BIV-36 ßu jN N 3.
\l\ s
»

3=S
Sl.
f" Iy 0,1587.10s In"

Iz «0,1163.106i.n*
\\ m
m
19m
\[
\\
1 2 3 r. 5 6 7 6
br \io \l \« \B \« \l5 \« \l7 \lB v9

J
C
0,1979.10'in*
0,3114 -10sLn*
\»9 \20 \îl
\28 \so
\,23 \x \?5 \26 \27
\j9 \ 2

\ai \s2 \33 \34 \a5 \S6

31 \3B \40 \39 \*2 \43 \45 \tt
\ \
55 \41
\*6 \47 » \49 S0 \51 \5Z \S3 \s
\55 \56 \57 \56 \59 \60 \61 \62 ^3
\et
nés \66 \67 \68 \69 \70 \7t V2
36 in

Fig. 10. — Quelques exemples réels de poutres-caissons préfa¬

V- 8x 2.375m
\
briquées. Fig. 11. — a) Pont biais formé de sept poutres-caissons
KL61-I8.
a) KL 61-18 (Tchécoslovaquie) pour les portées jusqu'à 20 m; b) Assemblage d'éléments finis (M 8) pour l'ana¬
b) B IV-36 (Etats-Unis) pour les portées jusqu'à 88 ft. lyse numérique de l'ouvrage.

178
Tableau 4

Pont biais constitué par sept poutres-caissons KL 61-18

Confrontation des résultats de différentes méthodes de calcul pour la charge répartie p 1 t/mz

Flèches w [cm] au milieu des poutres

I II III IV V VI VII
Méthn\.

Kop. [16] 1,51 1,55 1,57 1,58 1,57 1,55 1,51

El. finis 1,509 1,548 1,568 1,574 1,568 1,548 1,509

Moments fléchissants My [tm] aux quarts de la portée des poutres

X 0 0.25 0,50 0.75 1,00

1

R Jir. [9] Kop. [16] El. finis Jir. [9] Kop. [16] El. finis Jir. [9] Kop. [16] El. finis Jir. [9] Kop. [16] El. finis Jir. [9] Kop. [16] El. finis

ïî^ -7,08 -8,48 -8,45 15,37 15,59 15,23 31,41 31,64 31,55 29,57 29,02 29,56 -4,91 -7,67 -4,28
II -5,37 -6,30 -5,78 18,32 18,10 18,25 32,53 32,30 32,69 27,07 26,77 27,19 -2,94 -4,79 -2,73
m -3,81 -5,00 -3,89 20,59 20,44 20,60 33,13 32,96 33,29 24,52 24,39 24,84 -2,30 -4,05 -2,06
IV -2,76 -4,23 -2,96 22,50 22,39 22,68 33,33 33,20 33,48 22,50 22,39 22,68 -2,76 -4,23 -2,96
v -2,30 -4,05 -2,06 24,52 24,39 24,84 33,13 32,96 33,29 20,59 20,44 20,60 -3,81 -5,00 -3,89
VI -2,94 -4,79 -2,73 27,07 26,77 27,19 32,53 32,30 32,69 18,32 18,10 18,25 -5,37 -6,30 -5,78
vu -4,91 -7,67 -4,28 29,57 29,02 29,56 31,41 31,64 31,55 15,37 15,59 15,23 -7,08 -8,48 -8,45

z -29,17 -40,52 -30,15 157,94 156,70 158,35 227,47 227,00 228,54 157,94 156,70 158,35 -29,17 -40,52 -30,15
Ce problème a déjà été résolu par d'autres méthodes appro¬ la contribution de l'effort normal est
chées (voir [9] et [16]), dont les résultats sont confrontés
avec les nôtres dans le tableau 4 (pour M 16). On Nm 13/>!
remarquera que nos résultats concordent, en général, très
&N max —
710
0,018 ft
bien avec les deux autres solutions approchées.
Notons que les deux méthodes [9] et [16] considèrent, ce qui ne représente que 0,85 % de o~m Ajoutons
à priori, que les charnières entre les poutres se trouvent au que les résultats de notre méthode (déformations, efforts
niveau des centres de rotation S (ce qui est approximative¬ internes) concordent très bien avec les résultats du calcul
ment vrai pour les poutres-caissons KL 61-18) et négligent exact. Ainsi, les hypothèses simplificatrices adoptées se
le petit effet de la torsion fléchie. Dans l'application de trouvent pleinement justifiées.
notre méthode, qui est plus générale, nous avons par Notons encore au sujet des contraintes normales que,
conséquent posé zs 0 et C 0. On verra plus loin en plus de celles dues aux moments My et aux efforts nor¬
l'influence de ces approximations. maux N, d'autres contraintes normales viendront s'ajouter,
qui seront dues aux moments Mz et aux bimoments B. Par
notre méthode, on trouve (toujours pour la même section
et le même cas de charge)
Problème 3: Vérification des hypothèses simplifi¬
catrices 540/>!
Les deux problèmes précédents concernaient des ponts
°m. 116 300
18 0,084 a
avec charnières au niveau des centres de rotation. Dans
ces cas, l'influence des hypothèses simplificatrices, énon¬ Avec B max 45 fà et co 36 in2 (où | co
cées au paragraphe 3 de la première partie de ce travail,
| | | \

est la valeur absolue de la coordonnée sectorielle aux

„^ |

„^
n'apparaissait pas, puisque les poutres ne fléchissaient que angles de la section)
dans le plan vertical et le fait que leurs sections étaient
pratiquement symétriques par rapport au plan z — z8 a h "
¦—! max i
w
i

max
4J
1— 36 £
0,005»i
imposait des déplacements longitudinaux u automati¬
quement nuls au niveau des charnières.
Pour vérifier la validité de ces hypothèses simplifica¬ Les contraintes a dues aux bimoments B sont donc prati¬
trices, nous avons étudié le pont rectangulaire de la quement sans importance.
figure 12a, constitué par dix poutres-caissons américaines
de type B IV-36. Selon la figure 10b, la position théorique
des charnières sitrouve à peu prés à douze pouces au-
dessus des centres de rotation. Les inerties
/;,, Iz et le Problème 4: Influence de la position des charnières
moment d'inertie sectorielle C sont donnés dans la même
Dans le cas d'un véritable pont à poutres-caissons pré¬
figure. Etant donné que la conversion des unités améri¬
caines en cm (1"
fabriquées et assemblées par le remplissage des joints
1 in 2,54001 cm) donnerait des chif-
(« shear key »), il régnera probablement toujours une cer¬
fres peu commodes pom: le calcul, nous avons conservé les
taine imprécision quant à la position exacte des charnières
unités originales, les pouces.
Comme charge du pont, nous avons considéré la charge par lesquelles on remplace approximativement l'effet d'une
telle jonction, n est intéressant de savoir dans quelle
linéaire p Pi sin n - appliquée successivement : mesure la position des charnières influence l'interaction des
poutres.
a) sur la bordure extérieure (indice 0) de la poutre Nous avons examiné les deux ponts (rectangulaire et
I; biais) de la figure 12 et avons exécuté le calcul de chacun
b) au milieu de la largeur de la poutre I. d'eux pour deux positions zs différentes des charnières :

La raison de ce choix est que pour une charge sinusoïdale,

on trouve relativement facilement la solution exacte du
problème (voir [12]). Les équations sont applicables uni¬
quement pour des ponts rectangulaires à poutres de même
section. Pour faciliter leur application, nous les avons pro¬
grammées. 10x36'
Le tableau 5 résume les résultats des deux calculs, la
solution exacte et le calcul approché par notre méthode
(M — 16). En accord avec les hypothèses simplificatrices
du paragraphe 3, les résultats du calcul exact confirment
que les déplacements latéraux v des charnières sont vrai¬
ment très petits vis-à-vis des déplacements verticaux w.
Sous premier cas de charge, ce rapport est par exemple
pra 10888 pour la première poutre, et 19: 4477 pour la
dernière poutre du pont. D'autre part, le calcul des con¬
traintes normales confirme la très faible influence des efforts
normaux. Ainsi, par exemple, alors que dans la poutre I
le moment fléchissant My
max engendre la contrainte 1000

Fig. 12. — Pont rectangulaire et pont biais formés de dix

Cm. max
y
My max n 15 950
a 21 2,1 Pi poutres-caissons B IV-36. Dimensions en pouces (1* 1
pouce
158 700 2,54 cm).

180
BLEA

•< a

ft, 3 I •S •h for
H 5 tsä •8.X g o res-
&5 S<3 O1 «»

i
érie
u X sur
1O £>o XI x T3o
out
p. o

i
8g s.*
s oo|i S^ift*' •—'
(S
ô£ So o^ ô' "©^
•bÖ
K
oëjl sa""•/¦-S
h.
Exact Exact
«
finis finis finis inis
£1. £1.
UJ
x W 09
txJ act
m wX m en tu W mX u w §
CQ
*

O^NW^^^hOOOsO 00\hO\<-iVO-Mh-00
OOvoot^^o>n>n^-Tf'^tJTf
ooTfr<ic^irioofno\vOl^'st
OONOOt^vOlTiiri^-Tj-TfTt
oo oow-tr-vooOTtON^-iQOr-
_,hhï;>£Î>£8x
<nmNOO\oooo
ONOV<NOO^r*OWiWO Vi^^MMMmi
S \Ocncsoo\oooor-r--t**
vor-omr-r^i^mooTt"
O^fSoor^oo-HinNO
«n^Tj-Tj-ONOor-^oNin
ï-H "/ioono*oof>oN^"
Ttmoo-HinQifiO'ON
— «O'TrOMfiNNNNtS
Tj-iftooh-'Hhri^'OO
fS«nr^r*'-tr-(N'*oo
riocor-^TtcnrNM'-i
<OMfnOc7\CiOO\ON^ ~n'OOn'n^'nMO^
in^-nnn(S(NiNNM
t-Hrnoo^-or^^o^rnr^
i-«T*'00,*Ot--<o,«tmcNl 1—1

tf T—« m oooooomi-iO
m »— »s
n »—I
n 1

i/î
oo ON
o
1
o
1

m 1

on
^o
r^ oo^h- 1

vî
t^r^
tj-

m Il ti.

u Exact
^3ï
II

finis
><

Exact
I
S,w
s
finis
ft:
*s»
oO SjRen o3 I So

gj
h,"
oO £>S Et
K
O fe;» sa
**•/•**

El. El.
wx M va m Wx ¦g w a x W M u
m M '3 Wr| u U; <n
c

Hirifnn^oomwhift'n
i-«vif»c»»nÄrncKt*,,ftv>
©*-<C4tn',<tW»\Ol>00O\© O^OdON-1
00Tfiri\o>riCM»0N-"^O N^nnrlNrHrt"nONHt0\^'n'-'O
_,
a
¦*wrtOVimNO\0
"TOtna>^ootH\orlO *~~
0«>ûTt,os^o|nifiosN
'O'VMO^OOOOt^t^t^ vi'OooHinQ^ovOM
csl'*t-^-for^c^*-i'-t ,*oor^\o,^'cncsci^^
r^ooorr>t*-(SO<riooTi"
O
v~> ÛOO\tt 00 O\^O\^,©00V.'^tmfN
TtrJ-cncnc^rNfNtNfNCN
,"*OcS'-*T*av,*»n»-<Q *""0\0i-m<o«n>nmom
OTt-MTj-OOOVi^tnN
0\^,rHO\00
3 rn
£; g On ^ONOOONOOl'd;

o
> r- 1
o
ON
£ co
mNOn > *—< 0\ON^' 1

£ vo 1

!$ <N
co
x 0\ »-<
ti

181
a) zs — 12 in (charnières près de la surface supérieure Problème S : Influence du biais sur l'interaction des
des poutres d'après la figure 10b) ; poutres
b) zs 0 (charnières au niveau du centre de torsion Une opinion assez répandue en pratique admet que
des poutres). l'effet du biais sur l'interaction des poutres du pont est
négligeable. Ceci dit, on applique les mêmes coefficients de
La poutre I a été soumise à un groupe de trois charges répartition transversale pour un pont biais que pour un
concentrées : Pj 0,25 Q, P2 0,50 Q et P3 0,25 Q ouvrage rectangulaire. Pour vérifier dans quelle mesure se
appliquées respectivement aux 7/i6, 8/16 et 9/16 de la portée. justifie cette simplification, nous avons considéré les deux
(Leur effet correspond approximativement à une charge Q ponts de la figure 12a, b qui correspondent respectivement
uniformément répartie sur un tronçon central de 1/8 de la aux angles d'obliquité a 0° et a 40°.
portée. Notons toutefois que les charges nodales équiva¬ Dans un premier cas, nous avons supposé une charge
lentes à une telle charge répartie, au sens des travaux vir¬ uniforme p appliquée simultanément sur toutes les poutres
tuels, comprendraient encore des couples nodaux de forces. du pont. Dans un pont rectangulaire, elle provoque dans
Mais la différence pratique n'est pas très grande.) toutes les poutres les mêmes moments fléchissants avec
Les résultats du calcul pour les deux positions des char¬ maximum Mymax 0,125 pP. Les moments de tor¬
nières sont résumés dans le tableau 6. sion T sont nuls. Dans un pont biais, les moments de
| Le calcul a été exé¬
flexion My engendrés par la même charge sont plus petits,
mais il apparaît des moments de torsion T. En analysant
cuté avec les éléments finis de longueur A — I Pour le les diagrammes de My de la figure 13, où nous avons porté
16/
niveau inférieur des charnières, les flèches de la poutre
chargée et des poutres voisines accusent une légère aug¬
les résultats de notre calcul I pour M — 16

mentation, tandis que pour les poutres situées près de tate que les valeurs maxima des moments fléchissants ont
l'autre bord du pont, les flèches ont une légère tendance à diminué d'environ 21 %.
diminuer. Les différences ne sont pas toujours apparentes En pratique, la charge uniforme qui pourra être prise en
avec le nombre de décimales choisi. En ce qui concerne les considération pour la répartition transversale des effets dus
efforts internes, on constate sur la poutre chargée (respec¬ à l'interaction des poutres sera le plus souvent le poids
tivement sur la poutre voisine) une légère augmentation des propre du tablier. Par contre, le poids propre des poutres
moments fléchissants, ceci au profit des poutres plus éloi¬ sera porté indépendamment par chacune d'elles, étant
gnées où la différence reste toutefois trop petite
pour être donné que le remplissage des joints ne s'effectue qu'après
mise en évidence avec le nombre de décimales utilisé. Les la mise en place des poutres. Notons toutefois que chaque
moments de torsion T restent pratiquement inchangés. Il poutre isolée bénéficiera de l'effet favorable du biais, qui
est, par contre, naturel que la position des charnières se manifestera par une certaine diminution des moments
influence fortement les moments de flexion Mz. Mais, fléchissants (à peu près 18 % dans le cas envisagé). Ceci
comme on l'a déjà vu dans le problème 3, vis-à-vis des est dû au fait que l'effet du biais des appuis équivaut,
moments Mv, leur contribution aux contraintes normales a dans une certaine mesure, à un encastrement partiel. Pour
n'est pas très importante. ce qui est des surcharges, notons, par exemple, que les

-0,025

D.125 0.250 1375 0 50C 0.625 0.750 0.875 I.OOC

Tmax

F>ta
0,025
0,0285
0.0321
D347

363
037

My Pi

Fig. 13. — Diagrammes des moments fléchissants Mu et les valeurs maxima

des moments de torsion T pour les cinq premières
pou très du pont biais de la
figure 12 b. Poutres-caissons B IV-36, charge uniformément répartie p par
unité de longueur de chaque poutre. En traitillé : effet de la même charge sur
un pont rectangulaire.

182
Tablbau 6

Influence de la position zs des charnières sur l'interaction des poutres. Pont rectangulaire
I
et pont biais de la figure 12. Poutre chargée par un groupe de trois forces concentrées :
P XU Q, P2 V« Q et P3 1UQ, appliquées respectivement à 7/16, 8/16 et|»i6 de
la portée

Pont rectangulaire (« 0)

10*El warn : Q Muum : Ql Mutm : Ql T™* :Ql

R
z, 12" z, =0 z. 12' z. 0 z» 12" z. 0 z, 12' z. =0
I 0,2014 0,2016 0,0442 0,0444 0,0015 0 0,0080 0,0080
II 0,1815 0,1816 0,0365 0,0365 0,0023 0 0,0149 0,0149
m 0,1583 0,1584 0,0290 0,0289 0,0014 0 0,0128 0,0128
IV 0,1395 0,1395 0,0242 0,0242 0,0010 0 0,0109 0,0109
V 0,1244 0,1244 0,0209 0,0209 0,0007 0 0,0090 0,0090
VI 0,1124 0,1123 0,0185 0,0184 0,0005 0 0,0072 0,0072
VU 0,1032 0,1031 0,0167 0,0167 0,0004 0 0,0055 0,0055
VIII 0,0964 0,0964 0,0155 0,0155 0,0002 0 0,0039 0,0039
IX 0,0921 0,0920 0,0147 0,0147 0,0001 0 0,0023 0,0023
X 0,0899 0,0898 0,0143 0,0143 0,0000 0 0,0008 0,0008

Pont biais (a 40°)

WElw am Q '. Muai i:Ql Mzan ï-.QI Tm» ¦ Ql

R
z. 12' z, =0 z. 12' z. 0 z> 12* z. =0 z. | 12* z. 0

I 0,1802 0,1804 0,0408 0,0410 0,0012 0 0,0146 0,0146

II 0,1576 0,1577 0,0325 0,0326 0,0021 0 0,0200 0,0200
m 0,1309 0,1310 0,0241 0,0241 0,0010 0 0,0165 0,0166
IV 0,1086 0,1087 0,0187 0,0188 0,0007 0 0,0134 0,0135
v 0,0905 0,0905 0,0150 0,0150 0,0004 0 0,0108 0,0108
VI 0,0760 0,0760 0,0123 0,0123 0,0003 0 0,0086 0,0087
VII 0,0646 0,0646 0,0102 0,0102 0,0002 0 0,0067 0,0067
vm 0,0559 0,0559 0,0087 0,0087 0,0001 0 0,0052 0,0052
IX 0,0494 0,0494 0,0075 0,0075 0,0001 0 0,0038 0,0038
X 0,0446 0,0446 0,0066 0,0066 0,0000 0 0,0026 0,0026

normes françaises pour les ponts routiers prévoient comme si la charge Q se déplace le long du bord extérieur (indice 0)
surcharge du type A une charge uniformément répartie. de la poutre I, le moment My
m^ max apparaît quand la
Vis-à-vis d'un ouvrage rectangulaire, celle-ci bénéficierait charge se trouve à peu près à Vie / à gauche du centre. Le
donc pleinement de la réduction des moments de flexion My moment My max max ainsi obtenu vaut 488.10~4 Ql (contre
(21 % dans notre cas). 484.10-4 Ql pour la charge au milieu de portée). La diffé¬
Un second cas de charge a été envisagé sous la forme de rence n'est donc pas très importante.
charges Q constituées par des charges linéaires p Q : Al
réparties uniformément sur les petits tronçons Al= l: 16.
Nous les avons appliquées à tour de rôle au milieu des
Problème 6: Effet du tassement des appuis
joints r 0,l,2,3,4et5. Notons que par la combinaison
de telles charges, nous pouvons exprimer par exemple Comme dernière application, nous avons étudié l'effet
l'effet des roues d'un véhicule roulant sur le pont. Les du tassement d'appui sur le pont biais de la figure 12b.
résultats du calcul, qui a été exécuté avec des éléments Nous avons supposé que l'angle gauche de l'extrémité
finis de longueur A — 1:16, sont résumés dans le tableau 7. gauche du pont (x 0, y 0) a subi une dénivellation
En plus des moments fléchissants au milieu de la portée, égale à Z/1000. Ainsi, l'extrémité gauche du bord extérieur
nous avons confronté les valeurs approchées des moments de la poutre I est descendue de / : 1000, tandis que l'extré¬
maxima Mv max et les maxima des valeurs absolues des mité gauche du bord intérieur (attachée par le joint 1 à la
moments de torsion 7^ max- On constate que, vis-à-vis
I
poutre H) est restée en place. Pour obtenir l'effet de ce
du pont rectangulaire, les moments fléchissants du pont tassement, il suffit de substituer au déplacement nodal du
biais sont beaucoup plus petits et les moments de torsion nœud 1 la valeur Wi /: 1000.
plus grands. Le calcul a été exécuté avec des éléments finis de lon¬
Notons encore que sur un pont biais, la position longi¬ gueur A — l : 16. Les résultats (moments de flexion Mv et
tudinale la plus défavorable d'une charge Q (pour engen¬ moments de torsion T) sont représentés dans les figures 14
drer les plus grands moments possibles My max max) ne et 15. Quant aux moments fléchissants, on remarquera que
se trouve, en général, pas au milieu de la portée. En procé¬ le tassement affecte essentiellement la région autour du
dant par tâtonnements, nous avons trouvé par exemple que, nœud 1, tandis qu'ailleurs ces moments restent relative-

183
Tableau 7

Influence du biais sur l'interaction des poutres

Dix poutres-caissons B IV-36 de portée l 1000 in. Charge Q p-Al répartie sur Vie
de la portée et appliquée au centre du joint r 0,1, 2...
104- My : Q! au centre de la potrée

Charge Q appliquée sur le joint r

Poutre r - -0 1 2 3 4 5

a =0" 40° 0° 40° 0° 40° 0° 40° 0° 40° 0° 40°

I 525 484 443 406 326 277 266 212 225 169 197 137
II 362 308 407 371 384 345 285 239 237 187 204 151
m 290 229 302 255 367 328 355 315 264 218 222 172
IV 242 180 249 197 273 227 346 305 340 299 254 208
V 209 144 213 158 228 177 258 212 336 295 334 293
VI 184 117 188 129 198 145 218 168 252 206 334 293
VII 167 98 170 108 178 121 192 140 216 166 254 208
VIII 155 83 157 92 164 104 176 120 194 142 222 172
IX 147 71 149 80 155 91 166 106 182 125 204 151
X 143 62 145 70 151 81 161 94 176 113 197 137

104-Af« Ql

I 525 484 443 406 326 299 266 234 225 193 197 163
II 362 341 407 371 384 345 285 256 237 205 204 171
III 290 263 302 276 367 328 355 315 264 234 222 189
IV 242 204 249 219 273 244 346 305 340 299 254 224
V 209 177 213 182 228 196 258 228 336 295 334 293
VI 184 150 188 154 198 165 218 185 252 222 334 293
VII 167 130 170 133 178 143 192 158 216 183 254 224
VIII 155 116 157 120 164 127 176 141 194 160 222 189
IX 147 106 149 110 155 118 166 129 182 147 204 171
X 143 102 145 106 151 113 161 125 176 140 197 163

10».| r|max:ß/

I 170 240 24 111 15 71 11 59 10 50 9 43

II 149 199 149 202 46 112 34 82 29 70 27 61
III 128 164 129 168 130 175 60 119 48 92 45 81
IV 109 133 109 138 110 145 112 154 69 126 63 102
V 90 107 90 111 92 118 94 127 96 143 81 134
VI 72 85 72 89 74 95 76 104 78 114 81 134
VII 55 69 55 70 56 75 58 82 60 92 63 102
Vili 39 53 39 53 40 58 41 64 43 71 45 81
IX 23 39 23 39 24 43 24 47 26 54 27 61
X 8 25 8 27 8 29 8 33 9 38 9 43

ment faibles. Par contre, les moments de torsion sont cements et des rotations entre les éléments de la même
assez élevés et presque constants sur la totalité de la poutre. Dans le cas d'un pont biais, cette continuité n'est
poutre I. La torsion des autres poutres est peu importante. toutefois pas complète pour les dérivées des déplacements
horizontaux v. Cette petite imperfection n'intervient
d'ailleurs que si les charnières entre les poutres se trouvent
Conclusion au-delà du niveau des centres de rotation des sections.
Le présent travail propose une solution approchée du Etant donné la faible influence du niveau des charnières
calcul des ponts biais à poutres multiples sans entretoises. sur les résultats, son effet est minime.
Il se base sur quelques hypothèses simplificatrices qui sont, Mis à part cette petite imperfection, notre modèle est
en règle générale, bien remplies pour les ponts à poutres de
donc du type « conforme ». Comme l'a montré De
section pleine ou fermée (poutres-caissons). Le calcul par Veubeke1, un tel modèle est plus rigide et son énergie de
éléments finis se fonde sur le théorème du minimum de déformation est plus petite que la valeur exacte. Partant
l'énergie potentielle et sur la méthode des déformations. dans le cas d'un ouvrage rectangulaire (où la continuité
Les fonctions de déplacement choisies pour les éléments est parfaite), notre solution converge vers la solution exacte
finis comprennent tous les modes rigides de déplacement de la construction idéalisée, qui obéit aux hypothèses sim¬
et sont susceptibles de donner toutes les déformations cor¬ plificatrices admises.
respondant à une valeur constante des efforts internes.
Elles assurent aussi automatiquement la continuité des
1 Voir,
par exemple, « Displacement and Equilibrium Models
in the Finite Element Methods. Stress Analysts », J. Wiley and
déplacements entre les poutres et la continuité des dépla¬ Sons Ltd, Londres 1965.

184
10ST: El3

10'My El3

® Ol
«N
en
(^ 00 CO
^*"*--ö_ Ö" 0,1 0,1 Cf cäT^-
I 1 1 1 —

1 i i ' 1

\^o% m p oo J- en i
.^^cr^ O 0,1 0,1 o* O 1

fN" cî IM* (M" fN co

;
1 1 1

\ N
i
CM

1
CD O CM en t»
m
0/2
cf

00 Ö" cf o O -^"""T
+
\ \- - 1 1
| i 1

- N. J D,02\ i
CD*

I
j en CD CD
o 0,2
CM co' t£> ID CO ID en W* '
o cf cf cf o ¦"- f
Ci
O o
® 1
-¦ 1
O
I
1
Ö
1
Ö
1
o"
| ^=~
cf
_ \
en
-
¦

cy
m
o o
c»
o
en
(Ml
CM
o"
to'
1 1
o j o
o" CS" cf Ò"
—
rf
®
1
i CO

o" o
o o cf o" o !Jt-~~8~~
1

1
1
1

- 1 1 Cf
\
CD
O
O
o cf
iD
O
CD
o
cf
£2
cf
-_
(D o1

®
o~ CT 1

en «J- ¦rf Ci 00
j L+———k—
Ö" o o o
cf
t u
in ;
1

o
i 1
J_ - 1
1

o
^
¦rf O
en
o o
oo
o O
o O Ö* o" 0~ O o
o CM
El'O
"

en
Ö" Ö"
QO'O
1
cf
1
m W£M
' o
O" C?

¦ I i
- i 1 -

\ ;

o 0,03
m
o o
en
O
o"
CD
O
cf
in
O
cf
o cf
1

o"
1 Cf cf. ;
°<°9__
CM en -rf 5 (M î
¦
O o f o
i
i— 1—
o* o" Ö"

i
Ö*
\
°o" —*f
1

CM
o
co
o
cf
-rf
o
o"
8
o*
o
o"
CD
O
¦rf'
o
c>

o 3L
\* Ì 1

j i
»-
Q.

® I
I
rf
o in m
o ¦rf en
1 1
1 "

—W
«M
0,03
cf o
cf o
O
o
1
I o

I
lO'O Q.
i O ; k s. :
;

CO
o CM 3. O 1 1
•

cT o"
0,09
10,06
o en en en en
¦
1 1
1 -Jir-1 Q02 zo'o 0,03 o o o o ÊO'O

-Â^r-
1

CS o~ o"

1
Fig. 14. — Effet du tassement wi
1000
/ de l'angle supérieur Fig. 15. — Effet du tassement wi / de l'angle supérieur
gauche du pont biais de la figure 12 6 (dix poutres-caissons gauche du pont biais de la figure 12 6 (dix poutres-caissons
B IV-36 de portée / 1000 in). Diagrammes des moments B IV-36 de portée / 1000 in). Diagrammes des moments de
fléchissants. torsion.

Dans le cas d'un ouvrage biais, le problème de la conver¬ avons de bonnes raisons de penser que ce contrôle confir¬
gence est plus compliqué. On peut montrer que, dans la mera la validité de la méthode proposée.
mesure où l'on diminue la longueur des éléments et où
chaque élément s'approche de l'état de contrainte cons¬
tant, la petite discontinuité des dérivées des déplacements
horizontaux v diminue et disparaît à la limite. Par consé¬ Bibliographie
quent, la solution converge. Notons toutefois que le 1. Arya, A. S., Khachaturian, N., and Siess, C. P. : Lateral
calcul d'un ouvrage biais ne converge vraisemblablement Distribution of Concentrated Loads on Multibeam Highway
que vers des valeurs très proches de la solution exacte, et Bridges. Civil Eng. Studies Struct. Res. Series Rep. 213,
ceci pour la raison suivante : l'adoption de la fonction de Univ. of Illinois (1960).
déplacement sous la forme donnée par l'équation (4) 2. Arya, A. S., Lohtia, R. P. : Distribution of Wheel Loads
on Multibeam Highway Bridges. Indian Cone. Jour., 38,
signifie que dans la bande limitée par les bords longitudi¬ n° 4, p. 127-135 (1964).
naux d'une poutre, le déplacement vertical w est linéaire
le long des droites parallèles aux lignes d'appui. Dans le 3. Arya, A. S., Lohtia, R. P. : Distribution of Wheel Loads
on Multibeam Highway Bridges. Indian Cone. Jour., 39,
cas d'un pont biais, ceci revient à l'hypothèse des « sec¬ n° 4, p. 146-151 (1965).
tions biaises rigides ». Or, cette hypothèse n'est exacte
qu'au droit des appuis. Ailleurs elle n'est valable qu'ap¬
4. Brilla, J. : Statiche riesenie zaluziovych dôsdk (Calcul sta¬
tique des dalles articulées). Stavebnicky caso pis Slovens ké
proximativement, à moins que les poutres-caissons ne Akademie Vied Bratislava, 1963, n° 1-2.
soient dotées de diaphragmes rigides parallèles aux lignes 5. BÛCI, B. : Staticky vypocet montovanych mostov s klbovym
d'appui. spojenim prefabrikàtov v prieenom smere (Calcul statique
des ponts assemblés par une liaison transversale articulée
Nous prévoyons d'étudier encore expérimentalement sur entre les éléments préfabriqués). Inzenyrské stavby, 1963,
modèles réduits toutes nos hypothèses de travail ; nous n° 2, Bratislava.

185
6. Cejtlin, S. Ju. : 0 sovmestno) rabote panetnych perekrytij sans entretoises). Sbornik Vysokého Uéeni Technického
(traité sur l'interaction des éléments de toitures formés de Brno, 1968, n» 2.
panneaux). Sbornik statej VII, Gosstrojizdat, Moscou 17. KrajöI, L., Smrôina, M. : Spolupôsobenie nosnikov monto-
(1957). vanfch mostov (Interaction entre les poutres préfabriquées
7. Cusens, A. R., Pâma, R. P. : Design of Concrete Multibeam des ponts), Inzenyrské stavby, 1962, n° 5, Bratislava.
Bridge Decks. Proc. ASCE, 91, n° ST5, p. 255-278 (1965). 18. KWstek, V. : Vypoâet spolupäsobenl komûrkovfch nosnikû
8. Duberg, J. E., Khachaturian, N., Fradinger, R. E. : montovany'ch mostù (Calcul de l'interaction entre les poutres-
Method for Analysis of Multibeam Bridges. Proc. ASCE, 86, caissons des ponts). Inzenyrské stavby, 1962, n° 7, Bra¬
n" ST7, p. 109-138 (1960). tislava.
9. JiROUäEK, J. : Statickê vySetfovdni sikmych mostnîch ialu- 19. Nasser, K. W. : Design Procedure for Lateral Load Distri¬
ziovfch desek (Analyse statique des ponts obliques à dalle bution in Multibeam Bridges. Prestressed Cone. Inst. Jour.,
articulée). Sbornik Vysokého UCenI Technického Brno, 10, n° 4, p. 53-68 (1965).
1966, n° 2-3, p. 219-253. 20. Pama, R. P., Cusens, A. R. : Edge Beam Stiffening of
10. JmousEK, J. : ReSeni Sikmfch mostnîch îaluziovych desek Multibeam Bridges. Proc. ASCE, 93, n° ST2, p. 141-161
(Calcul des ponts obliques à dalle articulée). Inzenyrské (1967).
stavby, 1967, n° 2, Bratislava. 21. Pool, R. B. : An Investigation of Joint Forces in Multibeam
11. Jirousek, J. : Some contributions to the Finite Element Bridges. Unpubl. Ph.D. thesis, Univ. of Illinois (1963).
Methods of Stress Analysis. 10th Science Week University 22. Pool, R. B., et al. : Analysis of Multibeam Bridges with
of Damascus, Part 4, p. 21-76, Supreme Council of Sciences Beam Elements of Slab and Box Section. Eng. Exp. Sta.
Syrian Arab Republic (1970). Bull. 483, Univ. of Illinois (1965).
12. Khachaturian, N., Robinson, A. R., Pool, R. B. : 23. Ulicku, B. E. : Prostranstvennyi rascet bezdiafragmennych
Multibeam Bridges with Elements of Channel Section. Proc.
strojenij mostov (Calcul des ponts sans diaphragmes). Avto-
ASCE, 93, n° ST6, p. 161-187 (1967). transizdat, Moscou (1963).
13. Klimes, M. : Montované silnicni mosty (Ponts routiers 24. VItek, J. : Spolupùsobeni nosnikû montovanych mostù
formés d'éléments préfabriqués), Inzenyrské stavby, 1967,
(Interaction entre les poutres préfabriquées des ponts).
n° 12, Bratislava. Inzenyrské stavby 1961, n° 12, Bratislava.
14. Kopecky, J. : Statickê vySetfovâni montovan$chpravoûhlych 25. Watanabe, E. : Study of Load Distribution in Multibeam
desek-polorostù s kloubovym spojenim (Analyse statique des Highway Bridges, Unpubl. M.S. thesis, Iowa State Univ.
dalles articulées rectangulaires formées d'éléments préfa¬
(1968).
briqués), Stavebnicky îasopis Slovenské Akademie Vied
Bratislava, 1965, n° 6.
15. Kopecky, J. : Vfpoâet provoûhlych mostnîch poloroStù
metodou rozdëlovacich âisel (Calcul des ponts rectangulaires
à poutres multiples sans entretoises par la méthode des Adresse de l'auteur :
coefficients de répartition), Sbornik Vysokého UCeni Tech¬ Dr J. Jirouäek,
nického Brno, 1965, n° 3. Chaire de statistique et de résistance des matériaux, EPFL
16. Kopecky, J. : Statiche feSeni Sikmych polorostù mostniho Avenue Dapples 23
typu (Calcul statique des ponts obliques à poutres multiples 1006 Lausanne

Bibliographie Ecole polytechnique fédérale

The science of rock mechanics. Part 1 : The strength
de Lausanne
properties of rocks, par W. Dreyer, Dr-Ing. Professeur
à l'Université technique de Clausthal. « Series on Rock and Conférences
Soil Mechanics », Volume 1 (1972), No 2, janvier 1972.
Trans. Tech. Publications. P.O. Box 9787, Bay Village, L'Institut de la construction métallique nous prie d'an¬
Ohio 44140 (USA); D-3392 Clausthal — Zellerfeld, Adolf- noncer les conférences suivantes :
Ey-Str. 5 (Germany). — Un volume 15x21 cm, VIII- Le mercredi 17 mai 1972, à 14 h. 15, à la salle de confé-
501 pages, 137 figures, 86 tableaux, 200 références. Prix :
relié, 25 dollars US. (International Standard Book Number rence de l'Institut de la technique des transports, 9, ch.
(ISBËf: O-87849-002-7.) (Library of Congress Catalog des Délices', Lausanne : Behavior of Steel Frames Under
Card Number: 78-149276.)
Cyclically Applied Lateral Loads, par le professeur
Ce volume est le premier d'une série de cinq consacrée Le-Wu Lu, de l'Université Lehigh à Bethlehem (Penn-
à la « Science de la mécanique des roches ». Il traite plus -"$mm&me).
particulièrement des propriétés de résistance des roches.
Après avoir exposé quelques généralités sur les carac¬ Le lundi 5 juin 1972, à 10 h. 15, en salle B 304 : Accidents
téristiques mécaniques des minéraux (déformations, cons¬ constatés sur les ponts existants, par M. Ed. Rey,
tantes élastiques), l'auteur aborde la mécanique des roches adjoint du Service fédéral des routes et des digues.
proprement dite, aux points de vue à la fois théorique et Le mercredi 7 juin 1972, à 10 h. 15, à la salle de conférence
pratique. l'ITEP, 9, ch. des Délices : Problèmes relatifs à la
D montre en particulier quels sont les procédés de déter¬
construction en bois et en particulier en bois collé, par
mination des caractéristiques des roches sur échantillons
et donne des valeurs numériques de ces caractéristiques M. G. Kaempf, Holzbau AG, à Rupperswil.
pour différents types de roches. Il décrit également les
méthodes d'essai in situ, en tenant compte des progrès
les plus récents réalisés dans ce domaine. Quelques cha¬
pitres de son livre sont aussi consacrés aux essais sur Congrès
modèle en mécanique des roches.
Signalons, pour les lecteurs que la question intéresse, XXV« Congrès international CEBEDEAU
les titres des quatre autres volumes formant ce traité :
Liège et Gand, 16-19 mai 1972
Bins and Bunkers for Handling Bulk Materials. — Engi¬
neering Geology and its Application in Civil Engineering Le programme comprend des conférences, discussions
Projects. — State-of-Stress Measurements in Rock Masses. et visites d'étude, sous les thèmes Eau - Air - Corrosion,
— Slope Stability in Soil and Loose Rock. entre autres :

186