CRT Partie B
CRT Partie B
Lorsque le tablier d'un pont à poutres contient des entretoises intermédiaires entre les
poutres, la section transversale devient rigide. Dans ce cas, les efforts dans les entretoises
sont déterminées d'après la méthode de Courbon.
Mais, lorsque le tablier est dépourvue d'entretoises intermédiaires, qui est le cas
souvent rencontré de nos jours à raison de la préfabrication des poutres, le hourdis joue le
rôle de l'entretoisement. Ainsi, une flexion générale est née. Cette flexion représente la
flexion qu'elle qu'aurait endossée l'entretoise intermédiaire si elle existe. Cette flexion
globale (parfois on l'appelle aussi flexion générale) est déterminée d'après la méthode de
Guyon-Massonnet [1,2] en prenant 1 m.l. du hourdis comme étant une entretoise
intermédiaire. Les efforts dans le hourdis, sous l'effet de la flexion globale, sont ainsi
calculés comme dans le cas d'entretoise. ainsi, cette flexion globale est surtout importante
dans le cas des ponts à poutres sans entretoises intermédiaires.
D'après Guyon-Massonnet [1,2], le moment fléchissant dans une entretoise est donné
par:
∞
My(x,y) = b .∑ μn.q n.sin n.π.x (1)
8 n =1 Lc
Lr+2.Ltr LT
avec b: demi-largeur active = =2 ;
2
Lr: Largeur roulable; Ltr: largeur du trottoir; LT: largeur totale.
μ n : Coefficient de Guyon-Massonnet.
q n : Charge appliquée en forme de lame de couteau (développement en série de
Fourrier).
Lc : Longueur de calcul (dans le sens longitudinal).
Le moment maximum est au centre de la dalle ⇒x= Lc
2
⇒ sin n.π.x = sin n.π
Lc 2
Si n est impair sin n . π = ±1
2
Si n est pair sin n.π = 0
2
Ainsi, on ne retient donc que les harmoniques impairs ( n = 1, 3, 5, 7, …). Pour un
calcul manuel, on peut se contenter des deux premiers termes (à savoir n=1 et n=3). Ainsi,
8
(
My = b. μ1q1sin π +μ3q3sin 3π
2 2
)
b
My = . (μ 1 q1 − μ 3 q 3 ) (2)
8
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 112
sont réduites à la charge de type B [3]. Ainsi, pour l'étude de la flexion globale, on a à
considérer les mêmes charges prévues pour la flexion locale à savoir:
- Charge permanente, gper
- Charge de trottoir, qtr
- Charge de type B avec ses 3 systèmes: Bc, Bt et Br
- Charge militaire si elle existe.
Ces charges doivent être développées en série de Fourrier en une lame de couteau
(Fig. 1).
L'équivalence des charges usuelles est présentée, pour trois cas, comme suit [4]:
Lc
4q
qk= si k est impair
kπ
qk =0 si k est pair
4q
En particulier ici, k=1 ⇒ q1=
π
4q
k=3 ⇒ q3= (3)
3π
b) Charge uniformément répartie sur une longueur, 2c, centrée sur le point d'abscisse d
(t.q. les charges militaires Mc80 ou Mc 120).(Fig. 3)
2c
q
d
Lc
Figure 3: Charge uniformément répartie sur une longueur, 2c
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 113
.sin k.π.c .sin k.π.d
4q
k impair qk =
kπ Lc Lc
k pair qk = 0
d
Lc
Remarques
i
Lc Lc
3π.di
q3= ∑
2Pi
. sin (5')
i
Lc Lc
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 114
μ n = f (α, θ, y, e) est déterminé par les tables ou les formules de Guyon-Massonnet
[1,2]. Celles-ci sont analogues aux tables de K = f (α, θ, y, e) employées pour le calcul du
CRT [5] et le calcul du coefficient correspondant se fait de la même manière. En
particulier, cette fois on cherche les moments dans la fibre centrale (y=0). Ainsi, on utilise
uniquement que la première ligne (y=0) de chaque table. C'est pourquoi les tables à
employer (et les courbes obtenues par la suite) sont toujours symétriques par rapport à
e=0. Ces tables sont présentées dans l'annexe.
Le paramètre d'entretoisement est θ= n.b (6)
Lc
Pour le calcul de μ 1: θ1= b (6a)
Lc
Pour le calcul de μ 3 : θ3 = 3b (6b)
Lc
α, le paramètre de torsion, est le même que celui employé pour le calcul de K [5], à
savoir,
γ +γ
α= P E (7)
2 ρP.ρE
où γP, γE, ρP, ρE sont respectivement les rigidités de flexion et de torsion des poutres
et des entretoises par unité de longueur tels que définies dans le chapitre 3 pour le calcul
du CRT K.
De la même manière que pour le calcul de ce dernier (CRT K), l'interpolation sur α
fait intervenir μ 0 pour α = 0 et μ 1 pour α = 1 (à ne pas confondre ce μ 1 pour α = 1
avec μ 1 pour n=1). L'interpolation peut se faire d'après la relation de Guyon- Massonnet
[1,2], c.à.d.,
μ= μ0 + (μ1−μ0 ) α (8)
μ = μ 0 + (μ 1 - μ 0) α0,05 0 ≤ θ ≤ 0,1
μ = μ 0 + (μ 1 - μ 0) α(1-eθo)
0,065−θ
0,1 ≤ θ ≤ 1 avec θo= 0,663
μ = μ 0 + (μ 1 - μ 0) α θ>1
μ 0 et μ 1 sont données par les tables de Guyon-Massonnet en fonction de θ, e et y (voir
annexe). Il ne reste plus qu'à chercher la variation de μ 1= f(e) et de μ 3 = f(e) en
employant les tableaux de Guyon-Massonnet (ou les formules) présentés en annexe et
correspondant aux valeurs de α et de θ. On trace ces deux courbes comme dans le cas de K
(sur un même papier millimétrique). On charge transversalement par la charge permanente
et par les charges réglementaire (Bc, Bt, Br, Mc80, qtr) de la manière la plus défavorable.
On cherche μ 1 et μ 3 pour chaque charge.
l c = 15,36 m
bp=0,3
Déterminer les moments de la flexion globale dans le hourdis sous l'effet des charges
permanentes, des surcharges du trottoir et de celles de type B et Mc80.
Solution
1) Courbes de µ1 et de µ3 en fonction de e
Ainsi, on commence par la valeur de θ1=0,40. On copie les valeurs des tables de
Guyon-Massonnet dans les deux premières lignes et puis on interpole sur α. Les résultats
sont présentés dans le tableau N°1. Le tableau est symétrique par rapport à e =0.
θ1 = 0,40
e -b -3b/4 -b/2 - b/4 0 b/4 b/2 3b/4 b
µ (α=0).104 -2292 -1161 -20 1151 2372
µ (α=1).104 -1016 -617 -131 546 1563 Symétrique
µ 1.104 -1514 -828 -88 782 1878
Ensuite, on procède de la même manière pour la valeur de θ1 = 1,22. Mais cette fois,
on ajoute une interpolation sur θ. Les résultats sont présentés dans le tableau N°2.
θ3= 1,22
e -b -3b/4 -b/2 - b/4 0 b/4 b/2 3b/4 b
θ µ (α=0).104 -120 -190 -190 80 940
= µ (α=1).104 -77 -94 -101 8 657 Symétrique
1,2 µ 3.104 -94 -131 -135 36 767
θ µ (α=0).104 -4 -111 -174 -0,5 801
= µ (α=1).104 -41 -59 -81 -21 566 Symétrique
1,4 µ 3.104 -27 -79 -117 -13 658
µ 3.104 -87 -126 -133 31 756 Symétrique
Lch = Lr = 9,50 m;
Nv = 3 voies;
Lr ≤ 7m ⇒ Pont de la 1ère classe.
b = 6,25 m.
0,1878
μ1
0,0782
0,0756
μ3
0,0031
-0,0087
-0,0126 -0,0133 -0,0088
-0,0828
-0,1514
Figure 6: Courbes de µ 1 et de µ 3 en fonction de e.
a) Charge permanente
Transversalement
On charge toute la largeur transversale puisque cette charge existe toujours. Vu que
cette charge est uniformément répartie on détermine les coefficients µ 1 et µ 3 en prenant les
différentes surfaces positives et négatives (Fig. 7). On prend l'avantage de la symétrie en
traitant deux fois la moitié.
⎛1 -3b -b -b 1 ⎞ b
µ 1 = 2.⎜2 µ 1(e=-b) +.µ 1(e= 4 ) +.µ 1(e= 2 ) +.µ 1(e= 4 )+ 2 µ 1(e=0)⎟ . 4
⎝ ⎠
6,25
= (-0,1514 - 2.0,0828 - 2.0,0088 + 2.0,0782 + 0,1878). 4
= 0,015 m
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 118
⎡1 -3b -b -b 1 ⎤ b
µ 3 = 2.⎢2 µ 3(e=-b) + µ 3(e= 4 ) + µ 3(e= 2 ) +µ 3(e= 4 ) + 2 µ 3(e=0)⎥ . 4
⎣ ⎦
6,25
= (-0,0087 - 2.0,0126 - 2.0,0133 + 2.0,0031+ 0,0756). 4
= 0,033 m
gper
0,1878
qtr qtr
μ1
0,0782
0,0756
μ3
0,0031
-0,0087
-0,0126 -0,0133 -0,0088
-0,0828
-0,1514
Figure 7: Courbes de µ 1 et µ 3 en fonction de e et emplacement de la charge gper et qtr.
Longitudinalement
gper
Lc
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 119
4q 4.0,656
q1 = π = π = 0,835 t/m2
4q 4.0,656
q3 = 3π = 3π = 0,278 t/m2
Enfin, on obtient le moment de la flexion globale sous l'effet de la charge permanente:
My = b.(μ1q1−μ3q3 )
8
My = (6,25/8) . (0,015 .0,835 - 0,033 .0,278) = 0,003 t.m/ml.
b) Charge du trottoir
Longitudinalement
qtr
Lc
c) Charge Bc
Transversalement:
On place la charge Bc sur les courbes de manière la plus défavorable (tel qu'il a été fait
pour le calcul du CRT K et notamment pour la poutre centrale vue sa symétrie). Comme
on doit respecter la règle Nf ≤ Nv (=3), on charge une file, 2 files ou 3 files, symétrique
par rapport à l'axe transversal ou l'une des files de roues sur l'axe (cas non-symétrique). les
valeurs de µ sont lues directement sur la courbe (Fig.10).
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 120
1er cas: 1 file de Bc
1ère position: Symétrique
µ 1,i= 0,115 et µ 3,i = 0,020
1 1
µ 1 = 2 ∑ µ 1,i = 2 2.0,115 = 0,115
1 1
µ 3 = 2 ∑ µ 3,i = 2 2.0,020 = 0,020
2ème position: Non-symétrique
µ 1,1= 0,1878 µ 1,2= 0,050 µ 3,1 = 0,0756 µ 3,2= - 0,006
1 1
µ 1 = 2 ∑ µ 1,i = 2 (0,1878+0,050) = 0,119
1 1
µ 3 = 2 ∑ µ 3,i = 2 (0,0756-0,006) = 0,035
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 121
1 1
µ1 = 2 ∑ µ 1,i = 2 .2.(0,168+0,038) = 0,206
1 1
µ3 = 2 ∑ µ 3,i = 2 .2.(0,058-0,010) = 0,048
2ème position: Non-symétrique
µ 1,1= 0,023 µ 1,2= 0,150 µ 1,3= 0,1878 µ 1,4= 0,050
µ 3,1=-0,012 µ 3,2= 0,044 µ 3,2= 0,0756 µ 3,2=-0,006
1 1
µ 1 = 2 ∑ µ 1,i = 2 (0,023+0,15+0,1878+0,05) = 0,205
1 1
µ 3 = 2 ∑ µ 3,i = 2 (-0,012+0,044+0,0756-0,006) = 0,051
Longitudinalement
La position la plus défavorable est déterminée par le théorème de Barré (ou plus
simplement on emploie le tableau N°1 du chapitre 3 page 46).
Lc = 15,36 m 11,75 < Lc < 17,44 m
δ= 0,375 m 1,5 m
P P P P
4,5 m 4,5 m
2 2
d
1
d
2
d
3
d4
l c
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 122
Les distances di sont déterminées d'après la valeur de δ donnée dans le tableau.
15,36
d2 = Lc − δ = 2 − 0,375 = 7,305 m ⇒ P2 = P = 12t.
2
P
d1 = d2 − 4,5 = 7,305 − 4,5 = 2,805 m ⇒ P1 = 2 = 6t.
d3 = d2 + 1,5 = 7,305 + 1,5 = 8,805 m ⇒ P3 = P = 12t.
P
d4 = d3 + 4,5 = 8,805 + 4,5 = 13,305 m ⇒ P4 = 2 = 6t.
⎡ πd πd πd πd ⎤
q1 = 2P ⎢ 1 sin 1 + sin 2 + sin 3 + 1 sin 4 ⎥
lc ⎣2 lc lc lc 2 lc ⎦
⎡ πd πd πd πd ⎤
= 2P ⎢ 1 sin 1 + sin 2 + sin 3 + 1 sin 4 ⎥
lc ⎣2 lc lc lc 2 lc ⎦
2 . 12 ⎛1 1 ⎞
= 15,36 ⎜2 0,543 + 0,997 + 0,974 +2 0,408⎟
⎝ ⎠
= 3,823 t /m
2P 3πd1 2P2 3πd2 2P3 3πd3 2P4 3πd4
q3 = 1 sin + sin + sin + sin
lc lc lc lc lc lc lc lc
⎡ 3πd1 3πd2 3πd3 1 3πd4 ⎤
= 2P ⎢ 1 sin + sin + sin + sin
lc ⎣ 2 lc lc lc 2 l c ⎥⎦
2 . 12 ⎛1 1 ⎞
= 15,36 ⎜2 0,989 − 0,974 − 0,771 +2 0,952⎟
⎝ ⎠
= - 1,210 t /m
Enfin, les moments fléchissants sont obtenus et comparés pour en tirer la valeur
maximale. Comme les coefficients µ sont multipliés par les charges q, la comparaison
entre les différents cas doit se faire à ce niveau. C'est là aussi, qu'on tient compte du
coefficient bc pour comparer les différents cas (la comparaison se fait avec bc.My).
My = b .(μ1.q1−μ3 .q3 )
8
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 123
2ème position: Non-Symétrique
My = (6,25) . (0,18 .3,823 + 0,044 .1,210) = 0,579 t.m/ml.
bc My = 0,95.0,579 = 0,550 t.m/ml.
Donc le 2ème cas (2 files) à sa 1ère position symétrique représente le cas le plus
défavorable.
A retenir: My = 0,661 t.m/ml bc=1,1
On tiendra compte par la suite des différents coefficients, à savoir bc, γQ1 et δB.
bc = 1,1
γQ1: Coefficient de pondération du BAEL; γQ1=1,6 à l'ELU γQ1=1,2 à l'ELS.
δB : Coefficient de majoration dynamique (le même que pour la flexion locale).
d) Charge Bt
Transversalement:
Comme le pont est de la 1ère classe, on charge 1 ou 2 files de Bt, symétrique (1ère
position) ou non (2ème position). On emploie la figure 12 pour savoir les valeurs des
coefficients µ.
2,0
1cv Bt 2ème position
2,0
1cv Bt 1ère position
μ3
μ1
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 124
1er cas: 1 file de Bt
1ère position: Symétrique
µ 1,1 = 0,115 et µ 3,1 = 0,020
1 1
µ 1 = 2 ∑ µ 1,i = 2 .2.0,115 = 0,115
1 1
µ 3 = 2 ∑ µ 3,i = 2 .2.0,020 = 0,020
2ème position: Non-symétrique
µ 1,1= 0,1878 µ 1,2= 0,050 µ 3,1 = 0,0756 µ 3,2= - 0,006
1 1
µ 1 = 2 ∑ µ 1,i = 2 (0,1878+0,050) = 0,119
1 1
µ 3 = 2 ∑ µ 3,i = 2 (0,0756-0,006) = 0,035
2ème cas: 2 files de Bt
1ère position: Symétrique
µ 1,1= 0,150 µ 1,2= 0,023 µ 3,1 = 0,044 µ 3,2= - 0,012
1 1
µ 1 = 2 ∑ µ 1,i = 2.2.(0,15+0,023) = 0,173
1 1
µ 3 = 2 ∑ µ 3,i = 2.2.(0,044-0,012) = 0,032
2ème position: Non-symétrique
µ 1,1= 0,050 µ 1,2= 0,1878 µ 1,3= 0,115 µ 1,4=-0,003
µ 3,1= -0,006 µ 3,2= 0,0756 µ 3,2=0,020 µ 3,2=-0,013
1 1
µ 1 = 2 ∑ µ 1,i = 2 (0,050+0,1878+0,115-0,003) = 0,175
1 1
µ 3 = 2 ∑ µ 3,i = 2 (-0,006+0,0756+0,020-0,013) = 0,038
Longitudinalement
On place la charge Bt de la manière la plus défavorable pour un moment centrale. On
emploie les lignes d'influences et on fait avancer le tandem pour en tirer l'effet maximum.
δ= 0,675 1,35
P P
d1
d2
Lc
Figure 13: Cas le plus défavorable pour Bt dans le sens longitudinal
1,35
d1 = Lc - 2 = 7,005 m
2
1,35
d2 = Lc + 2 = 8,335 m
2
P=16t l c = 15,36 m.
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 125
2P1 πd1 2P2 πd2
q1= sin + sin
lc lc lc lc
⎛ πd πd ⎞
= 2P ⎜ sin 1 + sin 2 ⎟
lc ⎝ lc lc ⎠
2.16
= 15,36 (0,990+0,990)
= 4,125 t/m
2P 3πd1 2P2 3πd 2
q3 = 1 sin + sin
lc lc lc lc
⎛ 3πd1 3πd2 ⎞
= 2P ⎜ sin + sin
lc ⎝ lc l c ⎟⎠
2.16
= 15,36 ( - 0,915 - 0,915)
= - 3,812 t/m
Ainsi, on obtient les moments pour chaque cas. Comme le coefficient bt (ici bt=1) ne
dépend que de la classe du pont, il n'intervient pas dans la comparaison des résultats.
My = b.(μ1q1−μ3q3 )
8
1er cas: 1 file de Bt
1ère position: Symétrique
My = (6,25/8) . (0,115 .4,125 + 0,02 .3,812) = 0,430 t.m/ml.
2ème position: Non-Symétrique
My = (6,25/8) . (0,119 .4,125 + 0,04 .3,812) = 0,503 t.m/ml.
Donc le 2ème cas (2 files) à sa 2ème position non-symétrique est le cas le plus défavorable.
A retenir: My = 0,680 t.m/ml avec bt = 1.
On tiendra compte par la suite des différents coefficients, à savoir bt, γQ1 et δB.
γQ1: Coefficient de pondération du BAEL;
γQ1=1,6 à l'ELU et γQ1=1,2 à l'ELS.
δB : Coefficient de majoration dynamique.(le même que pour Bc).
e) Charge Br.
Le traitement de cette charge est plus simple puisqu'elle est représenté par une seule
roue isolée.
Transversalement
On lit directement les valeurs de µ sur la courbe (Fig.14).
µ 1 = 0,1878 et µ 3 =0,0756
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 126
0,1878
μ1
0,0756
μ3
Longitudinalement:
d
Lc
ℓc πd π
d =2 . sin =sin =1
lc 2
3πd 3π
sin =sin =−1
lc 2
P = 10 t.
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 127
πd π
q1= 2P sin = 2P sin = 2P = 2.10 = 1,302 t/m.
lc lc lc 2 l c 15,36
3πd 2P 3π 2P 2.10
q3 = 2P sin = sin =− =− = - 1,302 t/m.
lc lc lc 2 l c 15,36
f) Charge Mc80
Transversalement
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 128
1er cas: 2 chenilles symétriques
µ 1,1= 0,117 µ 1,2= 0,060 µ 3,1 = 0,022 µ 3,2= - 0,003
LMc 0,85
µ 1 = 4 ∑ µ 1,i = 4 .2.(0,117+0,060) = 0,075 m
LMc 0,85
µ 3 = 4 ∑ µ 3,i = 4 .2.(0,022-0,003) = 0,008 m
2ème cas: 1 chenille sur l'axe, l'autre à 1,95 m.
µ 1,1= 0,154 µ 1,2= 0,030 µ 1,3= - 0,012
µ 3,1 = 0,048 µ 3,2= - 0,010 µ 3,3= - 0,012
LMc 0,85
µ 1 = 4 ∑ µ 1,i = 4 (2.0,154+0,030-0,012) = 0,069 m
LMc 0,85
µ 3 = 4 ∑ µ 3,i = 4 (2.0,048-0,010-0,012) = 0,016 m
3ème cas: 1 chenille dont l'extrémité est sur l'axe, l'autre à 1,95 m.
µ 1,1= 0,1878 µ 1,2= 0,122 µ 1,3= 0,006 µ 1,4= - 0,038
µ 3,1 = 0,0756 µ 3,2= 0,026 µ 3,3= - 0,012 µ 3,3= - 0,013
LMc 0,85
µ 1 = 4 ∑ µ 1,i = 4 (0,1878+0,122+0,006-0,038) = 0,059 m
LMc 0,85
µ 3 = 4 ∑ µ 3,i = 4 (0,0756+0,026-0,012-0,013) = 0,016 m.
Longitudinalement
2c
qMc
d
Lc
Figure 17: Position la plus défavorable pour Mc80 dans le sens longitudinal.
ℓc πd π
d =2 . sin =sin =1
lc 2
3πd 3π
sin =sin =−1
lc 2
4,90
c = 2 = 2,45 m.
72
q = 4,9.0,85 = 17,29 t/m2.
4q πc πd
q1= sin sin
π lc lc
4.17,29 π.2,45
= sin .1
π 15,36
= 10,57 t/m2.
4q 3πc 3πd
q3 = sin sin
3π lc lc
4.17,29 3π.2,45
= sin .(-1)
3π 15,36
= - 7,32 t/m2.
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 129
Ainsi, on obtient les moments pour chaque cas.
My = b.(μ1q1−μ3q3 )
8
1er cas: 2 chenilles symétriques
My = (6,25/8) . (0,075 .10,57 + 0,008 .7,32) = 0,665 t.m/ml.
2ème cas: 1 chenille sur l'axe, l'autre à 1,95 m.
My = (6,25/8) . (0,069 .10,57 - 0,016 .7,32) = 0,661 t.m/ml.
3ème cas: 1 chenille dont l'extrémité est sur l'axe, l'autre à 1,95 m.
My = (6,25/8) . (0,059 .10,57 - 0,016 .7,32) = 0,579 t.m/ml.
Donc le 1er cas symétrique est le cas le plus défavorable.
A retenir: My = 0,665 t.m/ml
On tiendra compte par la suite des différents coefficients, à savoir γQ1 et δMc.
γQ1: Coefficient de pondération du BAEL;
γQ1=1,35 à l'ELU et γQ1=1,0 à l'ELS.
δMc : Coefficient de majoration dynamique pour Mc.
Tableau récapitulatif
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 130
Références:
[2] R. Bares et Ch. Massonnet, "le Calcul des Grillages de Poutres et Dalles
Orthotropes", Dunod, Paris 1966.
[3] M. Ben Ouézdou, "Cours d'Ouvrages d'Art", Polycopié ENIT, Première version,
1993. Code ENIT: 358 PE.
[4] J. A. Calgaro, "Calcul Pratique des Dalles Minces", Master Ouvrages d'Art, ENPC,
1987.
[5] M. Ben Ouézdou, "Etude de la Répartition Transversale des Charges sur les Ponts à
Poutres par la Méthode de Guyon-Massonnet", Polycopié ENIT, Première Edition,
Oct. 1992. Code ENIT: 352 PE.
______________________________________________________________________________________
M. Ben Ouézdou Chap 5, page 131