Pastorel Henri MScA 1977
Pastorel Henri MScA 1977
Pastorel Henri MScA 1977
Henri PASTOREL
.Sherbrooke, Québec
remorque.
de ce travail.
modeles miniatures.
REMERCIEMENTS
rement Ie Pr. Paul BOURASSA qui a ete pendant deux ans man
SOMMAIRE................................................ .ii I
REMERCIEMENTS............................................iv T
la
INTRODUCTION..............................................1 "
1. LE PNEU............................................. .3
6.4.1 Generalites........................53
CONCLUSION............................................90
VI 1
APPENDICES
BIBLIOGRAPHIE....................................... ..117
LISTE DES ILLUSTRATIONS, GRAPHIQUES ET TABLEAUX
Partle A
R .resistance au roulement
sur Ie sol.
v vite.sse du vehicule
a : angle de glissement
Partie B
remorque
remorque
m^ : masse placee en A
m9 S ^P PP
rrT : a/ T = pr T" = ^' V~ =
remorque)
xin
a 1'axe de la remorque
(angle de lacets)
remorque
a 1'axe du tracteur
m. masse du tracteur
m, masse de la remorque
AR1, BR1
charges respectives des trois essieux
et AR2
AKS1,
raideur statique laterale des pneus (ici
AKS2
AKS1 = AKS2
u vitesse de traction
Ie calcul.
BK : 0.13 - AK
point d'attelage
gravite de la remorque
du tracteur
AMU1 et
correspondent a p dans la partie theorique
AMU2
AKS1 et
correspondent a Ks darts la partie theorique
AKS2
AM 4 et
correspondent a m-, et m^ dans la theorie
AM5
H(J) et
correspondent a h et ^ dans la theorie
AL
AK et BK correspondent a a et b dans la theorie
All et
correspondent a I., et I., dans la theorie
AI2
AR1, BR1
correspondent a R-, , R^, et R., dans la theorie.
et AR2
INTRODUCTION
ete conduite avec des donnees (13) prises sur des vehicules
petits pneus.
LE PNEU
1.1 Generalites
angle de 40° n'est plus que 20° et les cordes ferment une
Ie commerce.
prieu, ceci est evidemment une supposition qui n'est pas en-
500
TrfjTn -->
^00
^
•50C3
Force
^/
-»-• -
t& lo 30 l\o So 60 70
$00 ^
^
a oo <•••
^
FoR.CE
idle .-
^ con '/c n»- ionnc \
dct^ ^
t1a d i«
(00
/
^
/
'/'
+ •+
5 6 5 I?- IS ? l\
25 T
t<^-
con^n^i ohncl
—- r4<j"id|
^ >s
c. o<-> p ><•
da^-
10
0-> 0-^f^D
forme:
ft C(
w h
^>
^
^
. Figure 2. Resistance au
,0- - /
roulement sur
1. un pneu ,
-^
^
d u C o h) Y'<^ &f
3^
IOOB y r c.si"$h4n<-e
io
<A\J rodVc.T'c.nt" C ol'o(^
Cod V> Ck<;nV-
»5
le.TL^i\'^^
d,u rc>u'c-^e-n,
\-0-^}^$/2'° pr&55top Je c^onFiqqe
IQ
^' (10,5^ J
n
^0
^^
^ - ^'
^e-'^
Vi'l-e^s c.
conventionnels.
pneu.
Notation:
M.A.T.
contruction.
chacune d'elles.
^
So-
F-^
\0
30
10
^ Fy F« = FpCoAX+F^k'r
ID
F<<
-<—I—h ^—I—I—I—»—>—>.
o ia.^i(6«-?83io<><
deformation
laterale zone de glissement
e lateral
ligne de deformation
c1 es bandes d ecoulement
(caoutchouc)
ll,gne de deformation
de 1 enveloppe
-> plan central de la
jante.
zone de contact
vehicule.
en virage. C'est sur elle que nous aliens nous pencher ici,
virage.
04 ' —T
D'apres les trois cas vus ici, on peut done donner comme
sement a. '
18
W = charge verticale
]j = coefficient de proportionnalite
a = angle de glissement
^ r ^ n
F^,.-^^ ='A/tvJ cl'°' ^ =-^L3'
Ls
Cette relation etablie pour des problemes de tenue de route
19
coefficient.
ture est infiniment moindre que celui d'un essai sur modele
grandeur nature.
vrais pneus.
pas reglable.
pn,eu,
disque en
rotation
tambour
Figure 11. Montage "type I" Figure 12.. Montage "Type II"
disque en rotation tambour en rotation
pneu
bande de roulement
/7/' / / / '/•/">"/ f / / / ^
85.
=. 7^2, ,0 ^ m
va s'appuyer sur les poutres Al. Sur ces poutres sont col-
lees des jauges qui vont servir a mesurer les forces late-
.. ^ S>-'.i.">." : • •/
;;;,..-i^<^< • - ..••
>s' - -,
^••1' '•..'••
's^-'^" ::;
."iir^w'..' "•:.' ..:• • .
' >^^:<!u't; '• •»• s>! . •' " ;' '
'.-•'•iv. i- -.•...•.•••
22 23
• ••'^^•'•.-.^.-'
SZ1
7X
•' _^..w
w 3-
^f
-^ ,—*..
k^,
^
-1-''. -I
^
7ZL
.(•'•' •<
Z.2.I
^'>-
^
.'t--'
'•• ••/:
/
7 /\.r ~/ /
'^^':
Iz-.r^—J-- ^r/ / //
;^.^..>. ' • i
l^^!^'':'i. I: :
••'/..'
'A';'''',-:l •-. ' S " Fieur? 15. Dessin du montage d'essai vu de face I
r •
;:..'•-.'.'s.; ' ;•• _•: ^^\ _'-. .....-'•I" .
28
la valeur de p.
W= 17^
^3^5
v/. lLi-65tV/i
^
p.
OQ
-^ ^= l'Z.65^ c
n
(D
I-1
^1
^- . 10.65^
-^4?^ 3 ro hrl
H- rt 0
b 1-{
H- (^ 0
&> f0 (D
rt
c \-t I-1
K M
,5.55^ G> p3. rt
p f0\
OQ ^
{U
Jv/= 8. 65 ^ I-1
(D 1-'
ro
0-
ro. (U
A-^5-1 p
OP
I-1 h-h
p- 0
en p
m 0
ro n-
9 H-
^ .^ (D 0
y 3
rt-
d.
Efl (D
c
K I-1
&>
2,.?T5 a
p 0
y
^ co
y 1~{
n? (N
£ (D
\-^s
0.32.5
•T'!1T
30
a 10
a 2 4 6 8 10 12
y
Fig. 19a 0. 15 0. 14 0.13 0. 11 0.094
p
Fig. 19b 0. 065 0. 065 0. 062 0.054
31
celle mesuree sur des pneus reels. On peut done dire qu'il
un vehicule.
suivantes:
32
5, 55
82
SLIPANGLi ^ • 1000
y:'
^
,'\
;;1' 's
\.\:< Y>*,1
•'•TIRE ^'' ;;;T.';'. 'J''.l
.;<
Z8 f'SI tNFLATtON'PREt y:
.li}
&-It.
''.?,. , ?,60'1'(JIRE;,'.,|'^ /',.>'''y
,';*v.li
11.:
' . ' .''' .^' " 'I •" 'I ^'1; ^ ^'..',,'1'1. ^^ 11/\'.'^ l.(
i fig.U-ILateral force voriiuaanglaanct normal forca /, \ i •,1'1' ';''.iTt<:. f',\'' •^ ^'•:^{
kl./';T.'!11' :••
^'''•;';,;r'; I ' ^_i\ ^ .'. . • ,',' .•'.:,'. I '^'^1.'.' ^: 'lt:'.''tl}'^./;{
6°<o<<<8c> VA/^-I.^^^
tant. Dans les essais effectues ici sur ce petit pneu, les \
sur une roue lors d'un vlrage ou d'un eca.rt du soit a. une
(Fig. 21), nous ont donne des resultats tout a fait sembla-
'f^ofc^ I at-e'rq^
^
Figure 20. Comportement de la force laterale
dans Ie cas de fortes charges
imposees au pneu
5.U
c< ^ S
^•^
{/H h- ^
I 0 rv zo ^' ^n ^ w
1<;Ki
Figure 21.- Force laterale en foriction de lrangle de glissement et de la charge
sur Ie pneu - essai sur un pneu miniature
t^
? or ec.
lo ^ rule- ^. '|r.C°( ^
2B.^i ^
3.55
\z FoR-cE nop.mdLe
11
Q f\ ^le- dc-
t'.S.S. C m C. n
5.5! ^
L
5.75 ^
L ^5
2.2,5 N
^
a\
r°T:rT"TT7"' '"T
B - ETUDE D'UN PROBLEME DE STABILITE
CHAP ITRE 4
ment tenue par un ressort, DIBOLL et HAGEN (3) ont trouve une
5.1 Generalites
appelee m^ et la seconde m.
40
remorque
remorque
m^ = masse placee en A ,
:2 = a, ^ = p, -^ = K*, ^ = I
•m'-i _ . 4
coordonnees
tion de liaison:
^ - Xl
Gr ^
^
D'autre part, cos 6 — 1 ce qui a pour consequence qu'on
central de 1'essieu.
2.4 de la partie A.
ge able.
vante:
v
et la tension du ressort: "V\ =. - k Xl
^ A 0 - ^Yx
on a:
r- L[ m ^
^
t^ \>t
x' +- ^
t T
'<! -^
•I- ^
T )^.
l)- ^ (^ IJ]
2- •.
^. •/'
-»- x
L
^ ^ r J_ ^1 ou u est r energie potentielle.
z
^T
^ i-( ^T ^ ^u a
Itr*^ ^ •^x ^\
ou Q est ici la force laterale
Sl s. = -Z f<s on obtient:
<W1 2,
, ^ £A= k h.r
r''
^ ^ - p '
^
^ -£ ^
•>(
et H- ^
c
on peut alors ecrire Ie systeme
)-k
^p;^(l-f) °
^fl;^'-P) 0 [/1+ Tl +1 H+ v 0
ks
H --0
-\
K'
5.5 Resultats
un amortissement. i:
mention plus haut ont ete effectues avec une remorque de-
tion 1^ < ^.
CHAP ITRE 6
6.1 Introduction
systeme.
nees.
remorque).
J=J5
-^
L ^ 3>J> J
^- -^- -X." ~~ -9\
I
t ov.''*^~r cA -C^HC-ICXE^C.
I
I
-L-
G-r
^
^ -.^
f\
tr[^ y
><l ^^
-^.) . , f^)
rf^
U3
r"'''!'i"i''-''----''-'--l-'rar;a"::"iT'T''
50
a 1'axe de la remorque
(angle de lacets)
remorque
a I*axe du tracteur
cette partie).
m^ : masse du tracteur
equations de liaisons:
tions:
6.4.1 Generalites
deux cas. Un premier cas oil les charges supportees par les
du tracteur
notation:
T
^~\ u des pneus avants du
-s.
tracteur
I e^ = deflection laterale
^ (•<.),
I
des pneus arrieres
- J. - ->•<- ^- — ^ du tracteur
1
'f'^)
la charge est identique
sur chaque essieu et
Figure 25. Vecteurs vitesses
egale a R-,
auxquels est soumis
Ie vehicule tracteur
vecteur resultant.
u'
Figure 26a. Vecteurs vitesses sur Figure 26b. Vecteurs vitesses sur
les pneus avants du les pneus arrieres du
tracteur tracteur
0</ =
\/ + c,r
d
u u
0
V - l.r C.-L
et que ^ r
(J u
55
notation:
e^ = deflection laterale
du pneu de la remorque.
V/-^r' a,
^ -2» ^ "— - c-3
u/ u/
V . ^Y ^^-^rj - 6^ -t)
0
Csachant d'autre part que r' = r - 4Q
rale statique d'un pneu. (Dans notre modele les six pneus
'</-L>r i^
. ^cl(V^(J^+I,-t , ^«+^ ^R,(^r- ^.')-(t>-^)^l(/
-----C5)
^^c^^-^^x^-^j^^c^^0-^
& 1 -v. ei,
_ yvM^^UrJ 4- ^>^TL. - ^ - ^^)
c) Etude de stabilite
.c ^3
on obtient Ie systeme suivant:
.... (:Zo)
. _
L
^62. ^.Z 4-^6^^.^ 4- Cc^^ + C^^^?.+C63»^-a> 4-C6fc=<-€
o .-_-.-. c i^)
dans 1'Appendice 5.
[<]S^[^J
possible.
R; P.I hR.
^v.
^1
/^ /-/ / ,7
soulager cet essieu (si Ie. centre de gravite G^, est en arriere
ce chapitre.
essieux du tracteur
des roues avants. Nous avons ici traite ces deux cas, car
gravite du tracteur.
61
{,(.,). ^,t^_'-£-]
^,) - ^^<^_ i^j
{^)-^C^ ^-^-c^^^).^.]
b) Reecriture du systeme
c) Etude de stabilite
lors du calcul
arrieres du tracteur
RESULTATS 1.
/isr^C-^
0.1 +
1
KCO)
K IB11
r^C>}
10,l|
s
^
'.7
\6
i-a h
2, I
RcCA)
/fa
^
t.'
(o,(f
AN' = loo rp m
^(^)
^\
*.6
?. 1
R<,)
,6
ft ^ = 5oorpn">
l°,ll
-!»U]
|0/tl
,0.1
3 1^ '
^(^)
A.
1°1
/^r^z iooor^"i
20 rpm (33a), 100 rpm C33b) , 500 rpm (33c), 1000 rpm (33d) .
talement) .
sinusoidal.
71
RESULTATS 2.
reels negatifs.
w/y
5 •,
0.1
^^}
->
0.1
f\^ -s ^.oa^pm
0.1
&.1 .e.J.
^e.CA)
"T:T(T"T
74
.0. Pl
o.ooT
i>rti
^eC^|
o.75
o.7
o.6
ftS
->
0.0»<> o-*»»8
0.&Z7 o.o3<. G.o<(5 rr> PO
75
Ml = 0.6 Kg).
RESULTATS 3.
stabilite du systeme.
76
la realite.
^^<
^0,11
r^i>)
.7
0-t
.. -2
O.|
P.e<^
A ^ =• tiQ f \3 r^
t0,(<
•^
A8
7
\
\<,
-,s
0. t
-°-^
I I 5 3 A
^
^
A^= lobf ^>»yi
A
<y
9/
(0,If
Figure 37b. Comportement de la valeur propre critique dans Ie plan complexe pour chaque *^
position de la masse de chargement de la remorque a differentes vitesses, co
mais dans Ie cas de charges inegales sur les roues avants et arrieres du
tracteur
io,(!
o./ 4
|l2.3 0.1
^•C.n
i?
1& ^=5oo rp^,
1^
l°,l(
Figure 37c. Comportement de la valeur propre critique dans Ie plan complexe pour
chaque position de la masse de chargement de la remorque a differentes -^I
\0
vitesses, mais dans Ie cas de charges inegales sur les roues avants et
arrieres du tracteur
•wm-.r 'r":-T
^-CM.CI).
iVi
s
^
il3
CA)
co
0
"T!-i"(']T"fl"
81
dice 7).
conclusions theoriques.
ensemble "tracteur-remorque"
meme que celle utilisee pour les essais sur les pneus
co
en
T""B"
86
6.7.3 Conclusions
ces vitesses mais nous avons note que des qu''un mouvement
de depart.
R
CONCLUSION
remorque.
reels.
APPENDICE 1
"•«-'-t'.T"~" ^••-r/-yw.';./--
^ -•T"T •
k ...~ ..
f • £ •
V ."• —-
tt •
\<
;\ '^
c
7. w
D-
<
r-i
31
z
r^-
0
VI
0
>
-a
tn -a •
X r-
0
30 c
0 n
m
6 VI
Tf
chMtai
l.^;?!?'?
11 -3 ^ ; -^ I 7
•
'.^ •t •',' •
; ; i., v-: •!<! ^ -' ,, •^
j?!'rl:^i^;; '''
,r2 ^H^?ii;!- •.<"•
^ i 2-i ^ • • •
^!^;r!
-1
^ ')iri5|^ , ^ ^:'^;J^''^^;^>\^';;^'.f^^:^^ ;
M" '•< ^—-^, •.', ,"•• • n——n ;•''.', .•'"-'.' I '.;-''."•' ''" ' ''.'.'•t. '•';:"'.' •'. ^'''':' ^•'•''•••^'.; •••;>,"\ •'...:-'•••..:./•;^''^^.;;-.'^ : :'
' /<-
•^
4-J ^;:'.'; .''.^^''^ •< • i:-•••'-• "••.•/l.^,^":'<' .; •;
•<• >
^
^V77
i
—--<->' ,.t.._r,:r. ;..j ._'f-uL... ].j»s l.^vi.i.
^-L///.
T • •
.'•V-^ r— ^»r-.«,©
<
.©
't^
\ '"-.?.>,/,
I. s
I. ttS
S 13
1. /E5
t'.l
.T":.l"T"^ ^•—.<^.-^0»
t 100-M ^ - ,0"
7.^rc. W ^ C^
^
i-
'o'jt-r.- 0-t75 :.,
<it3/IC -^:lt-<
-E^
C/.Q,1«^.. .US
I. 1.1 I . 0..g I,
i-
?-
M fTl ". i f(L> £Hf(o»'( :»Citt
-^
»-• ^ . S/l-
y
\-
< ..
"rrTrr'r FT
APPENDICE 2
largeur b = 9 mm
epaisseur h = 3 nun
longueur 9. = 105 mm
mation:
si £f = deformation
Or = contrainte
Mp = moment flechissant
on a:
^-^ — ^^t
M
96
or n rt{
1- -.
= P.
r- -T
1 &* = longueur flechie = 100 mm
y =1.5 mm
^ F. -"- £.\
rt '
Ie caldul de ^— nous donne ici E^— = 9.25103 N
A*y ""-" --""- --- j^^y
obtient alors:
^.Z5 lo3 k
"^ou a est 1'angle de glissement.
(0.04 N).
APPENDICE 3
x
Si x represente Ie vecteur X. Ie systeme
Ie
d'equations (5) peut s'ecrire sour la forme:
^
t -1 €> ~\
fi^f ^(i-F)
^ °
fn (^<' o^o-o o j<| . T I -I -\ <^ 0
0i .1 1 I J I ^
K*
\t
si x = X e"*" pour que ce systeme ait une solution il faut
que:
suivante: S ^ \^ ^ \2.
QoA 4.0|^4.Gi^ A4- c?^ ^ + 1^ \ +^5 ^
0 u- ^ c^ o = ^i.fy^
k-s
q I = •<(! -<T ^z
«•^ --
•^
^s
['- ^) + t7< [•-^ T
Q 2» ^ (k:(1^ ^) + v ^(1- ^
a kr
<a 5 =
k' k
98
Q t C(
^1 qo
c\ ( A^ ^
/ 3. r
^3 02.
C»(
^^ <=1Z
Q5 ^ q a
^t q<> o 0 to
^1 «to 0 0
°^ t<L ^1 <1° °
<*3 ^t. 0, I <<G
As - q5 ^ c<2. C<L °
a5 /
^5 a<( ^2.
0 0 C(^ ^(f &
0 0 & 5 t(
o o o o Q5 J
22
A-|> 0 est toujours vrai puisque a^ = a(l-p)"m".
m2
ex etant un scalaire positif (a = :—) , a-, sera toujours
positif'7 p.
0-f)<l- ^)>0
^
.. -. -.-.0)
K*
est positif, done (1 - ^-) > 0 pour que A^ > 0, il faut done
que 1 - p > 0
^
(/i-f)>o <=^> ^j-< i => ^i.<; i (^.^^ f=Ll
p.
^
on a done t^ < i
99
a 1 qo o
^ A 2. - Q 5 ^3 ^<- C( o > °
QS a.^ ^z.
relation suivante;
I + ^f'-l 4. f >0
:eui3q-sAs n^eAnou
•o =£x= ^x = Tx
-jeqo uOinb nA) q-Tq-ed ^ q-3 q-usq-suoo q-ss n ' seo e.zq.ou su^.a
hi
A^^>^ - As-^A+^7^ r^Jf^v'l13?f^'A) -/h^^n) + C/2-3- -
^03 C-i ~J)CJ^ -A} ~" A t^3' (JP' ^) '" ^ tA^ (,3-- T-)n^ ^?o:3n] t <M
TOT
102
1-iR^ }JR,
Si AX1
u
al2 = m^
a22 = -I2
bl3 im + m.
b23 m-, . jl b33 = ^d
c^ =-a.AXl c^ = 0
c^ = b.AXl . c = 0
•13
(AX2 + 2.AX1) C44 = Ks
'23
(h.AX2 - 2JI.AX1
AXl(b-a-2d) c^ = -Ks
33
'43
-AXl
-AX1
'53
'63
AX 2
b 11 '12 '13 b
14
b
15 '16 L12
A= 0 0 0 b44 •o 0 0
0 0 0 0 b55 0 0
0 b62 0 0 0 b66 0
0 1 0 0 0 0 0
et
B =- C31 0 C33 0 0 0 0
C61 C&Z
^ 0
0
0
&
CC6
0
0
0
APPENDICE 6
R.
06»^ L - ^ xt^
05 3. = <-^?<l^
inchanges.
APPENDICE 7
a) Le vehicule tracteur
b) La remorque
a 17 cm.
^ GP.
i— ^
^e. (7-c
^I
2> 5
1 \J ^
-^
N/degre = ijR
AKS1 = AKS 2 ) ;
u vitesse de traction;
Ie calcul;
BK . : 0.13 - AK.
point d'attelage;
gravite de la remorque;
riere du tracteur;
^M^.<< +(^HL(-AHI^^O^€.S- ^)
/} H ^
ft K = o./3_<sK
b K ^ o.l^ - ^ H
b) Calcul de H(J)
j!_ALi ////
ft H5
de 1'attelage
on a: K^ =. - AH5
. HC7\
'0
. ^
Ill
fi
l-'C3-) ?-^ c^H-c^e. ^-
^
IG.R (^^
G-T
[h^-W)^ ( ^
•7—r^
1\ HI.I
^Q\^S -AH^ ^
HK 0 . I ~1- a r,
3>-&
^TQ.C-t-(
H C^
A fi^Z ^ 4^5-^f/
0- '?•
Le calcul nous.donne:
arriere du tracteur.
e) Si AR1 = BRl
On a alors:
a) Generalites
donnee.
b) Organigramme general
Donnees generates
du montage
FIN
•HATCH wATi-rv GTunrNT ^^4
\.,'. '•} -....;; ^ !;';^'.vi^1^^;1;',; If' .;':', !,< I, >' 'w1 ;.f!.. ,'^.^!';":;' <l' ly !/,i;; ;.};•?;;.;( I ^ )^-i< "(),.;• ri<'.j( ',)r,,l )j'. •1^;;^! ,•>( ^ »+.; •J(l i<
6 AKSi-=-?or'n«
7 i\K ^°=AK S 1
R ^OQO PC:AD( '-•, , 1 '-!3 ) AN, AMI , AM?
y ' in? rnpM^T("' ("Q* ?.)}
1 0 U=0 . '"I 7 •' AN,-.? . ,7. 1 4/^0 •
11 . AM^ = A M'l + 0 • 7 1.
.?..L...»,-:^--1..^<.RK.=<AM}:-. ^K +c • ° 65':il..o.«.7,.l.,)./.(.c..*.7.lAA.M.1,.,)..^...,.;-
Z ? ^K = ( '•> . 1 ^ - n K )
2? D= ( ." .0^-^+RK )
?A AT ?-=( p .?(5^''.'C . ^5i.O .35.,? ,/6. ) +0. 51.!H( J ) .. H( J) +A ^^•.' ( ( ( ^.f! /AM-^ ).i-:H( J )
.;;AM1 :{ ( ( Q .^(L)5 ;;0 .71 -nK'^.O * 77 ) /0 .71 ) ;,• ( ( P . (.65 .' f> .71- ^K >;<-» .71 ) /r; .71) )
?.~! AL="l. 1 7~K J )
2M \\ =<\MU1 ..AP 1 /U
29 X?= AMD? .'.AP P/I.J
3 0 .'' AMI = \ ^ ^
3lAMP=&M=i
3 ?. A1 ?=4-&M? :;AL
~^~>i A??^-A,T0
34 - f'-n 1 =- A M ') ;• ( P + A L )
35 m2=-X7.<^*17
^R R15=^X1
39 R 1 6 = -^ ?
40 n ? 1 ^ A j •?
^ ^ n;?? ^+11 ( j )', x? •:r, •' 7
^ ? I'L'3 :^- '\V|^ ; Al..
''-* 3 11?.^ =X I • A |.
^5 f^r^-^o :-.H( j )
^6 ' mi =M i
^7 m3-= <vii ^n
/*<4 in/4 -( ft y, < p ) ' xl
r> /» <' 1 1 " ( A ^ 1 < A M o ) , | n.x '^ I,, f) < ^ -| , ( UK - A !•. ) < V ;^ ^; D . ^ ?
I . I
r^'1 '-11( ^) ;1. .v;> I A| 1 (1<|< ..^l. ) >< 1 | A nt A I IH ( .| t( ,| ) ,l . 1 7 ;• y'~> )
57 C4 1 =-AK ;(Xl • ' 115
5P C51 -=HK ;'X1 . •
5Q C.6\=-X?;: (D+0.17)
60 C5 ??=-U'A'X'3 .
66 c^^-^-X^
67 C.r\ 3=- Xl.
e.p, c^=+x?
6^ CA4=A|<S1 ' . I.
70 C5CS.=<\KSl E
71 C^6=-AKCi? . p
72 DC") 1 ^ 1=1, 7
73 • DD 11 K=1,7
7 4 A ( T ,K ) = ) . ' E~
75 . t-i( I ,K)=") . . • ^ . - S:
76 • 11 CONTINUE • L
77 10 CnNT!NUtr
78 ^ ( ?, < 1) =n 1 ?
79 A(l,2)=R^?.
BO A ( \, ~n = A 1 ?
81 A ( 5 , ^ ) = n 1 /t
8 2 <\ ( ? , ^ ) = n I 5
83 A ( 3 , 6 ) = 'U 6
84 AC?,1) ^= H ?1 • • .
R ^ A( ,?,-?) =P??
R6 A(?,-l)-A?7
87 A( °, 4 ) =(3° A
9 A .A(?»5)=H25
R 9 <\( 2 , r-- ) =1^6
q Q A ( 3 , ? ) ^1, .
91 A ( 4 , ^ ) - R ^
^?. A( 5 , ^) =PR^
n~? <\ ( ^ t ? ) =t.</:- "1
Q 4 A( (^ .^)^U6^
9^. n( l , l }-=-c\ \
96 R(l,?)=-a2
97 l-i ( ? , 1. )=-€??
98 M( 2 , ,?) =-€??
C5C» n( 3,^ ) =+ i, •
I 00 (-)( /i ,1 ) =-C^l
101 P( <4 . 4 ) •=-0.^^
1 C ^ H ( f> , c, ) = - O'i 6
107 A(1,7)=^1^
1 0 H A ( ?. ,7 ) = Pi ? 3
\0^ A(7,7)=n^~\
110 A(7,1)=H-"1
T 1 5 /\( 7 , /i ) ^m/i
11^' A( 7 ,'-,) -11 ••"-•;
1 1 (I I • ( 7 , /)-.-(- '
1 '-I
XM ,K.) -:" .
1 n
r
-1
Y( I ,1<) -.',( T..^ )
1 ? <'» Y(I,K)^H(T ,K )
1 e:
c; "»1 c'' ^''T T ^itj""
1 ^ -^ CPMTTMU"
< ?. 7 CALL v< TMv/r'r ( x ,7,7 , r'l- T , I'- 1- , T r.; , T C )
1 ^" (- ^ L L ' ^ '- ^ r" P ( X, -;' , ^ , Y , 7 , 7,7. -^ , 7 , '.* , )
? o n C A L 1.. -" T ^ V A. I ( < , 7 , / , r ';:, V/ L'-' , r <.V Al I , I I •!') J <'. , ";'. V M \ T , \>i '~\-' \r' c , ^ " )
1 -(
l1 pr:. T v.) T -, , J , r ( J ) ,AM,U,^ '•''• 1 , ^ '/i ,""
1 ^ n --) ^ T -: 1 , -
1 -I
Tr-( T-^ln)T^( T ) .1 T. < ) '~,^ TH Z
1 f. np T \T ::i
1.
-}
,C2 ! <" H MT T 'H I"
1 -3
Q ^ rnr: '^T( • ry , i-!,^,(r-^ - , :^) )
\ ^ ',-> p prr-»/ ^T ( I, ^X , 'MnT f'ni P'l"' ' , ''-' )< , ' '•.I FT r- r\ J ND '• )
1 6. •3.
r;n Tn i, i : "i
1 Zi u r- N.n
1;C-N|- ^Y
BIBLIOGRAPHIE
May 1977. .
1968.
Aout 1971.
1968.