GOODYEAR AEROSPACE ~UASI-CB- 120261) USERS IIANUAL: DY YALIICS U74-29309 OF THO BODIES COUYECTSD dY AN ELASTIC rETBEI, SIX EEGBAiLS OF FAiiZDOn FOREBODY ,192 FIVE DEGREES OF (Goodyear Aerospace Unclas :ore.) 139 p Hc S10,UO -- SSCL 20K 63/32 16556
GOODYEAR AEROSPACE
~ U A S I - C B - 120261) USERS IIANUAL: DY YALIICS U74-29309 OF THO BODIES COUYECTSD dY A N ELASTIC r E T B E I , S I X EEGBAiLS OF F A i i Z D O n FOREBODY ,192 F I V E DEGREES OF (Goodyear Aerospace Unclas :ore.) 139 p Hc S10,UO
--
SSCL 20K 63/32 16556
CODE IDENT NO. 25500
GOODYEAR AEROSPACE CORPORATION
AKRON 15. OHIO
USERS MANUAL
DYNAMICS O F TWO BODIES CONNECTED
BY AN ELASTIC TETHER - SIX DEGREES OF FREEDOM FOREBODY
AND F I V E DEGREES OF FREEDOM DECELERATOR
(REF. NASA CONTRACT NAS~-29144 s / A ~ )
GEORGE R. DOYLE, JR. &
JAMES W. BURBICK
APRIL 1974
ABSTRACT
One important a s p e c t t o recover ing a body f a l l i n g through the atmosphere, is t o
d e c e l e r a t e and s t a b i l i z e i t . This is u s u a l l y accomplished by means o f a parachute .
The d e s i g n of t h e recovery system n e c e s s i t a t e s a knowledge of t h e dynamics and
loads dur ing parachute deployment and i n f l a t i o n . I n many c a s e s , a p i t c h p lane
a n a l y s i s w i l l p rovide adequate informat ion. However, i f t h e body is i n e g e n e r a l
tumbling motion, i t is necessary t o analyze i t s motion i n +Free dimensions.
T h i s r e p o r t c o n t a i n s t h e equa t ions of m t i o n and a computer program f o r t h e dy-
namics of a s i x degree of freedom body joined t o a f i v e degree of freedom body
by a q u a s i l i n e a r e l a s t i c t e t h e r . The forebody is assumed t o be a completely
g e n e r a l r i g i d body w i t h s i x degrees of freedom; the d e c e l e r a t o r is a l s o assumed t o
be r i g i d , bu t w i t h only f i v e degrees of freedom (symmetric about i t s l o n g i t u d i n a l
a x i s ) . The t e t h e r is represen ted by a s p r i n g and dashpot i n p a r a l l e l , where t h e
s p r i n g c o n s t a n t i s a f u n c t i o n of t e t h e r e longa t ion . Lagrange's equa t ion is used
t o d e r i v e t h e equa t ions of motion w i t h the Lagrange m u l t i p l i e r technique used t o
express the c o n s t r a i n t provided by the t e t h e r . A computer program is included
which provides a time h i s t o r y of the dynamics of both bodies and t h e t e n s i o n i n
t h e t e t h e r .
...
GER- 16047 .* 9 ,
TABLE OF CONTENTS
Abstract
L i s t of Figures
Naecnc la t u r e
Chapter
I Introduct ion
I I Derivation of Equations of Motion
1. Coordinate System
2. Euler Angle Transformation
3. Kinet ic Energy
4. Po t en t i a l Energy
5. Rayleigh's Diss ipa t ion Function
6. Lagrange ' s Equation
7. General Equations of Motion
8. Simplified Equations of Motion
9. Generalized Forces - Aerodynamics
I11 Computer Program
1. Features of Computer Program
2. Input
3. Output
4. Numerical Solut ion
5. P l o t t i n g Routine
6. English t o Metric Conversion
7. Conclusions and Recornendat ions
IV Program L i s t i ng and Sample Computer Run
Page
1
iii
References
>- .f
GER- 16047 k
- M
Figure No.
LIST OF FIGURES
T i t l e - Coordinate Systems
Euler Angle Rotations
Aerodynamics Coordinate System
Input Parameters f ~ r Parachute Added Air Mass
Parameters for Program Calculated Parachute Drag Area Time History
Page
The fol lowing is a l i s t of variables used i n the computer program and i n the
derivat ion of the equations a s discussed i n t h i s report. A br i e f descri , , ion
and associated units are included.
GER- 16047
METRIC ENCLISH FORTRAN STANDARD DESCRIPTION UNITS UNITS
A a Distance along t h e l o n g i t u d i n a l a x i s of -4
AALPDE (8)
AALPPE (16)
AALPME ( 16)
AALPPE (8)
AAM(8)
AAI@ (8)
AERATO
t h e forebody (5) from t h e i n t e r s e c t i o n
of t h e body axes t o t h e te ther- torebody
confluence p o i n t , p o s i t i v e towards the
nose (m) f t
Dummy v a r i a b l e s use3 t o express incre- (mlsec) f t l s e c
a e n t a l v e l o c i t i e s of t h e forebody i n
t h e Runge-Kutta i n t e g r a t i o n
An a r r a y of e i g h t va rLab les s i g n i f y i n g
angle-of-a t tack of t h e forebody used
wi th damping c o e f f i c i e n t s
An a r r a y of s i x t e e n - a r i a b l e s s i g n i f y i n g
angle-of-a t tack of t h e forebody used
wi th fo rce c o e f f i c i e n t s
An a r r a y of s i x t e e n v a r i a b l e s s i g n i f y i n g
angle-of - a t tack of the forebody used
wi th moment c o e f f i c i e n t s
An a r r a y of e i g h t v a r i a b l e s s i g n i f y i n g
angles of a t t a c k of t h e d e c e l e r a t o r
An a r r a y of e i g h t v a r i a b l e s s i g n i f y i n g
Mach number of t h e forebody used wi th
f o r c e and moment coe f f i c i e n t s
An a r r a y of e i g h t v a r i a b l e s s i g n i f y i n g
Mach number of t h e f o r c ~ o d y used wi th
damping coef f i d e n t s
An a r r a y of e i g h t v a r i a b l e s s i g n i f y i n g
Mach number of t h e d e c e l e r a t o r
Suspension Line AE Rat io
(AERATO = AE/AENy 1
GER- 16047 IY.
t . FORTRAN STANDARD DESCRIPTION - AIPHI
AIPHID
AJALPD
A JALPF
A JALPM
AKAM
ALPE 0(
Number of elements i n PPHIE a r r a y
( = 2 t o 8)
Number o f elements i n PPHIDE a r r a y
( = 2 t o 8)
Number of elements i n AALPDE a r r a y
( = 8)
Number of elements i n AALPFE a r r a y
( = 8 o r 16)
Number of elements i n AALPME a r r a y
( = 8 o r 16)
Number of elements I n AAM a r r a y
( = 2 t o 8)
Number of elements i n AAMD a r r a y
( = 2 t o 8)
Angle-of-attack of t h e forebody
ALPPE aP Angle-of - a t t a c k of the d e c e l e r a t o r
AM Mach number of t h e forebody
AMAX l Larger o r d i n a t e o f two points on the
I-- -0 M Y - AMAY 1
l o n g i t u d i n a l added mass versus
Do l o g log p l o t
METRIC ENGLISH UNITS UNITS -
kg s l u g
Smaller o r d i n a t e of two po in t s on the
1ongI tudinal added mass ve rsus
Do l o g log p l o t kg s lug
Larger o r d i n a t e of two po in t s on thc
l a t e r a l added mass versus
Do log log ~ l 0 t
-vi-
kg s lug
COODVEAR AEROSPACE C 0 I . 0 1 . 1 6 0 .
GER- 16074 - -. --
I .* r J3zmau STANDARD DESCRIPTION
AMAY 2 Smaller o r d i n a t e o f two p o i n t s on the
l a t e r a l added mass versus
METRIC ENGLISH UNITS UNITS
AMP
kP
Do log l o g p l o t
Mach number of t h e d e c e l e r a t o r
Dis tance a long t h e l o n g i t u d i n a l a x i s
of t h e d e c e l e r a t o r (X ) from t h e c.g. pb
t o t h e t e t h f r - d e c e l e r a t o r con£ luence
po in t
Exponent of l o n g i t u d i n a l added mass
equa t ion
(MPAL = RHOOO * BX * DS ** AX)
Exponent of l a t e r a l added mass equa t ion
(WAS = RHOOO * BY * DS * AY)
Dis tance a long the l a t e r a l a x i s of t h e
forebody (Yb) from t h e i n t e r s e c t i o n of
t h e body axes t o t h e te ther- forebody
confluence p o i n t , p o s i t i v e towards t h e
l e f t wiqg
Dummy v a r i a b l e s used t o express inc re -
mental v e l o c i t i e s of tile d e c e l e r a t o r i n
t h t Runge-Kutta i n t e g r a t i o n
C o e f f i c i e n t of l o n g i t u d t n a l added mass
equa t ion (MPAL " RHOOO * BX * DS ** AX)
C o e f f i c i e n t of l a t e r ? ? added mass
eque t ion (MPAS = RHOOO * BY * DS *" AY)
Dis tance a long the v e r t i c a l a x i s of t h e
forebody (Z,,) from the i n t e r s e c t i o n of
t h e body axes t o t h e te ther- forebody
conf luence p o i n t , p o s i t i v e up
kg s l u g
GER- 1604 7 4.P.
$ I METRIC ENGLISH FORTRAN STANDARD DESCRIPTION UNITS UNITS - CA c~ Axial f o r c e c o e f f i c i e n t of forebody
CAP f 'AP
Axial f o r c e c o e f f i c i e n t of d e c e l e r a t o r
CC(3,3) C i j
Elements of t r ans fo rmat ion mat r ix from
i n e r t i a l coord ina tes t o body coord ina tes
of the f orebody
CCAP (8,8)
A t h r e e dimensional a r r a y of v a r i a b l e s
s i g n i f y i n g a x i a l f o r c e c o e f f i c i e n t s
of the forebody corresponding t o
AAM(8) , AALPFE ( 16) , and PPHIE (8)
A two dimensional a r r a y of 64 v a r i a b l e s
s i g n i f y i n g a x i a l f o r c e c o e f f i c i e n t s of t h e
d e c e l e r a t o r vFth r e s p e c t t o ang le of
a t t a c k corresponding t o AAMP(1) th ru
. .- AAMP (8)
B - . , CCHI c x c o s ( x
CCLL(8,16,8) A t h r e e dimensional a r r a y of v a r i a b l e s
s i g n i f y i n g r o l ' : ng moment c o e f f i c i e n t s
of the forebody c o r ~ e s p o n d i n g t o
AAM(8) , AALPME (16) , and PPHIE (8)
CCLLP (8,8,8) A t h r e e dimensional a r r a y of v a r i a b l e s
s i g n f f y i n g r o l l damping c o e f f i c i e n t s
of t h e forebody corresponding t o
AAMD(8) , AALPDE (8) , and PPHIDE (8)
A t h r e e dimensional a r r a y o f v a r i a b l e s
s i g n i f y i n g p i t c h i n g moment c o e f f i c i e n t s
of the f orebody corresponding t o
AAM(8) , AALPME (16) , and PPHIE (8)
GOODVEAR AEROSPACE CO.CO.4TI0U
GE R- 1604 7 *, . * I METRIC ENGLISH
FORTRAN STANDARD DESCRIPTION U N I T S U N I T S
C C W (8 , a , 8) A t h r e e dimensiona 1 a r r a y of v a r i a b l e s
t : s i g n i f y i n g p i t c h damping c o e f f i c i e n t s A,
of the forebody corresponding t o
CCLNR ( 8 , 5 , 8 )
CCNP (8,8)
AAMD (8) , AALPDE (8) , and PPHIDE f 8)
A t h r e e dimensional a r r a y of v a r i a b l e s
s i g n i f y i n g yawing moment c o e f f i c i e n t s
of the forebody corresponding t o
AAM(8) , AALPME (16) , and YPHIE (8)
A t h r e e dimensional a r r a y of v a r i a b l e s
s i g n i f y i n g yaw damping c o e f f i c i e n t s of
t h e foreoody corresponding t o
AAbi(8) , AALPDE (8) , and PPHIDE (8)
A t h r e e dimensional a r r a y o f v a r i a b l e s
s i g n i f y i n g normal f o r c e c o e f f i c i e n t s of
t h e forebody corresponding t o A A M ( ~ ) ,
AALPFE (16) , P P H I E (8)
A two dimensional a r r a y of 64 v a r i a b l e s
s i g n i f y i n g t h e p i t c h i n g moment c o e f f i c i e n t s
o f t h e d e c e l e r a t o r wi th r e s p e c t co ang le
of a t t a c k corresponding t o AAMP(1) t h r u
AAMP (8)
A two diinensjonal a r r a y of 64 var iaLles
s i g n i f y i n g t h e normal f o r c e c o e f f i c i e n t
of the d e c e l e r a t o r wi th r e s p e c t t o angle-
of a t t a c k corresponding t o AAMP(1) t h r u
AAMP (8)
Elements of t r ans fo rmat ion mat r ix from
i n e r t i a 1 coord ina tes t o body coord ina tes
of t h e d e c e l e r a t o r
GOODVEAR AEROSPACE CO.CO14110m
GER- 16047
METRIC ENGLISi FORTRAN STANDARD DESCRIPTION UNIT S UNITS - CCRIT c l c r = damping r a t i o :: 0.06
4
CCY(8,16,8) A t h r e e dimensional a r r a y of v a r i a b l e s
CDAP
CGAMP C 'fp
CHIE 7(
CHIPE X P
CLL C1l
CLLP
CLM C m
CLN 'n
CLNR
s i g n i f y i n g s i d e fo rce c o e f f i c i e n t s of
t h e forebody corresponding t o AAM(8),
AALPFE (16) , and PPHIE (8)
Drag a r e a o f d e c e l e r a t o r
Center of g r a v i t y
Cos ('6 )
cos ( '6
F l i g h t path angle o f forebody i n
h o r i z o n t a l p lane, measured from X
a x i s toward Y a x i s
F l i g h t path angle of d e c e l e r a t o r I n
h o r i z o n t a l p lane, measured from X
a x i s toward Y a x i s
Rol l ing moment c o e f f i c i e n t of the
f orebody
Rol l ing damping c o e f f i c i e n t of the
f orebody
P i t c h i n g moment c o e f f i c i e n t o f the
f orebody
P i t c h damping c o e f f i c i e n t of t h e
fo re body
Yawing moment c o e f f i c i e n t of the
f orebody
Yaw damping coefficient of t h e forebody
METRIC ENGLISH UNITS UNlTS - -I_--
FORTRAN STANDARD
CN C~
CNP C N P
COM(20)
DESCRIPTION - Normal f o r c e c o e f f i c i e n t of the forebody
Normal f o r c e c o e f f i c i e n t of t h e d e c e l e r a t o r
Input v a r i a b l e used t o d e f i n e computer
s imula t ion - up t o e i g h t y f i g u r e s
CPH I C b Cos (4)
CPHII cdi
CPHIPI Cbpi
CPSI c
CPSIP C y P
CS Cs
CSIGP
P
Damping c o e f f i c i e n t of t e t h e r
Cosine of one h a l f t h e apex angle of
t h e cone formed by the suspension l i n e s
CTHE CQ
CTHEP CQ P
CY CY D d
Cos (8 ) P
Side fo rce c o e f f i c i e n t o f forebody
Aerodynamic re fe rence l eng th of
forebody (m) f t
A two dimensional a r r a y of v a r i a b l e s
s i g n i f y i n g t h e c o e f f i c i e n t s of the (kg) s lug second d e r i v a t i v e s i n t h e equat ions of
01: , cnotions (kg- mL) s lug- f t
Z
-0 -a 7"' PI- -0 m u -
DDP (3,3) A two dimensional a r r a y of v a r i a b l e s (kg) s l u g s
s i g n i f y i n g t h e c o e f f i c i e n t s of the o r (kg-m2) s l u g - f t 2 second d e r i v a t i v e s i n t h e equat ions of
motion cf the d e c e l e r a t o r
T o t a l s t ~ s p e n s i o n l i n e d e f l . . t ion a r r a y m f t DELSX
GE R- 16047
HLTRIC ENGLISH FORTRAN STANDARD DESCRIPTION - UNES UNITS
DELTX Tc ta l t e t h e r l i n e de i l ec t i on a r r ay m f t
DLTO I n i t i a l elongat ion of t other beyond
DLTX
a ". 3 * I
DPR
DSP
'0 - Inw -
n -. DSY 1 -I. ' Lu I d E
unstretched length. DLTO is negative
i f the forebody and dece ldra tor con-
fluence points a r e c loser together
than LTO (m) f t
Tether de f l ec t i on component i n array
element (DLX(1) associated with load
PX(1) m f t
Ef fec t ive sp r ing def 1ectic.t a r r ay m f t
Rate of change of longi tudinal added
mass kglsec s luglsec
Aerodynamic reference length of
dece le ra tor (m) f t
Degress per radian - 57.2957795
Parachute diameter associated with SP m f t
Parachute projected diameter associated
witn DS m f t
Larger absc issa of two points cn the
longi tudinal added mass versus
Do log log p lo t m
Smaller abscissa o f two points on the
longi tudinal added mass versus
Do log log p lo t m
Larger abscissa of two points on the
l a t e r a l added mass versus
D3 log log p lo t m
COOOVEAR AEROSWCE CO.CO..~IO.
GER- 16047
METRIC ENGLISH FORTUS I
STANDARD DESCRIPTION UNITS UNITS - DSY 2 Smaller absc issa of two points on the
l a t e r a l added mass versus
Do log log p lo t
DT In t eg ra t i on increment
m1
DTT
DYPR q
DYPRP P
EE ( 6 )
EPL
EPS
EPSI
-I,. ' W W K EPT
Number of i n t eg ra t i ons betweek; da ta
output
Number of i n t eg ra t i ons between da ta
output when T b TDTC
f t
s e c
In t eg ra t i on increment when T < TDTC sec
Estimated parachute system period112 sec sec
Dynamic pressure of forebody (~lrn*) l b f / f t 2
Dynamic pressure of dece le ra tor (Il/m2) lb f 1 f t 2
An a r r ay s ign i fy ing the nonhomo- (Nlm) f t - l b f
geneous terms i n the s i x equations or
of motion of the forebody (N) Lb f
Suspension l i n e s t r a i n a r ray d m f t l f t
Number used t o check f o r incons,:;tent - 13
equations i n PIVERT Subroutine, !O
Number used t o check i f 8 i s approaching 2n+l
a s ingular point 8 = - 2 T. I f 0 is
approaching a s ingu la r po in t , t he
acce le ra t ions a r e kept fixed u n t i l t h i s
region is passed. EPSI = 0.00003bl
f reezes the acce le ra t ions i f 8 i s within
0 . 2 ~ of a s ingular point.
Tether l i ne s t r a i n a r r ay m/m f t l f t
COODVEAR comeonateom A- f i
GER- 16047 i
METRIC ENGLISH FORTRAN - STANDARD DESCRIPTION UNITS UNITS - ETA1 Number which con t ro l s D!I i f 0 i s near a L
s ingu la r point. ETAI = 0.00061 s e t s
DT = "'15, i f 0 is within 2 O of a
~. s i n g u l a r i t y
An a r r ay s ign i fy ing the acce le ra t ions (m/sec) f t / s e c o r
of the dece le ra tor rad/sec
Estimated psrachute sys tem frequer :y l / s e c I / sec $ 3 FREQP
Ultimate design f a c t o r of s a f e t y f o r
parachute
Generalized force on forebody i n
X d i r e c t i o n
Body force i n d i r e c t i o n o f 1 6 due t o
aerodynamics
FXP Generalized force on dece le ra tor i n
X d i r ec t i on
Body force i n d i r ec t i on of X pb
Generalized force on forebody i n
Y d i r e c t i o n
Body force i n d i r ec t i on of Y dce b
t o aerodynamics
FYB
FYP G a e r a l i zed force on dece le ra tor i n
Y d i r e c t i o n
FYPB
FZ
Body force i n d i r ec t i on of Y P"
Generalized force on forebody i n
Z d i r e c t i o n
FZB Body force i n d i r ec t i on of Zb due t o
aerodynamics
7- . . . , - . . .
.. .- .-I- f $?
,L f COOOVEAR AEROSPKE
C O I C O . 4 l l O I
GER- 16047 ;*. .
h fl METRIC ENGLISH FO-TRAN STANDARD DESCRIPTION UNITS UNITS
IXPB I Apparent moment o f i n e r t i a about xpb
X a x i s 2 (kg-m) s l u g - f t
2 pb
IXY B I Product of i n e r t i a a s s o c i a t e d w i t h xyb
X,, and Y axes (kg-m2) s lug- f t 2 b
IXZB 'xzb Product of i n e r t i a a s s o c i a t e d wi th
X,, and Z axes 2 (kg-m ) s lug-f t
2 b . ..
.R, k: IYB I Moment of i n e r t i a about Y a x i s 2
(kg-m ! s lug- f t 2
#- yb b
IYPB I Apparent moment of inertia about Y P ~
Y a x i s 2 (kg-m ) s l u g - f t 2
pb
IYZB I yzb
Product of i n e r t i a a s s o c i a t e d wi th t Yb and Z axes 2 a (kg-m j s lug- f t
2
ST b ir
IZB Moment of i n e r t i a about Zb a x i s 2 (kg-m) s l u g - f t 2
d I z b
KS
LS
LSCL
LTD
LTO
M
MP
MPAL
Ks Te ther s p r i n g cons tan t
Suspension l i n e l eng th
Dis tance a long parachute c e n t e r l i n e
between t h e confluence p o i n t and the
p ro jec ted diameter p lane
5 Tether l eng th - d i s t a n c e between
confluence po in t s
i Time r a t e o f change o f t e t h e r length
5 0 Unstretched t e t h e r length
m Mass of forebody
m Real mass of d e c e l e r a t o r P
Added mass o f t h e d e c e l e r a t o r a long X pb
axis
(m) f t
(mlsec) f t / s e c
(m) f t
(kg) s l u g s
(kg) s l u g s
GER- 16047
METRIC ENGLISH FORTRAN STANDARD DESCRIPTION UNITS UNITS - WAS Added mass of the d e c e l e r a t o r a long Y
~b o r Z a x i s
P" (kg) s l u g s
MPL m Apparent l o n g i t u d i n a l (X ) mass o f PI P"
d e c e l e r a t o r (kg) s l u g s
r MPS
-- . 3 ' . OHYBE
0mBE /
OPAM
OPDA
m Apparent s i d e (Y o r 2 ) mass o f P s pb pb
d e c e l e r a t o r (kg) s l u g s
Number of parachute suspension l i n e s
Nuruber o f t e t h e r l i n e s
Option v a r i a b l e : i f OMETRC = l., Inpu t
and Output a r e i n t h e m e t r i c system.
I f OMETRC = 0.0 Inpu t and Output ar.2 i n
t h e Engl ish system.
Wxb Angular v e l o c i t y about Xb a x i s
y ~ b Angular v e l o c i t y about Yb a x i s
Wzb Angular v e l o c i t y about Zb a x i s
Option v a r i a b l e : i f OPAM = l . , added
mass of t h e d e c e l e r a t o r f 0; i f
OPAM = O., added mass of d e c e l e r a t o r
= 0
Option v a r i a b l e : i f OPDA = l . , damping
moment c o e f f i c i e n t s of t h e forebody a r e
read i n a s a r r a y s ; i f OPDA = U, damping
moment c o e f f i c i e n t s a r e read i n a s
cons tan t s
Option v a r i a b l e : i f OPOS = l., a t l e a s t
one of t h e c.g. o f f s e t s o r products o f
i n e r t i a o f t h e forebody f 0.; i f OPOS = O.,
a l l c.g. o f f s e t s and products of i n e r t i a = 0.
deg l sec
deg l sec
deg l sec
WOOVEACl A E R O W CO.CO.A1IOm
GER- 16047
METRIC ENGLISH UNITS UNITS - FORTRAN STANDARD DESCRIPTION
Option v a r i a b l e : i f 0PPZI)T = I . , a p l o t
t ape can be made; i f 0PPII)T =@., no p l o t
t ape i s made.
OPPRIN Option v a r i a b l e : i f OPPRIN = l., a l l aero-
dynamic c o e f f i c i e n t a r r a y s a r e p r i n t e d ou t ;
i f OPPRIN = Q, no aerodynamic c o e f f i c i e n t
a r r a y s a r e p r i n t e d o u t
OPDT Option f o r automat ic DT determinat ion
(OPDT = 1)
Option f o r automatic parachute a r e a
c a l c u l a t i o n s (OPSP = lJ
OP SP
OPSYM Option v a r i a b l e : i f OPSYM = l., t h e f o r e -
body is aerodynamically symmetric such t h a t
Cy% =Cm=O; i f OPSm = 0, t h e forebody i s
no t symmetric
Parachute o v e r i n f l a t i o n a t r ee fed
s t a g e ( I ) . (percent/100)
PCTO 1 Parachute over in f l a t i o n a t reefed
s t a g e 1. (percent/100)
Parachute o v e r i n f l a t ion a t r ee fed
s t a g e 2. (percent/100)
Parachute over inf l a t i o n a t reefed
s t a g e 3. (percent / 100)
PHIAE Aerodynamic r o l l ang le of forebody,
0 S PHIAE \( 180' . a
PHIDDE 0 Angular a c c e l e r a t i o n about Xb a x i s . PHIDE 8 Angular v e l o c i t y about X a x i s b
GER- 16047 - - ~~
METRIC ENGLISH FORTRAN STANDARD DESCRIPTION - UNITS UNITS
PHIE 0 Euler angle ro t a t i on about % ax i s deg
PHIIE Ii Aerodynamic r o l l angle o f forebody,
-180' * tii ,( 180" deg
PHIPI 8 Aerodynamic r o l l angle of dece le ra tor , o i
POROS Parachute porosity. Use POROS = 0.15
PPHIDE (8 ) An a r r ay of e igh t var iab les s ign i fy ing
forebody r o l l angle used with damping
coe f f i c i en t s deg
PPHIE (8) An ar ray of e igh t var iab les s ign i fy ing
forebody r o l l angle used with fc.~*ce
and moment coe f f i c i en t s de8
PSIDDE Angular acce l e r a t i on of forebody
about -2 a x i s deglsec 2
PSIDE Angr?!sr -:c?eclty of fo i~bodj i about
-Z ax i s deglsec
PSIE Euler angle r o t a t i o n of forebody about
-Z a x i s de8
PSIPDE Angula r ve loc i ty of dece le ra tor about
-Z axis
Angular r o t a t i o n of dece le ra tor about
-2 ax i s
PSIPE
Angular acce l e r a t i on of dece le ra tor about
-2 axis deg/sec 2
PSPDDE
PULAN Angle between t e t h e r and forebody
cen t e r l i ne del3 de8
Suspension l i n e load a r r ay N l b f
GER- 16047 A. -4 - '* *
FORTRAN STANDARD DESCRIPTION METRIC ENGLISH UNITS UNITS -
I PSX
m
Tota l suspension l i n e load a r r ay N lb f
Tether l i n e load a r r ay N lb f
PTX N d
Tota l t e t h e r l i n e a r ray lb f
x a Ef f cc t i ve spr ing load array N
lbf
Parachute load due t o r a t e of change
of mass of the parachute t i m e s the
r e l a t i v e ve loc i ty , XPBDR N lb f
QPHI Q0
Generalized force about X,, a x i s (m-N) f t - l b f
QPSI Q Y Generalized force of forebody about
-2 ax i s (m-N) f t - l b f
QPSIP Q ~ P
Generalized f o r of dece le ra tor about
-2 a x i s fm/N) f t / l b f
QTHE Q~ Generalized force of ,orebody about
negative l i n e of modes (m/N) f t / l b f
QTHEp Generalized force of dece le ra tor about Qe, negative l i n e of modes
RATIO Nondimensional length used i n the
dece le ra tor ' s body torque expressions
(m) f t Radius of e a r t h - 20,926,435. :A
RHO Density of atmosphere a t 2 (1962
Standard)
A i r dens i ty r a t i o (RHO/RHOO)
Aerodynamic reference a r ea of forebody
sin ( X
s i n ( Xp)
s i n ( X 7 - SGAM
GER- 16047 :I ! a ! f - METRIC ENGLISH 4
\ / FORTRAN STANDARD DESCRIPTION UNITS UNITS 3 - 2
SP
SPD
Parachute drag a r ea
Time r a t e of change of parachute
drag area 2 2 3
m I s ec f t I s ec
SPHI S0 s i n (0)
SPHII sOi
SPHIPI SQpi
s i n (gi)
s i n (0 .) P l
I n i t i a l parachute area m 2
f t 2
F i r s t reefed s tage parachute drag a rea 2 f t 2
Second reefed s tage parachute drag a rea m 2 f t 2
Third reefed s tage parachute drag a rea m 2 f t 2
SPRQ
SPRl
SPRL Parachute drag a rea associated with
reefed s t age (1-1) m 2
f t 2
SPRU Parachute drag area associated with
reefed s t age ( I ) m 2 f t 2
SPSI s 'I SPSIP st',
SSP (16) An ar ray of s ix t een var iab les s ign i fy ing
aerodynamic reference a rea of the
decelerator corresponding t o TTIP (16) (2) f t 2
STHE SO s i n (0) - 0 - 40
7" STHEP SO so m w
P s i n (0 )
P - 0 .. r T F l igh t time - LL
s e c
lbf
I W UP
TENS Tension i n t e t h e r
-c THEDDE '0' Angular acce l e r a t i on of forebody about
negative l i n e of nodes
COOWEAR A€- COmCOm.TIO*
GER- 16047
METRIC ENGLISH UNITS UNITS - FORTRILV STANDARD DESCRIPTION
THEDE Angular v e l o c i t y of forebody about
negat ive l i n e of nodes deglsec 2
THEE Euler angle r o t a t i o n of forebody about
negat ive l i n e o f nodes deg
deglsec
deg
deglsec 2
Angular v e l o c i t y of d e c e l e r a t o r about
nega t ive l i n e of nodes
THEPDE
THEPE Euler angle r o t a t i o n of d e c e l e r a t o r about
negat ive l i n e of nodes
Angular a c c e l e r a t i o n of d e c e l e r a t o r about
r ~ e g a t i v e l i n e of nodes
TPPDDE
TFI Parachute i n f l ~ t i o n time from Stage ( I )
t o Stage ( I + 1) s e c
TINF Time when i n f l a t e d a rea f i r s t
equa l s SPRU s e c s e c
Tin.? a t s t a r t of i t f l a t i o n of
I cefed s t a g e ( I )
T o t a l Time spen t i n region where
(1-ABS,(SIN(THE)). LT .EYSI) s e c
TNINY
s e c
s e c I n i t i a l t ime s e c
One ' a ? f t h e time spen t i n t h e o w r -
i n f l a t i o n of s t a g e I scc s e c
Time a t s t a r t of i n f l a t i o n o f f i r s t
reefed s t a g e
TOTRO
s e c s e c
s e c s e c
s e c s e c
TOTRl Time a t s t a r t o f i n f l a t i o n of
second reefed s t a g e
Time a t s t a r t of i n f l a t i o n of
t h i r d reefed s t a g e
TPD T o t a l t e t h e r load ( t e n s i o n + damping) N lb,
GOOUVEAR AERO- C O n C O n A I I O I
GER- 16347
TPDRB
1
TPDXB
TPDYB
TPDZS
" ". m
I -, METRIC ENGLISH 1
FORTRAN STANDARD DESCRIPTION - UNITS UNITS ! 1
TOTR3 Time a t end . t h i r d ree fed s t a g e s e c s e c , I ,IC
TRO
TTIP (16)
TXB Txb
TTT
TYB
TYPB
TZB
Component of t e t h e r load normal t o
f orebody c e n t e r l i n e N l b f
Tether load component along f orebody
XI3 a x i s N l b f
Tether load component a long forebody N lb f
YB a x i s
Te ther lozd component along forebody
ZB a x i s N lb f
Time from "TO" t o s t a r t of i n £ l a t i o n
of f i r s t s t a g e s e c s e c
Time from "TO" t o s t a r t of i n f l a t i o n
of second s t a g e s e c s e c
Time from "TO" t o s t a r t of i n f l a t i o n
of t h i r d s t a g e s e c s e c
Time from "TO" t o end of t h i r d
s t a g e (TR3 ) TTT) s e c s e c
An a r r a y of s i x t e e n v a r i a b l e s s i g n i f y i n g
i n £ l a t i o n time even t s s e c
T o r q u e a b o u t X a x i s d u e toaerodynamics (m-N) i t - l b b f
F l i g h t time a t which s imula t ion i s ended s e c
Torque about Yb a x i s due t o aerodynamics (m-N) f t - l b f
Torque about X a x i s due t o aero- pb
dynamics (m-N) i t - l b f
Torque about Zb a x i s due t o aerodynamics (M-N) f t - l b f
GER- 16047 . -
t t r $ 8
FORTRAN STANDARD DESCRIPTION METRIC ENGLISH UNITS UNITS --
4 4 i TZPB T Torque about Z a x i s due t o aero-
zpb P" dynami cs (m-N) f t - l b f
v
VP
VS
WT
WTC
WTCM
WTL
WTLM
WTP
X
Tota l ve loc i ty of forebody (mlsec) f t l s e c
f t l s e ~
f t l s e c
l b
l b
s lug
l b
s lug
l b
Tota l v e l c c i t y of dece le ra tor (mlsec)
Speed of sound a t Z (mlsec)
Weight of forebody (N)
Parachute canopy weight i;,o*m/sec 2
Mass of parachute canopy kg
Parachute suspension l i ne s weight kg*rr./sec 2
Mass of parachute suspension l i n e s kg
Weight of dece le ra tor (N)
Down range i n e r t i a l ax i s or d i sp l acemnt
of forebody (m)
Longitudinal body ax i s o r displace-
ment of forebody (m)
c.g. o f f s e t alocg 5 a x i s (m)
Xb body ax i s ve loc i ty (mlsec)
XBAR
XBD
Down ranbz ve loc i ty of forebody (mlsec; X,
XDD
XP
Down range acce le ra t ion of forebody 2
(mlsec )
Down range displacement of
dece le ra tor
X body ax i s v e l o r 4 t y of dece le ra tor (mlsec) P"
XPBD
XPBDR Velocity of a i r en t e r ing or ex i t i ng
the parachute r e l a t i v e t o the
parachute ve loc i ty rnlsec
COOWEAR AEROSPACE C O . C O ~ . l I O I
GER- 16067
METRIC ENGLISB UNITS - UNITS STANDARD DESGRIPTION
XPBDI Para-hute ve loc i ty , XPBD, a t T TINT.
It i s used t o ca lcu la te f i l l time m/sec f t /sec
XPD 4 Down range veloci ty of dece lera tor (m/sec) P
f t / s e c . . XPDD X Down range acc , l e r a t i o n of decelerqtor (m/sec ) f t l s e c 2 2
P
Cross range i n e r t i a l ax is o r displace-
ment of forebody (m)
Letera l body ax i s or displacement of
f orebody (m)
I - YBAR Y
5 v
c.g. o f f s e t along Yb ax is (m)
1' YBD Yb
t ! YD ; . . YDD Y
Yb body ax is ve loc i ty (m/sec)
Cross range ve loc i ty of forebody (m/sec) n
Cross range acce lera t ion of forebody (m/sec )
Cross range i n e r t i a l displacement of
decelerator (m)
YPBD Y body ax is veloci ty of decelerator (m/sec) P b
YPD 4 P
Cross range ve loc i ty of decelerator (m/sec)
Z Cross range acce lera t ion of dece lera tor ( m / ~ ? c ) YPDD
z Vert ica l i n e r t i a l ax is o r displacement
of forebody (m)
Ver t ica l body ax i s or displacement of
forebody (m)
GER- 16047
METRIC ENGLISH FORTRAN STANDARi) DESCRIPTION - 'UNITS UNITS -
i ZBAR 2 c.g. o f f se t along Zb axis (m) f t $ 4
!
ZBD Zb Zb body axis veloci ty
Vert ical ~ ~ 1 0 1 : i t y i ~ f forebody (mlsec) f t / sec 0. 2 2
ZDD 2 Vertical accelerat ion of forebody (m/sec ) ft./sec
Z P 2 Vert ical i n e r t i a l displacement of P
decelerator (m) f t
ZPBD 2 Z body axis veloci ty of decelerator (mlsec) f t / s e c pb P"
t P D Z Vert ical veloci ty of decelerator (m/sec) f t / sec . P
2 2
ZPDD Vert ical accelerat ion of decelerator (m/sec ) f t l s e c 2
P
CHAPTER I - INTRODUCTION
The system t o be simulated is two r i g i d bodies joined by an e l a s t i c t e the r . The
forebody may have a completely general shape and mass charac ter is t ics , and w i l l be
free t o move with s i x degrees of freedom (three t rans la t ional , three ro ta t ional ) .
The decelerator is assumed t o be symmetric i n shape and mass cha rac te r i s t i c s about
its longitxdina: ( r o l l ) sx f s , and sill be f r e e t o move with ffve degrees of freedom
(three trscslat iona ' l , two ro ta t ional ) . A f r i c t ion less swipe1 is assumed a t the
decelerator-tether confluence poict . Thus the r o l l motions of t h e forebody w i l l not
couple with the decel.erator. The t e the r i s simulated by a spring and dashpot i n
para l le l . Damping coefficfents fo r t e the r l i n e s are d i f f f c u i t t o obtain; but spring
constants f o r a te ther cqn be f o ~ i d fro^ experimental s t r e s s s t r a i n curves. Conse-
quently; the damping coeTficient is a s s u e d constant, while the spring constant is
assumed t o be a functfon of clangs-Lion i n tile computer program, thereby introducing
2 qunsilinear spring.
CHAPTER I1 DERIVATION OF EQUATIONS OF MOTION
SECTION 1 - COOPJ)INATE SYSTEM
Figure 1 shows t h e d i f f e r e n t coordinate s y s t e r s used t o de r ive t h e
equations of motion. X y Z is an i n e r t i a l orthogonal coordinate
system at tached t o a f l a t non-rot;.ting ea r th . XYZ and X y 2 a r e P P P
orthogonal axes f i xed t o t he forebody and dece l e ra t e r a t "On and
"0 " respect ively . Coardinate systems XYZ and X Y 2 t r a n s l a t e P P P P
with the bodies but do not r o t a t e , always remaining p a r a l l e l t o
corresponding i n e r t i a l axes. The displacements X, Y , 2 , Xp, Y P'
and 2 a s measured Prom t h e o r i g i n of HY 2 , a r e t h e s i x P'
t r a n s l a t i o n a l degrees of freedom of t h e two bodies. The reference
forebody body axes, longi tud ina l (Xb) , l a t e r a l ( Y b ) , and v e r t i c a l
(2 ) , i n t e r s e c t a t "On, the o r i g i n of t he aerodynamics load system b of t h e forebody. The re fe rence dece le ra tor body axes, longi tud ina l
(X ) , l a t e r a l (Ypb) , and v e r t i c a l ( 2 ) intersect a t '0 ", t he pb pb - P
c.g. of t h e dece le ra tor , The va r i ab l e s X , P, a r e t he d i s tances
frozr "a" tu the c o g . of t h e forebody measured p o s i t i v e l y i n t h e
d i r e c t i o n of the p o s i t i v e body axes Xb8 Yb, Zb r espec t ive ly . For
o r i e n t a t i o n purposes, t h e reader should pos i t i on himself a s a
p i l o t i n an a i rplane. In t h i s pos i t i on , Xb is p o s i t i v e toward t h e
nose, Yb is p o s i t i v e toward the l e f t wing and Zb is p o s i t i v e up. +
is the vector d i s tance from t h e i n t e r s e c t i o n of t h e longi tud ina l ,
l a t e r a l , and v e r t i c a l axes of the Lorebody ('On) t o t he t e t h e r +
confluence point of the forebody. r2 i s t h e vector d i s tance from
t h e c.4. of t he dece le ra tor ("0 " 1 t o the t e t h e r confluence po in t -+ P
of the dece le ra tor ; r2 lies aloqg X pb'
8 p
COOOVEAR AEROSPWE CO.~O.A~lO.
GER-16047 --.
f 3 . , SECTION 2 - EULER ANGLE TRANSFORMATION
I n order t o spec i fy t he angular o r i e n t a t i o n of a body w i t h r e fe rence
t o a non-rotat ing coordinate system (X, Y , Z ) , t h r e e succesr'-r-
r o t a t i o n s a r e made a s shown i n Figure 2. The f i r s t r o t a t i o n i s . i n
t h e d i r ec t ion , -02, such t h a t OX and OY a r e r o t a t e d tnrouy. an angle
$ i n t o Oa and ON respec-ively. The second r o t a t i o n is i n t h e I d i r e c t i o n , -ON, such t h a t Oa and OZ are ~ o t a t e d t h ~ o u g h an angle 0
i n t o OXb and Ob respec t ive ly . The finl: r o t a t i o n i s about O X , i : 3 .
such t h a t ON and Ob a r e ro t a t ed through an angle 4 i n t o OYb and 3 # -
i- s OZb respect ively . The t h r e e angular ro ta t io l i s ($, 0 , 4 ) spec i fy [: \
t h e o r i e n t a t i o n of t h e body axes (Xb , Yb, Z b ) with respec t t o t h e e F 5 -
i n e r t i a l axes ( X , y , ). Again, from a p i l o t s viewpoint, a 1 'r
[. t p o s i t i v e is a nose t o t he r i g h t yaw; a p o s i t i v e 0 is a nose up p i t ch ; ?' [ and a p o s i t i v e @ is a r i g h t wing down r o l l . g ,--- I
. i
The transformation matrix between the body axes and i n e r t i a l axes
is now found by considering one r o t a t i o n a t a time and then combining. ,~ -
The f i r s t r o t a t i o n i s given by:
where S$ = s i n $ and C$ = cos$.
' -0 The second r o t a t i o n is: - Q
, I.0
a
3
The f i n a l r o t a t i o n is: a---
By s u b s t i t u t i n g Equation (1) i n t o ( 2 ) and (2 ) i n t o ( 3 ) , t h e
transformation matrix [Cl is formed
- Since [CI i s a l i n e a r orthogonal ( 2 Cij Cik - 6 jk: j , k = 1.2.3) i=l
transformation, i t s inverse i s equal t o i t s transpose. . Therefore . cqce -c$ses$+s+c$
s$ses$+c$c$
ces4 cec4
For t h e dece le ra tor , t h e r e i s no r o t a t i o n about tho longi tud ina l
axis. Consequer.tly, t h e transformation matr ix i n ( 4 ) i s s imp l i f i ed
by 1e;lting $ = 0, The r e s u l t i s [C 1. P
The t o t a l angular v e l o c i t y of the forebody i s given by: + .+
( 7 )
1 From the inverse of ( 3 ) : i
t i
From the inverse o f (2'
i; = se?,+ cet,
subst i tut ing ( 8 ) , (91 , and (10) i n t o (7) :
The components of angular v e l o c i t y for the forebody are:
Likewise, f o r the decelerator , the angular v e l o c i t i e s are:
SECTION 3 - KINETIC ENERGY
The k i n e t i c energy of t h e system is due t o t h e t r a n s l a t i o n a l and
r o t a t i o n a l v e l o c i t i e s of t he forebody and t h e dece le ra tor . The
forebody is completely general i n shape, and products of i n e r t i a
and c.g. o f f s e t s w i l l e f f e c t t h e k i n e t i c energy. On t h e o t h e r hand, the dece le ra tor is assumed t o be symmetric about t h e longi tud ina l
a x i s and t h e aerodynamic loads a r e referenced t o t h e cog. Therefore,
a l l products of i n e r t i a and cog. o f f s e t s a r e zero. The expression
f o r k i n e t i c energy is: ( 2 )
I n Equation (18) , mp fi and m include d i r e c t i o n a l m a s s terms due i PS >
i ,
ts t he a i r enclosed i n t h e canopy, Ixpb, I ypb, and 1 zpb are apparent mass moments of i n e r t i a .
SECTION 4 - POTENTIAL ENERGY
The p o t e n t i a i energy of t h e system i s due t o t h e g r a v i t a t i o n a l
p o t e n t i a l of both bodies and t h e e l a s t i c p o t e n t i a l o f t h e t e t h e r .
L ~ o is t h e unst re tched length of t h e t e t h e r ; and LT i s the
s t r e t ched length of t h e t e t h e r a s given by t h e geometry of t h e
system. Referring t o Figure 1:
Sl and d2 a r e t h e vec tors from t h e i n e r t i a l coordinate system
( X 7 2 ) t o t he confluence po in t s of t h e forebody and
dece le ra tor respect ively . For t h e forebody,
a, b, and c a r e measured along p o s i t i v e body axes Xb8 Yb8
and Zb r espec t ive ly . Using t h e coordinate transformation
matrix ( 4 ) :
Similarly f o r the decelerator:
Substituting (25) i n t o (24) and using matrix Equation ( 6 ) , I
1 j
Define the variables Z, E, and F such that: i
Further on i n the derivat ion it w i l l be necessary t o know 7
the t o t a l t i m e der ivat ive o f L, and the p a r t i a l der ivat ives 5 3
.L
o f x, g, and F with respect t o the generalized coordinates. 1
SECTION 5 - RALEI'3I'S DISSIPATION FUNCTION
If the viscous damping force is proportional to the velocity
of the particle at which the force acts, an expression analogous
to the potential energy of a spring may be used. This fufiction,
F, 2s known as Rayleigh's di~sipation function, and is defined as (1 1
For this problem, Rayleigh's dampirt : r :maidered only in
the tether.
SECTION 6 - LAGRANGE'S EQUATION L he Lagrange equation for a non-conservative (aerodynamic forces)
system with a holonomic (can be expressed as an algebraic
c% expression), scleronomous (independent cf time) constraint and b- -0 S W
Rayleigh's dissipation function (dampir-g in the elastic tether) I 0 .. -LL can be written as:
(1 1 , ;Y
I n Equation (54), t h e t e r m X 2: expresses t h e general ized A
force exer ted by t he t e t h e r on t h e ith degree of freedom.
The zons t r a in t equation is:
Qi is t n e general ized fo rce due t o t h e aerodynamics.
is t h e fo rce due t o damping i n t h e tether. Ti Tne Lagrangian is equal t o *.e t o t a l k i n e t i c energy of t h e
system (Equation (18)) ini3us t h e t o t a l p o t e n t i a l energy of
t h e system (Equation (19)) . With s u b s t i t u t i o n s from Equations
( 4 ) , (6) and (12) t o (17), t h e Lagrangian is:
-{I,,, r - i tces4)-i ( ~ $ 1 I r-j, (ceco)+i ( ~ 4 ) l+rXzb 1-i (se1+611-$ (c~c+)+; (s+)I
f . ~ .. Note: I Y P ~
= I due t o decelerator symmetry. '% zpb
SECTION 7 - GENERAL EQUATIONS OF MOTION ".
I . Equation ( 5 6 ) displays a l l of the generalized coordinates - i
e x p l i c i t l y except those appearing r? IT. The terms t o be
1 subst i tuted i n t o Equation (54 ) are now developed, 5 *
X Equation
3
i d a~ .. aF = x{rnl+${m [P(-s~ce)+P(slysesg+clycg~+Z(slysec~-c~s~ I 1
I ax
@ Equation
+" ''xb (-ce) + [ I ~ ~ ~ (S~P +I,,, (cm 1 (-s e 1 1
+$* #-
IIXyb (CQ)-lxzb (SO) I (Ce) }
. . LT Equation
k z Substituting Equations (101) t o (104) i r t o (54) y i e lds :
? = -[Ks($-LTO)+csiT~ (105)
f
i e The value o f A i n Equation (54) i s now defined and is expressed r !- i n "srms o f the e leven generalized coordinator and the i r t i m e f
5 der. datives. The eleven simultaneous, nonlinear,
1 coupled d i f f e r e n t i a l equations o f motion are then writ ten as:
$ 1. X Equation
I
2. Y Equation
Y {m)+T{-m [ Z C ~ ~ + % ? C ~ ~ + T C ~ ~ I 1
+'B'{m [X(s$se)+Y(s$ces~) +~(s)cecg) ] 1+${n [ F c ~ ~ - Z ~ ~ ~ 1 1
= $4 {-2m[X[c+se)+Y(c+ceso) +T(c$cecm) ] ~ + $ ~ { ~ : ~ [ ! T C ~ ~ - ' Z C ~ ~ ] 1
+ B ~ I - ~ ~ [ P ( s $ c ~ c ~ ) -H(sqces@) I ~+$2~m[~~c,2+~~22+'Zcc,,~ 1
+B2{m [~c~~+~(s~ses~)+T(s)secg) 1 1
3. Z Equation
. ~ { m } + i i { m [y(ce)+~(-seso+~(-secg) I ~ + i ' ~ m I ~ c ~ ~ - % ~ ~ I
= 8){2m[Y(~ec~)-T(sesg) 1 I + B ~ { ~ [ X C ~ ~ + ! T C ~ ~ + Z C ~ ~ I 1
8. 0 Equation f
Y{m [F(-c$s~) +F(-c$c~s$) +Z(-C~CCIC~) 1 I
- Q 9 . I$ Equation -
P;" tz Inw - P { ~ ~ ~ Y E ~ ~ - R ~ ~ ] }+f {m [ ! F c ~ ~ - Z C ~ ~ ] }+i'{m [iiC jr-Z~23 ] 1 0 ; td
W K . . +${-1xb~13+~xyb~23+~xzb Cj3}+B {Ixyb (CQ1-I (SO) 1+;{1&}
-C xzb
, ,-
-33-
6 F
-0 - m *el
PI- -0 Y)w - 0 .- -I. ' W
w a - . .
SECTION 8 - SIMPLIFIED EQUATIONS Or MOTION
muations (106) to (116) are for the most general situation
possible, and as a result, are quite lengthy. Under some
circumstances, these equations can be simplified. If this can.
be accomplished, a significant decrease in computer time will
be realized. The first simplification occurs if the forebody's
aerodynamic (body) reference axes are principal axes. In this
case = P = 'i = 0 and I - xl'b
= I - Ixzb yzb = 0. Equations (1061,
(107), (1081, (1121, (1131, and (114) become:
The above s i x equations of motion have not only been shortened,
bu t they a l s o have been uncoupled i n t he t r a n s l a t i o n a l
acce le ra t ions making them e a s i e r t o solve. The second
s imp l i f i ca t ion involves t h e dece l e ra to r degrees of freedom.
I f t he added masses of t he dece l e ra t c r a r e ignored (mpe= mps
= m p ) , E q u a t i ~ n s (109). ( 1 1 0 ) . (ill), (115). and (116) become:
Like t h e s i m p l i f i e d equat ion f o r t h e forebody, t h e d e c e l e r a t o r
equat ions have a l s o shortened. Furthermore, they have
completely uncoupled i n t h e second d e r i v a t i v e s making numerical
i n t e g r a t i o n easy.
SECTION 9 - GENERALIZED FORCES - AERODYNAMICS
The nonconservat ive f o r c e s a c t i n g on t h e forebody a r e due t o
aerodynamics. The aerodynamics and t h e convention used i n t h i s
r e p o r t apply t o t h e Space S h u t t l e S o l i d Rocket Booster (S.R.B..).
If a d i f f e r e n t body i s t o be s imula ted , t h e aerodynamic c o e f f i c i e 3 t s
and p o s s i b l y t h e convention used t o d e f i n e them, would change.
For t h e S.R.B., t h e aerodynamics a r e a f u n c t i o n o f r o l l angle ,
angle-of-at tack, and Mach number. The angle-of-a t tack i s
measured from t h e t o t a l v e l o c i t y v e c t o r t o t h e p o s i t i v e l o n g i t u d i n a l
a x i s (Xb) a s shown i n F igure 3 .
The normal f o r c e c o e f f i c i e n t , CN, i s i n t h e p lane formed by
t h e v e l o c i t y v e c t o r and t h e l o n g i t u d i n a l a x i s , and is perpend icu la r
t o t h e l o n g i t u d i n a l a x i s ( X b ) . The r o l l ang le , i s then
measured from t h e normal f o r c e c o e f f i c i e n t t o t h e Zbbody a x i s .
The a x i a l f o r c e c o e f f i c i e n t i s de f ined a s u s ~ . a l , p o s i t i v e i n
t h e nega t ive Xb d i r e c t i c n . F i n a l l y , t h e s i d e f o r c e c o e f f i c i e n t
i s perpendicular t o t h e Xb body a x i s and t o t h e normal f o r c e ,
such t h a t t h e d i r e c t i o n s o f CA, CN, Cy form a right-handed
or thogonal coord ina te system. Mathematically, t h e aercdynamics
r o l l ang le i s given by:
-1 Oi = t a n [-Pb/-ib]
i
GOODVEAR AEROSPACE CO.COn.l IO~
GER-16047
The positive directions of the moment coefficients are shown
in Figure 3 as double arrows. Dampirig . z n t coefficients
are about the body axes (Xb, Yb, Zb). Aerodynamic body axes
forces are given as:
The body axes forccs are converted to inertial a;?s force
using the elzments of LC], Equation 4.
-S are: Body axes torqu-
The body axis torques are transformed to generalized
torques using Equations (5) and (1) and iAotlng s i gn
conventions.
GOOWEAR AEROSPACE C O I . O I . l I C .
GER-16047
The aerodynamics of 3 dece l e ra to r (parachute) a r e not veil
known because a parachute i s not a r i g i d body, ana does not
lend i t s e l f t o e a s i l y obta inable t e s t daLa, espec iaJ ly
under dyna;.ic condit ions. Consequently the aerodynamics
of a symmetric c 'ocelerator tend tc be re la t ive ly simple
due t o a lack of bet'er understanding r a t h e r than the
i n a b i l i t y t o use ava i l ab ' t ineormation. I f b e t t e r
aerodynamic da t a i s a t t a i n a b l e , it i s a simple matter t c j
a l t e r the body forces and torques appropr ia te ly .
For t h i s r epo r t , t he decelerator body forces and t ~ r q ~ e e are :
The general ized iorces f o r the dece lerator are:
F xp
= F + F C xpb 'pll ypb p21 + Fzpb C p31
F YP
= F xpb 'pl2 + F ypb $22 + Fzpb $32
F z P
= F xpb Cp13 + 'zpb 'p33 - mpg
- Q * ~ - - Tzpb 'p33
= -T Y P ~
GER- 16047
CHAPTER 111 COMPUTER PROGRAM
SECTION 1 - FEATURES OF THE COMPFCER PRIGRAM
The computer program con ta ins t h e fo l lowing f e a t u r e s .
1. The progran has many opt ions which s i m p l i f y t h e inpu t of d a t a o r dec rease the
program run time. Use o f the op t ions a r e conta ined i n the f i s t i n g of t h e
program a s comment ca rds . These op t ions a r e :
a. Opt i )ns a r e included which change t h e dimensions of t h e aerodynamic
c o e f f i c i e n t a r r a y s a s d i c t a t e d by inpu t requirements.
b. An op t ion is provided (OPDT = 1.0) which au tomat ica l ly determines t h e
magnitude of t h e i n t e g r a t i o n t i m e i n t e r v a l , DT.
c. An op t ion is provided (OPSP = 1.0) which c a l c u l a t e s t h e parachute d rag
a r e a (SP) time h i s t o r y . I f t h i s op t ion is r o t used the d rag a r e a versus
time i s inpu t i n t o the program i n t h e form of look-up a r r a y s .
d. An o p t i o n f o r including l o n g i t u d i n a l and l a t e r a l added a i r mass e f f e c t s
on t h e parachute (OPAM = 1.0) i s included i n t h e program.
e . A p rov i s ion i s made t o use s i m p l i f i e d equa t ions of motion (OPOS = 0.0) t o
reduce run t ime, i f a l l t h e forebody products of i n e r t i a and c e n t e r 3f
mass o f f s e t s a r e equal t o zero.
f . An o p t i n ( o P P L ~ = 1.0) f o r making a p l o t tape is a v a i l a b l e .
g. English o r me t r i c systems may be used f o r d a t a inpu t and out by equat ing
OMETRC t o 0.0 o r 1.0 r e s p e c t i v e l y .
2. A l l aerodynamic c o e f f i c i e n t s a r e read i n t o the program a s func t ions of angle
of a t t a c k , r o l l ang le , and mach number i n t h e form of t h r e e dimensional look-up
a r rays .
-EAR A- CO.CO..IIOrn
GER- 16047
3. The i n i t i a l s t a r t cond i t ions f o r t h e forebody and a f t body a r e completely
genera 1.
4. The s t a c k i n g of des ign c a s e s is poss ib le .
5. The at tachment l o c a t i o n of t h e t e t h e r t o t h e forebody is completely genera l .
6. The t e t h e r load and t h e angle i t makes wi th t h e c e n t e r l i n e of t h e forebody a r e
program ou tpu t s .
7. -411 load and t r a j e c t o r y d a t a a r e output a t p re - se lec ted t imes.
8. Termination of a des ign case occurs et a predetermined time o r a l t i t u d e .
9. The program c a l c u l a t e s t h e e f f e c t i v e system s p r i n g constant .
10. The program c a l c u l a t e s t h e parachute phys ica l p r o p e r t i e s a s t h e perachute
i n f l a t e s a s a func t ion of time.
11. The parachute may have t h r e e s t a g e s of r e e f i n g , i f t h e automat ic d r a g a r e a
versus t ime op t ion is chosen.
12. As t h e parachute i n f l a t e s , t h e d rag a r e a ve r sus time follows a second degree 2
curve (y = ax ) .
SECTION 2 - INPm
GER- 16047
Except f o r t h e v a r i a b l e COM, a l l i n p u t s a r e r e a d i n under t h e format s t a t e m e n t
8F10.0. COH i s a n 80 column h e a d e r ca rd . A l l of the fo l lowing v a r i a b l e s a r e
d e f i n e d i n t h e nomenclature.
INPUT ITEM
a
VARIABLE
AIPHI, AIPHID, AJALPF, AJALPM, AKAM, AKAMD,
OPSYM, OPDA
PPHIE
AALPFE
AALPME
AAM
CCA
CCN
CCLM
CCY
CCLL
CCLN
CLLP, CLMQ, CLNR
PPHIDE
AALPDE
AAMD
CCLLP
CCLMQ
CCLNR
AALPPE
AAMP
CCAP
CCNP
CCMP
NUMBER OF CARDS
1 c a r d
1 c a r d
1 o r 2 c a r d s
1 o r 2 c a r d s
1 c a r d
4 t o 128 c a r d s
4 t o 128 c a r d s
4 t o 128 c a r d s
0 o r 4 t o 128 c a r d s
0 o r 4 t o 128 c a r d s
0 o r 4 t o 128 c a r d s
0 o r 1 c a r d
0 o r 1 c a r d
0 o r 1 c a r d
0 o r 1 ca rd
0 o r i t o 64 c a r d s
3 o r 4 t o 64 c a r d s
0 o r 4 t o 64 c a r d s 1
1 c a r d
1 ca rd
2 t o 8 c a r d s
2 t o 8 c a r d s
2 t o 8 c a r d s
W E A R AEROSPKE CO.CO.4llOm
GER- 1604 7
INPUT ITEM
x)
Y)
2)
aa)
bb)
VARIABLE NUMBER OF CARDS
PS, PT, EPL, EPT 1 c a r d
X, Y , X, V , GAME, CHIE, EPSI, ETA1 1 c a r d
WT, IXB, I Y B , IZB, IXYB, IXZB, IYZB 1 c a r d
S , D, XBAR, YBXR, ZBAR, OPPRIN, OPPIIOT, OPOS 1 c a r d
PSIE, THEE, PHIE, OMXBE, OMYBE, OMZBE 1 ca rd
A , B, C, OPAM, OMETRC 1 c a r d
PSIPE, THEPE, PSIPDE, THEPDE, VP, GAMPE, CHIPE 1 ca rd
LS, LTO, DLTO, NS, KT, DP, CCRIT 1 c a r d
AMAX1, AMAX2, DSX1, DSX2, AMAY1, AMAY2, DSYl, DSY2 1 ca rd
WTC, WTL, OPSP, OPDT
I F (OPSP .EQ.O .O) GO TO ITEM kk)
TRO , TR1, TR2, TR3 1 ca rd
SPRO, SPR1, SPR2, SPR3, PCTOl, PCTO2, PCT03, POROS 1 ca rd
T r I P , SSP 4 c a r d s
I F (OPDT.EQ.O.0) GO TO ITEM m)
GWAD, FSULT, AERATO, TO, DTP1, TTT, HHH 1 ca rd
DT1, TO, DTP1, TTT, HHH 1 c a r d
COM 1 ca rd
The va lues of t h e v a r i a b l e s r ead i n i n p u t i t e m "a" de termine , i n p a r t , t h e s i z e s
of t h e aerodynamic a r r a y s . The a x i a l and m m e n t c o e f f i c i e n t s have t h e o p t i o n o f
u s i n g e i t h e r e i g h t o r s i x t e e n a n g l e s - o f - a t t a c k (one o r two ca rds ) . I f , f o r example,
f i v e o r e l e v e n a n g l e s - o f - a t t a c k a r e needed, one o r two c a r d s a r e needed r e spec -
t i v e l y . The r o l l and Mach number a r r a y s may va ry from two t o e i g h t . A s a n
example cons ide r t h e a r r a y CCA where t h e va lue of CA depends on f i v e r o l l a n g l e s f
-0 - (D (di), e l e v e n a n g l e s - o f - a t t a c k (q) a l d seven Mach numbers (AM). The a r r a y
I-m
C 6 PPHIE would be r e a d i n on one ca rd c o n t a i n i n g f i v e d i s t i n c t r o l l a n g l e s , t h e l a s t m w - D .- t h r e e f i e l d s o f t e n d i g i t s would he b lanks . The a r r a y AALPFE would be r ead i n on -LL ' W i 4 nz two cards . The f i r s t c a r d ~ o u l d c o n t a i n e i g h t d i s t i n c t a n g l e s - o f - a t t a c k , and t h e
, 5 second cerd would c o n t a i n t h r e e d i s t i n c t a n g l e s - o f - a t t a c k snd f i v e b lank f i e l d s c f i: - r: t e n d i g i t s . The a r r a y AAM would be read i n on one c a r d c o n t a i n i n g seven d i s t i n c t
t Mach numbers and one b lank f i e l d of t e n d i g i t s . The f i r s t e lement i n each o f t h e
above a r r a y s shou ld s t a r t a t z e r o and i n c r e a s e numer i ca l ly u n t i l t r h i g h e s t
GER- 16047
. i r poss ib le va lue expected t o be encountered i s s p e c i f i e d . I n t h i s p a r t i c u l a r example,
t h e a r r a y s i z e used w i l l be CCA(5, 16, 7). The proper read sequence i s t o f i r s t
read two c a r d s c o n t a i n i n g t h e va lues of C a t e l e v e n a n g l ~ - o f - a t t a c k , t h e i n i t i a l A r o l l a n g l e (zero) and t h e i n i t i a l Mach number (zero). These ca rds a r e followed by
two ca rds con ta in ing values o f C a t e l e v e n ang les -o f -a t t ack , t h e second r o l l ang le A and t h e i n i t i a l Mach number. T h i s i s continued f o r f i v e r o l l a n g l e s a t t h e i . n i t i a l
Mach number. A f t e r these t e n c a r d s , t h e same procedure i s followed f o r t h e s :ond
Mach number, and t h e t h i r d , e t c . up t o seven s e t s o f t e n cards.
A l l t h e aerodynamic c o e f f i c i e n t a r r a y s a r e read s i m i l a r l y . However, n o t i c e t h a t
t h e angle-of-a t tack a r r a y a s s o c i a t e d wi th t h e moment c o e f f i c i e n t s i s d i f f e r e n t than
t h a t a s s o c i a t e d wi th the f o r c e c o e f f i c i e n t s . Also, t h e damping moment c o e f f i c i e n t
a r r a y s ( input i tems "p", "q" , and "r") may no t be read i n a t a l l , depending on the
v a l u e o f OPDA. I n s t e a d , inpu t i t e a "1" can be used i f the damping c o e f f i c i e n t s a r e
cons tan t . F i n a l l y , t h e damping c o e f f i c i e n t s correspovd t o the a r r a y s read i n inpu t items llm1l, " 11 n , and "0".
I Figure 4 he lps io c l a r i f y t h e meaning of t h e inpu t parameters a s soc ia ted wi th t h e
added a i r mass on t h e parachute (Ref. Inpu t I tem f f ) .
Po
3
n Po 1 r(
cn V
(I) D
S $4 -4 u. TI a w P 4:
Parachute Reference Diameter ( f t ) ( log)
Figure 4. Input Parameters f o r Parachute Added A i r Mass
Figure 5 he lps t o c l a r i f y the meaning o f the paramters associated with OPSO = 1.0
which d i r e c t s the program t o ca lcu la te the parachute area time h i s tory as time
ad lances (Reference Input Items hh and j j) .
GER- 16047
SECTION 3 - OUTPUT
A l l o u t p u t v a r i a b l e s a r e d e f i n e d i n t h e nomenclature. Before beg inn ing t h e simu-
l a t i o n , t h e fo l lowing v a r i a b l e s , s p e c i f y i n g t h e c h a r a c t e r i s t i c s o f t h e r i g i d body
and i n i t i a l pa rame te r s , a r e p r i n t e d ou t .
L ine 1. COM
Line 2. IXB, IXYB, XBAR, S, CLIP, OPPRIN, OPSYM, AIPHI, AIPHID, DT1 EPSI
Line 3. IYB, JXZB, YBAR, D, CLMQ, OPPLOT, OPDA, AJALPF, AJALPM, AJALPD,
TTT, ETA1
Line 4. IZB, I Y Z B , ZBAR, WT, CLNR, OPOS, OMETRC, AKAM, AKAMD, HHH
I f CLLP, CLMQ, and CLNR a r e c o n s t a n t s f o r t h e s i m u l a t i o n , t h e i r va lues a r e p r i n t e d
o u t i n t h e a p p r o p r i a t e p l ace . I f t h e damping c o e f f i c i e n t s a r e found from i n t e r p o -
l a t i o n o f t h r e e d imensional a r r a y s , CLLP, CLMQ, and CLNR a r e set e q u a l t o z e r o f o r
t h i s p r i n t o u t only . S e v e r a l v a r i a b l e s d e a l i n g w i t h t h e d e c e l e r a t o r a r e t h e n
p r i n t e d o u t .
Line 5. A , LTO, LS, AMAX1, AMAX2, AMAY1, AMAY2, AP, GWAD, FREQP, OPAM,
PCTO 1
Line 6. B, NT, NS, DSK1, DSK2, DSY1, DSY2, CHIPE, FSULT, POROS, OPDT,
PCT02
Line 7. C, DLTO, DP, LTC, WTL, WTP, CCRIT, VP, AERATO, TO OPSO, PCT03
Line 8. TRO, TRl, TR2, TR3, SPRO, SPR1, SPR2, SPR3
The parachute suspens ion l i n e load and s t r a i n a r r a y s a r e p r i n t e d o u t nex t on L ines 9
anr! 10.
L ine 9. PS (1)
Line 10. EPL(1)
The t e t h e r l i n e load and s t r a i n a r r a y s a r e p r i n t e d o u t next on L ines 11 and 12.
Line 11. PT (1)
Line 12. EPT (1)
GOOOVEAR AEROSPACE C O I . O I A l l O m
GE'-.- 16047
The parachute i n f l a t i o n time h i s t o r y a r r a y and drag a rea a r r a y a r e x i n t e d ou t
next . I f t h e op t ion (OPSP = 1.0) t h e a r r a y s a r e s c t equa l t o ze ro LCI:ause they
a r e not known before i n i t i a l time TO.
Line 13. TTIP ( I )
Line l!. SSP ( I )
The computer program then checks t h e op t ion v a r i a b l e OPPRIN. I f OPPRIN = l.,
a l l the aerodynamic d a t a i s l i s t e d a s follows:
PPHIE ( I )
AALPFZ (J)
AALPME (J )
AAM (K)
I f OPSYM = O . , t h e fo l lowing iierodynamic d a t a i s l i s t e d
CCY(I,J,K)
CCLL (1, J , K)
CCLN(I,J,K)
I n t h e above aerodynamic c o e f f i c i e n t a r r a y s , AALPFE(J) i s a s s o c i a t e d with CCA,
CCN, and CCY; AALPME(Jj i s a s s o c i a t e d with CCLM, CCLL, and CCLN.
I f OPDA = I., the damping aerodynamics i s l i s t e d .
PPHIDE ( I )
AALPDE (J )
AAMD (K)
CCLLP ( I , J ,K)
CCmQ(Iy JyK)
CCLNR ( I , J , K)
GOODVEAR AEROSPACE C O I C O . A I I O I
GER- 16047
The aerodynamic a r r a y s a s s o c i a t e d w i t h t h e d e c e l e r a t o r t h e n fo l l ow .
AALPPE (I)
(1)
CCAP (I)
CCNP ( I )
ccm (I)
A f t e r 'he l i s t i n g of t h e i n p u t d a t a , t h e computer program b e g i n s n u m e r i c a l l y i n t e -
g r a t i n g . A t T = T and a t p r ede t e rmined t ime i n c r e m e n t s , t h e f o l l o w i n g d a t a i s 0
p r i n t e d o u t .
L i n e 1. T , X, XD, XDD, FX, CAY V , TENS, XP, XPQ, XPED, F a , CDAP, CMP
L i n e 2. TXB, Y , YD, YDD, FY, CN, AM, LT, YP, YPD, YPDD, FYP, CNP, AMP
L i n e 3. TYB, Z , Z D , ZDD, FZ, CY, DYPR, TPD, ZP, ZYD, ZPDD, FZP, TYPB, DYPRP
L i n e 4. TZB, PSIE, PSIDE, PSIDDE, QPSI, CLN, ALPE, O W E , PSIPE, PSIPDE,
PSPDJIE , QPSIP , TZPB , ALPPE
L i n e 5. GAME . THEE, THEDE, THEDDE, QTHE, CLM, PHIIE , OMYBE, THEPE , THEPDE, THPDDE, QTHEP, TPDXB, GAMPE
L i n e 6. CHIE, PHIE, PHIDE, PHdDE, QPHI, CLL, PHIAE, OMZBE, KS, CLLP,
CLMQ, CLNR, TYPDRB, PULAN
L i n e 7. MPAL, WAS, DMD, QMAXPB, IXPB, IYPB, SPD, SPY SPRU, SPRL, TINT,
TNINY, T F I , XPBDI
Whel-I t h e s i m u l a t i o n r e a c h e s HHH o r TTT, t h e computer w i l l write o u t "RUN ENDED
BY CONSTRAINTS". It w i l l t h e n a t t e m p t t o r e ad i n more d a t a c a r d s , t o i n i t i a l i z e
f o r a n o t h e r s i m u l a t i o n , s t a r t i n g w i t h i n p u t item "y". I f t h e r e a r e no d a t a
c a r d s a v a i l a b l e , t h e program w i 11 CALL EXIT.
GOOWEAR C O ~ P O ~ ~ T I O ~ AEROSPACE
SECTION 4 - NUMERICAL SOLUTION
For the most general type r i g i d body, t he re a r e s i x second order d i f f e r e n t i a 1
equations, coupled i n the acce lera t ion terms. These s i x equations can be wr i t t en
and solved simultaneously using the PIVEP,T subroutine. PIVERT uses Gauss elimina-
t i o n with complete pivot ing t o obtain the l a r g e s t diagonal elements. After solving . . f o r the acce lera t ions (ui i n equation (153)), the r e s u l t s a r e numerically i n t e -
(3) grated using Runge-Kutta , fou r th order techniques . If the forebody has the proper t ies t h a t H = P = Z = Ixyb = Ixzb - - I Y Z ~ = O*, the
equations of motion g rea t ly s implify f o r the forebody. I n t h e case of i n t eg ra t ing
t h e Euler angles , three equations remain coupled i n t h e acce lera t ion terms, end a r e
separated using PIVERT. The th ree t r a n s l a t i o n a l acce lera t ions a r e already i n a
su i t ab l e form t o in tegra te immediately. A simpler s i t u a t i o n occurs i f the added
Tssses of the dece lera tor a r e neglected. A l l f i v e equations of motion a r e un-
coupled i n the second de r iva t ive and a r e e a s i l y in tegra ted by 4th order Runge-
Kutta .
GOODVEAR AEROSPACE C O . C O ~ . l I O Y
k t
GER- 16047
P i
SECTION 5 - PLOTTING ROUTINE
i I f OPPIllT 1. , e leven v a r i a b l e s a r e saved i n a r r a y s . A t t h e end of t h e simula-
t i o n , any o r a l l of t h e s e v a r i a b l e s ;re p l o t t e d by c a l l i n g PLTRAJ and s e t t i n g the
a p p r o p r i a t e arguments. PLTRAJ was o r i g i n a l l y w r i t t e n f o r use on a CALCOMP 563 i t plo ' ter and 750 t ape d r i v e . It has been modified f o r use a t M.S.F.C. where 8
SC 4020 p l o t t e r i s t h e p r e f e r r e d p l o t t e r . The o r i g i n a l PLTRAJ w i l l p l o t up t o
4 v a r i a b l e s ve r sus time on one graph f o r each c a l l t o PLTRAJ. The modified PLTRAJ
f o r t h e SC 4020 p l o t t e r p l o t s only one v a r i a b l e versus time per p l o t ; t n e r e f o r e
f o u r p l o t s w i l l be made i n s t e a d of one f o r each c a l l t o PLTL4J. Two hundred d s t a
p o i n t s a r e p l o t t e d on each graph p e r v a r i a b l e .
GOODVEAR AEROSPACE C O I C O . . l I O ~
SECTION 6 - ENGLISH TO METRIC C O I ~ S l O N
The computer program operates i n e i t h e r Fnglish or Metric un i t s . Fhe program input
and output i s i n Englisk un i t s unless t h e op t icn paraae te r , OME;TRC, i s s e t equal t c ,
1. I f OMETRC = 1. the input and o ~ t p u t i s i n t h e Metric Systom. A conversion
t ab l e from Engiish t o Metric !.s given b e l a ! f o r comonlg used engineering parameters.
ENGLISH TO METRIC CCINVERSION
REFERENCE NASA SP 7012
*EXACT
FJRCE
LENGTH
MASS
SPEED
PRESSURE
volume
ffim
ACCELERATION
INERTIA
TORQUE
DENSITY
v i scos i t y
SPRING CONSTAFT
(LR) X 4.4482216152605* = (1 ) MEWTON
(FT) X .304802~0 * = (1) METER
(SLUG) x 14.5939029 = (1) KILOGRAM
(FT/SEC) X ,3048 = ( 1 ) METERS SEC
(r.a/$ ) x 4: .882e58 = ( 1 ) N E " V I T O N / ~ ~ E R ~
( F T ~ ) x .02831684659* = METERS^ (F?) X .09290304* = MFTERS~
(FT/SEC~) x .3048+ = METER/SEC 2
( SL~G-I?? ) X 1.35581794 5= KILOGRAM- METER^ (FT-LB) x 1.355817948 = METER - NEWTON
(SLIG/FJ!~) X 515.379 = KILC)GRAM/METJB~
(SI~UG/FP-SEC) X 47.880258 = NLurON SEC/ME'TET~~
(LB/FT ) X 14.59390293 = NEWTON/METER
* ESract Numbers - No rolmd o f f s
mjsec
a/m2 01 3 2
m
m/sec 2 2
kg-m
GOODVEAR AEFICSPACE c O n P O n . l l O l
CER- 1604 7
SECTION 7 - CONCLUSIONS AND RECOMMEN'DATIONS
The i o l l o w i n g conc lus ions and recommendations a r e made f o r u s e of t h e 6+3 P.0.F.
computer progr3m.
1. The (6+5,\ DOF l oads a s se s smen t computer program s h o u l d be used p r i m a r i l y t o
a n a l y z e t h e loads induced on a wobbling o r s p i n n i n g body when t h e body i s
s t a b i l i z e d by t h e deployment o f a drogue pa rachu te .
2. A f t e r t h e body h a s been s t a b i l i z e d by t h e drogue , f u r t h e r pa rachu te deployments
(main c h u t e s ) shou ld bc ana lyzed u s i n g tt p l a n e r (3+3) DOF computer program.
The (3+3; program shou ld be used because o f t h e f o l l o w i n g r ea sons :
a . The (3+3) program i s f a s t e r and e a s i e r t o u s e t han t h e (6+5) program. 0
b. Terminal d e s c e n t w i t h ,he forebody p i t c h a i ~ g l e e q u a l t o 290. p r e s e n t s
no mathemat ica l s o l u t i n n problent u s i n g t h e (3-t3) WF program.
3. It should be noted h e r e t h a t t h e (6+5) DOF p r o g r a s h a s a mathemat ica l s i n g u l a r i t y
p o i n t a t a forebody p i t c h a n g l e o f 'go0. To permi t passage through t h i s p o i n t
t h e s i x forebody a c c e l e r a t i o n s a r e f r o z e n a t t h e i r l a s t va lu2 when t h e p i t c h
a n g l e i s i n t h e r e g i o n o f 89.8O ( < 90.2"' T h i s ' ~ d u c e s some e r r o r i n t h e
t r a n s l a t i o n a l c o o r d i n a t e s and a t t i t u d e of t h e forebody, b u t i t has been shown
t o be s m a l l f o r normal v e l o c i t y pas se s througn t h i s p o i n t . A t ime count.:
(TNINY) f o r t h e t ime s p e n t i n t h i s r e g i o n i s a program ~ ' ~ t p u t .
CHAPTER I V - PRNRAM LISTING AND SAMPLE COMPUTER RUN
The fol lowing prc:ram l i s t i n g i s f o r t h e Univac 1108 a t M.S.F.C. and adapted from
che IBM 360 l i s t i n g used by Goodyear Aerospace Corporation.
The sample problem, SRB S t a b i l i z a t i o n by 54' Drogue Palachute , r e p r e s e n t s t h e
depJoyment of a 54' drogue f r o r a SRB which i s wobbling and flyiklg broadside
t o t h e wino vector. The t r a j e c t o r y of t h e SBB is near ly v e r t i c a l . The drogue
is s t a r t i n g t o i n f l a t e and i s s t r e t c h e d ou t normal t o the SRB c e n t e r l i n e . The
drogue has one s t s g e o f r e e f i n g (0.82 of f u l l open a rea ) .
Some of tile xiore important i n i t i a l cond i t ions a r e given i n t h e t a b l e below.
- d - 4) Pl7 The output from the samp1.n problem s t a r t s on Page 57, and s e l e c t e d por t ions of
I A t l t i t u d e F t
Vel city i:t/Sec
A n g l e o f A t t s c k Deg.
F l i g h t Path Angle 3eg.
Body Axis Rates
P i t c h DegISec
Yaw Deg/Se c
Rol l DegISec
r-- -0 Ins
t h e t r a j e c t o r y a r e found s t a r t i n g on page 137. -
19,000.
553.
90.
- 85.
- 1.0
0 .O
0.0
w a x W O O - m u o x x a a m o c a O N - 0 C C r C t C W I X X * * * * N
w w n - a s x a 2 3 w w a o c r c C O O 0 - 0 r ~ a 3 ~ 1 0 c r ~ a a - ~ U ~ ~ * t r t C b - I l ~ u * * * N N
a a a a a a a a a o o a a a o o
. 3 3 r . . r - W W O O O O O U J W J J J J Z Z U U U U U W J Z Z Z Z ~ a a w w w w u J J X K X I C I C L r I . b . W W W W Y W J L L L L L L C a z w w w r ) w >
X T X Y ~ X
Z Z C C C C I - O W W J z z z z Z r f f W W W W W n W W X I L X S O - C C W W . . a W a A W W ~ t c c . - C C c J * c * * - * e u z ir. ~ C C ~ C ~ J *
IncnInmcnzPI W J Z U a r w w w a I w 0 ~ r r L a ( L a P U A C B - > > > > *
c r c c r -u a a L z o o x a a a a 0 -
000000Y) 2 z 2 0 0 Z Z Z t Z 0 - - 0 0 0 0 0 x ( n m - - - - - * w z z ~ m i n m n x 4 J - L Z z Z x C f X W W W W W C Y - - r x x I z - C o a - - - - - I
0 3 5 0 0 u e e * -
- - - I - 0 c Z Z 0 m - - z z z z z m - ----- 3
- 0 m - L L * - 0s 0 * - a 0 H n - a - - w - L U - d L m - * C I( - - f f L - - L J)Q
0 - a u - s t r m O W 0 0 a - a x = - - a - In - - A * N - A C e W m e - - - - u t - 0 uan0.u - - *a111 - 0 C - - x - 3 * - - * ul-3 - n P x o - * - m * - 3 4 - - -C -. 2 - a - s -- - - L A n a a m - ~ e - a - J - L ~ w t
m a x - * a U O P a C C - a - - * - a - 8 . -r - - a * - - 3 = - 0 * A < - - x - 0 - - - 0 0 - a m o L x - ~ e s a u x - -- 0 0 x a - - a - a * - + - X.l) L' N L A f e Q C - -U - - f - t 1 - a a - "( - ~ u r n m - r - u a a - -L -
m u rL - a - Q - - - * -- x - 0 - t * J a % - a - 4 - - - J o r , o * - . , x - w - * - W L O Y ) -L - u r n = N - - ~ m - r m x - - a -u - K O - w a - 3 a r w e - r - % C * - m c - = - I & - L d U & U a - - 3 - -0 - 3- - J a * 3 z r - _ I L U C u a W a 3 o o x o w X , O o - = a + - , .=a -. - - 3 - N O - c u o ~ a - a - ~ - m - - m a - A - - - x - - - - a & - - - a * a a o w m * - - L O O - n o a c n o a = - - * - 0 - - --- - 0 0 - a * O - a * O C ~ L ~ * u m L m o - c n a c - C - N W -- 3 J J O 3 1 - a - u - c - a x u L
m a - - w x - 0 - - N U - - 6 - r - c n - - N - -*t 0 - X C A - C z L - W - - Q x - O L - P C - - - 0- c - 2 0 3 ~ - P L Q L X U - - m a - 3 O I 3 N O - 0 u a - * * w a - I O U - N O - - - x -r - u o a - x a 0 L r A " J d a - m x 0..
n r u - - ~ I U - - z - a n - - + a - U C - a e m a - u r n - 1 0 - II)-A - r a x o r - u - - x r a ~ X L - P a -* I L Q _ U O I X L - - Q - 0 - - e s w - L O - - - r - H A * - a - - - a m - - 3 P W O 0 - N O = * l l - 6 - 4 Z O C - -::?
m x u x n O - - 0 3 - - x ~ m - u - - a m m o - * - h a * N U - L * x a - u a - X I - - m I n - - 0 - - - g x - O W - - w & * a x O - - L O - - - N X D - - a - a - -L - - - u r a s r -rdInr * - - a u m r w m - L - O W -L - O W m - a - O - * -- - - * I 2 - - a - - 0 - A
N A ~ - - - - r -z - m w r a - * z - - ~ o m u - a A z O - C I - t a n u o a - 0 I ) e - A I - a x e 0 - x - o a t : * - - - a -*to - - - - u r a a ~ r - 9 -err= - n , . - u l r - r r a a o - a u a - x as- * ( n x ( n - * - o z x - u - x - x x n - - - r x o -- - < - - a - a r * c b a - - - ~ r e -*-- I - o * - - P U I U U C X - C . - l u a t - ~ C C ~ V ~ C U Y I O ~ X
m r z w x m rc - r a w s - - - - u - - u x r -m- n - - a * - Z * Z Y ) C A - N ( ~ P - U - 4 0 - % ( + a -- ,VIP - - 0 - o n . - * - x z a m u - - -- -- - - 3 0 f - - - a m * ( L Y b Q - a m 0 ..-
L a - w - - m o m z - r r - * - 3 2 - a - - a r m s ~ a ~ a z o ~ c m o - - - W - o o - a r m c -
d - U D P - W N W X L X - O P - 9 X S a Z t - 0 % ' r r r - X E - f m a a - n f ~ ~ - ~ * n r m
W L a o a * - * - O a O m * -r r r o e u u e o e a a - - O N - O * O U N u - - U Q C U - U U 0 -
- N - N - - w n t a - - N O : ~ Q C ~
. .. a. . . b I U U uuuuuuu ------- I C C C C C C * n m n u l u l n m 0 0 0 0 0 0 0 ~~~zzzz......~...................... ~ O O u Q O u - N O l m 4 h a ~ O - h n * m * h ~ * O - ~ ~ t m * C I m * - - - - - - - - - - - - - - - - - h h r u h h h k h h - ------- P a O O O o O . . . . . . .
no
t le
i "C
l lo
P
OI 1
0
PO
! 1
1
Po
l 17
0
01
13
c
rOll
3
-01
19
P
OI
I3
C
CI
I3
P
ol 1
3
?O
t 1
3
PO
1 1
3
PC
11
3
PC
1 I
3
no
11
3
'?0
11
3
PO
I I
3
OC
I I
3
PC
11
3
f'C
II3
n
ot
13
n
o1
13
P
OI I*
Fottr
?e
ll*
P
C!
14
n
o1
I*
30
1 lY
n
ot
I*.
eo
tl
r
PO
! :*
P
@I I*
PO
I I*
PP
I I*
ro
t I
9
Po
l I*
n
o1
I*
"@
I
lY
'?0
11
* C
olt
*
00
1 1')
P
OI
14
?
01
1*
PO
I 1
q
no
1 1
4
00
1 1
s
"11
6
30
1 1
7
00
12
0
t'01)
l P
o1
13
C
+o
l3q
0
0t3
q
PO
1 3
4
00
13
7
AIP
H1
IS
TH
E
NO
* 0
1 E
LE
ME
NT
S
IN
PP
Mlt
A
RR
AY
(2
1
TO
6,)
A
tPn
lD
IS
TH
C
NO
* O
F E
LL
ME
NIS
IN
PP
HlD
E
rnW
AY
ti
?*
TO
a,)
A
JAL
CC
IS
Tn
t N
O*
OC
EL
EM
EN
TS
IN
AA
LP
Ft
rUR
Av
(
6,
O
n
16
1)
AJA
LC
M
IS
TH
C
NO
* O
f C
LL
ME
NT
S
IN A
ALP
MC
A
HR
AY
(Y
e
OR
14
,)
AU
AM
I
S T
HC
N
O.
OF
EL
EM
EN
TS
IN
T
nC
A
rb
AR
KA
V
12
1
TO
8
.)
AU
AM
D
IS
TH
E
NO
* O
C
EL
EM
tNlS
It
1 T
HE
~
Ah
0
AR
RA
Y
TO
Be
)
OP
SV
M
EQ
UA
L I. 1
1 n
oov
1s
ST
MM
ET
RIC
S
UC
H
Tn
AT
C
Yrc
LL
rcL
Nm
oe
O
CS
VW
CYO
IL O
r I
F B
OD
Y ,
IS N
OT
S
YM
ME
TR
IC
OC
OA
E
QU
AL
1.
TO
R
EA
D
AM
PIN
G
CO
CC
~I
C~
~N
TS
AS
, A
C
UN
CT
ION
0
1
RO
LL
AN
GL
E,
AN
GL
E-0
1-A
TT
AC
U,
AN
D n
rcn
NU
MB
ER
O
CD
A
EQ
UA
L
O*
To
R
cA
D
AH
PIN
G C
OE
~rI
ClC
NT
S A
S
CO
NS
TA
NT
S
&
CA
LL
P
MA
RG
t*6
6'e
*Y
I,*
2*
) I
IN
5
IOU
T.
6
IUS
CO
. 0
RC
4O
(IIN
*I)
L
~~
H~
B~
~~
~~
D~
AJ
IL
~I
,L
J~
L?
H,
AK
AM
~A
KA
P
,DP
SY
M,O
?Dh
ICI(
I~A
ICH
1
1P
w1
om
rlrn
to
JAL
PF
.AJA
LC
C
JIL
PM
-AJ
AL
CM
A
JALC
O.8
e
u u --- S z
0 0 a 0 X X L - -- z <*
- - x < - r r - Y Y ~ L f r a u r e - w - - - - a * = 2 5 5 = --- - - - - W . .
x . L Y - - m a -.. Q - - - I m x . Y X - - a x u - C x - - C - r - 0 0 0 0 0
0-0
a - - z -en x x x --I - z r a k a S X L * X Q I * n 0-
a & - * L - - - - - - . - - ..- x --- 0 -- --- e m . - - m Y ) - . a &fl --- . m-x m - - r - - - - 0 r W -. - - W --- " - - J 0 J O --- - - x u - x x w saw - m a L I a m LrnL c A 3 A - - - a % > z L A > z * < a --- -- > a - - 3 2 - - a - 2 ?2 m . .
( rx - * - a * d 7 n 0 0 --- - L a 3 - 3 - 3 - Z -L - 0 0 - - > A ---a 0- - 0 X J - . .. --- x a < - - m z - m m z a - o -- 3-,3- - w re. a a 37-• m > 7 - - 7 1 - - - - - - - 7 - 7 0 - z a - m a - - - - - - a - - - - a - - - w a y I - - y x y x - W . l * - . a x - - = w e - X Y w a m m - - - a m - 7 3 7 m 7 > ' Y Y - S e x - - B h U ->--b73 - - - - - - - -- - - > * - 3 3 C * - 0 -. - a - - . .. - - - - 7 - 7 - w 0 - y - - - --• 3 3 - - s c 0 - m o r - 7 ---- ---.- n-- --.- --- b - 3 - - - & O R ----- - w ~ - - - L ~ ~ - A z - u w u x a f z y - - - - U L t Y ~ Z d C U ~ A A C O - 0 0 - A A A a - a % L - a & - u u u - uuu - w a - a o u u u L a ( = = x > ~ s u u u > - u u u A X A L U U U a a - a - - - q 0 - - - * - A & a * - - -
> S U V U = a u u u - 3 ru & - - - - - a ( - - . .
0 . - - - - - r - a * - - - g - ----u- - - - - - a - - - L a - C W o - L C r I----- - - - - a
a o - - - - - - - -10- 0 - - - - - o m ----- ( - - a r r r - r r - - o m - ~ X ~ ~ ~ Z Z Z Z Z O L Z Z Z W O r z ~ ~ z z t z z U X Z Z L Z ( . r - - - - - r . - , & * - - - D . .me--- - - - r r .=,---m- .I*------.- (n--.-Z a - - - - - - -003t - - - - . - Q * . - - - - - - - C & - - - - C L - O - - - - - - m . . - - - - - - P m C 0 9 0 0 a 0 0 ( n 0 o 0 0 C ~ 0 0 C O O O a n o a a ~ C O O O O O ~ r x - ~ - a - ~ ~ w - - ( a z u - - a ( - * a - a x z z a . - - a ~ * a u W Y w w u Y C ~ ~ u U O C L W O W W w W W Y ~ d a O w u W W U - x a r a a l ~ e e e - = e m u - ~ ~ ~ ~ a a a a w u u u ~ ~ a a a
x a
- a *o ln a - &
* 0 W C e x a-
a n
a*) P O
u W c a 3. x U r n * * P * - a a P * rn w - C I - a - 3 n U - a m * C u - - r .m - x - U L o m t r a
0 -- 0 . . a I m m as. 0 0
0 u
P l n u 3 m 3 Y a
U * a a m
C m - c z - 3 Z 3 x m
u- - 3 a u tz u Y x
b b.- -
0 r * - 0 . .
U U e m C F U Y 1C x 0 0
- - - - - - - - . ~ - - - 0.- - c - -u P - a r L N O - u u n c n ~ n o z err z w c o t
+ - m a - x ~ t c a - W a a W - t Y * * - - x o - - - a W - b L w n - - Z Y O Id - -a& -u 0 - ens * -m - * * - ! 0 1 a - x m * - ~ % m - n a b - m - C d L - C
2 - - w - w 2 - o w o e - Y O 3 P. w 0 0
a - - a u ~ a f n r - ( D C C R a O C H I, OI t m t
w o - - f o e 0 -2 -r cn t -29fDr - - A - C * C - n ", r o - o - - o
c c , r x - u u c - . * c y e - c u e C x t w 0 - - w a 2 - C O L - a - a- - , , -. - 0 - m - - t u a - .E 2 - 0 - r U A - m a c - m - m * r c* i ~ c ~ a c a I r o e t 0
r t r - a r * * a a m a r a n c za - w CI w u U - w n u3 - w Q O D - a o o + a a a e 0 0 - - - - - - . L L - - - W . - - - o m ... & 0 - - - C C C * Y O ---------- ,7-- 0 . - a- -000 2 * * * * U W
0 - r - 9 r - - Q y o r - 0 0 - 0000 C U - r . 0 0 0000 I t .c)
3 Z Z Z Z Z Z Z l Z 2 0 O Z X U W N Z 0 0 0 0 0 0 • - U O . Z t Z O O O U Z a C b - 6 - t C O C * - - - - - - - . - - -LQ- - . - - . 3 ) - 0 0 0 0 0 • - 0 0 0 0 3 ~ c - ~ ~ e 2 - c ( L O ---------- m - . --z - r 0 . . a r , 9 9 . . . .z a- 0 . . O Z - = O m 8 . a f a.0 - - - - - - - - - - b a - - t - - - g - o o o o m I u - ~ n m - + a - 0 - c c * o - r u n O ~ O C 0 ~ o o ~ a ~ a ~ o w . o 0 ~ o a - ~ ~ 0 u ~ u 1 0 - ( ~ n 0 0 0 0 ~ c n 0 o ~ o ~ ~ ~ o a a a a C H - a u - - L Z * o - N c ) o P L ~ ~ * c o ~ ~ u - - - C .azrl en u c w c c Z U L O w ~ w Y w u ~ ~ w ~ c * , ~ ~ ~ c ~ ~ o ~ Y ( ~ . : ~ ~ & ~ & u ~ ~ o o c L ~ o w - I ~ ~ w o w u c 0 0 0 0 ( c - v
1 : e a ~ a ~ e w t m = - a ~ o c n u u : , . : c c . - t m m m n r ~ n o u - a w a u ~ ~ u t z r c t c t t I
?J W x * - m x n w 2 1 ; * C * e *+ * * C ' I
, , i V u u *
U u u - - 1 .I
C C C n m I) 0 0 0 * . . * . * . . . ..zo.......~**oo*oo*.o*ooI.o.ooo.o.oo~o.oo.*zo*
O - n ~ t a n * ~ a e o ~ - ~ ~ ~ ~ e ~ ~ ~ ~ - n c ) t m ~ ~ ~ o . ~ - ~ ~ a t m a ~ e ~ o - ~ n ~ m a ~ ~ ~ o ~ - u t t t t t s 8 t t t l l i ~ n ~ r r , n u r m 1 n n m a ~ ~ ~ ~ ~ ~ ~ e z r w r ~ ~ r n w n r w c e c c r c t r a r ~ a c - r ~ --.----------.-----.-----.--------------.----------.-.---
P 0 a
* c r r , w O L L - ~ n m m n n n m o - n t ~ - ~ u ~ ) t n ~ c o - h . n ~ ~ ~ , ~ ~ ~ - ~ - ~ n ~ r n o - m m w m * * * f i w - s m * ~ - ~ m w m e h b o o - n m n n ~ ~ n t t r t~ t w r ~ ~ ~ ~ ~ ~ n ~ r ~ ~ ~ w ~ c ~ o o - - - - - - - - - - - r u l l Y 8 l l i * ~ @ S W ~ 4 4 0 C * Q 4 Q 4 C 4 ~ < 4 ~ * 4 * * 4 * C U 9 Z 9 Q Z * h h * w ~ * W w * * * ~ * ~ ~ w O O O U ~ V C , C O O O C ' O L , O C I C , L , O O O Q O O Q ~ ~ O O O O O C ~ ~ O O C t , ~ O O O O G O O O O O O C I O O O C , O C C . C . C L C C C C C C C C L C . C C C C . C ' C C . C C t . C C C 9 0 . L C . C C . G C t c t C . C C C O C C C C . O C C . C L t o
193.
4
0
OT
Ta
~OY~SS*~~R~(~LS.~~PI~~GLOAO*WT~I~ULT~ACRATO~~
IV
Y*
rR
I?O
P-
Ir/t
IZ*
*D
TT
l 1-
5.
77
- 0
11
l9
6*
I D
O
19
70
b
I
NN
- Y
Y
I re
* I~
~N
N.G
C.O
G
O to
6
2
199.
1-
1.1
20
0.
VY.
YV
.10.
2
01
. GO
7
0
61
I0
20
b
Z
DT
I-
NN
/lO
r**
I 2
O3
* 3
CO
NT
INU
E
20
q.
GO
-32
. 1
7
2J5
. R
H0
0~
.00
23
78
2
06
. ~
Em
iOV
1b
Y3
50
~
DI~
GN
oS
TIC
O
tnt
TCS
Y
roR
C
QU
ALI
TY
B
ET
ILC
N'N
ON
-INT
~GC
RS
M
AY
NO
T
BL
MEA
N~NG
IUL,
20
70
If
fOn
ET
RC
rEO
*le
) 6
0 T
O
1.0
20
8.
GO
TO
14
1
20
9.
I40
IC
m6
37
13
77
. 2
10
. ~
~o
o=
r,
zr
ss
7
21
1.
60
09
e8
05
YIb
0
21
2.
14
1 O
PR
- 5
7.2
95
77
VS
2
13.
n
=a
l/G
o
21
4.
mP
.*TP
/GO
P
IS*
W
lCM
- IT
C/3
2.1
7
21
60
, IT
LM
- w
TL
/SZ
II~
2
17
. L
SC
LO
' SQRT(L5*LS'r10176.0r.DP)
2l8
*
AP
m
(~
~C
*(
.IS
~~
*D
~*
LS
CL
O)
*
~S
*L
SC
LO
*W
TL
J/~
U~
C+
UIL
) 2
19.
AX
- ~
~.
LO
G~
~~
A~
~X
I~
-A
L~
G~
O~
A~
AX
~~
)/
~A
LO
G~
O~
DS
X~
~-
AL
OG
~O
~D
SX
~~
~
27
0.
0x
0
AH
AX
I/(O
SX
lo*A
XJ
2
2 1
A
rm
~
AL
OG
~~
~A
~~
Y~
J~
AL
OG
IO
~A
~A
Y~
J~
/~
AL
OG
~O
~D
~Y
~)
~A
LO
~I
O~
DS
~~
~~
2
72
. B
Y-
AM
AY
l/fO
SV
Ie*A
Y)
27
-0
PS
I=P
SlE
/DP
R
27
Y
TH
Cm
TH
EE
/DP
R
12
5.
Pn
l-P
HlE
/OP
u
22
60
P
SlP
-PS
IPE
/DP
R
22
7.
1H
f.P
-TH
EP
LID
PW
2
28
. ' -
OM
XB
=O
MX
BE
/DP
R
27
9.
OM
YB
OO
HY
BE
/OP
R
23
0.
OI'Z
B=
OM
ZB
L /D
PR
2
31
. T
~E
PO
=T
HE
PO
L/O
?R
2
32.
PS
IPO
=P
SIP
OC
/D?
R
2 3
3 G
AY
-GA
MC
/OP
R
13
4.
CW
1-C
UlE
/DP
R
2 35
. G
AM
P-L
AM
PE
/DP
R
23
60
C
n)P
-Cu
lPE
/OP
R
23
7.
SG
A~
.S
~N
~~
~M
)
23
0.
CG
AM
~c
OS
tGA
Ml
23
9.
SC
HI=
SIN
(CH
II
2*0
* C
C~
I-C
Os
lCw
l~
2.1
* S
GA
MP
-SIN
IGA
MP
) 2
42
. C
;AM
P-C
OS
~G
AM
PI
2.3.
SC
HIP
~S
INIC
~IP
I ~
YY
*
CC
WIP
-cO
S~
CH
JP
) 2.5
. X
O-V
*CG
AM
.CC
~I
246.
.
YD
mV
*cG
AM
~S
CH
l -
24
70
Z
O-V
*SG
lM
n-n
- w e 0 . 0 --- --- a a a uuu UUU 0 . . -a-
a a a a ( ( + * * 0 0 0 C C C . a 2 2 0 0 0 Y
x 0 0 0 C C C C 'J 2 3 2 2 --- # @ # - - - x - - -an a
x - I- - - - U U a s a x n n m x • un m a - - - c m .u .m U U U U W - * -. .. uuu\ L V)- m - . . . - - 0 a m a m uuu- -+ m Q I n & + + + X X - + U D U - - - * a x - + - # ' N O U I I I L - - - - x - x - - 0 - 0 * H - - - - a x a x P N N N ~ ~ .
- - - p p - - - a p H a u a a W - - - Y N O o n o u w x u u u r x *
a a - - - - W - - - - w w W * W . s c u u u o o x -- - w - m r - w - m r w s I * - - x u - c u . e . + * o X I m x x a t m x a - o c x c e x T + x x u . m m m - - - u u L C D - - a C a - - C u H - - & r n C a . . L + + * x 1 . u m - - - A ~ - - - A A U . r w . a o m u ~ - a - - - ~ a w
A A A W W A : A U U . - - 0 . - . o - n w - h n ~ u x o P ~ ~ ~ ~ ~ Z Z ~ O ~ Z Z - ~ W ~ - W ~ - W ~ Q X - - 0 0 o c N ~ x r z z z v ~ z z z n m ~ a x a w x n y , x a w c - - - m m m - a ~ ~ H w ~ - - m m m - - a m C U a c U a C u D w - - - * - e 0 * w ~ ~ ) - - - z z - - - ~ ~ u ~ n ~ m u ~ m u m m m u u u ~ x m - c - U U b ~ z Z z - - n n m O o . . . . . . O . - - - u u u o o u 0.- . m ~ - - - n m ~ o o u u - - - - - - - - - - n r n ~ ~ ~ - e x a - - m a a a m m n m rn . u u m m - n o - n , m - ~ m - - - * - a D - 0 - m c m o a ~ ~ > ~ . m n a ~ r n m a ~ - - - - - - . - - - - - - + ~ + D D I F ~ ~ G D c - D D D - w - - w - w - - w - - - ~ N N ~ ~ ~ - - - x ~ N O O O ~ Z ~ ~ ~ U ~ D - O C C ~ X X ~ X ~ X X Y X - - - - - - - - - ~ P P m - l u - w - a x a z ~ m a ~ a a C a Q C a ~ a P C u u u u u u u U u U u u a a ~ V ) 1 x t z C * C 3 7 ~ n ~ ~ m w m m m u ~ u u u ~ u u u u ~ ~ ~ ~ ~ ~ ~ x ~ ~ ~ ~ a ~ ~ ~ - ~ u ~ ~
0
0 - - W 0 - C C Z u O r U - \ * -c u a -- - Q C 2 L a-0 -* z - a t a n x - 0 \ 0 - * a a 0 u \ - - u a - a - - C.0 a 3 a x 0 0 - - e
L W S N A f - - 3 • x - - a r - u I* 0 c - - a - 4 a r - w 01110- - - x o I- A & - 1 - 0 a n
* r - % m ~ o i x o a c a m - * N - - x 0 - * - - - 0 o m - ,
x - - - - x * u G S - - x x r * x - W W \ 0 - m t Q - - 0 7 r y - r - r - r a 3 x N X O - - - - m a E f 0 - C - Q L l t X * C U U - - ( r l a N * - m - a Q - - x A * K * - 0 E - a - m - - * - a - 0
r u e - n - x x - 7 - - a m - , n a n - - - m - a - 0 - w - 0 0 - 0 - 0x0 n,x 0 -
C - o e x t - - C - r - t * - 3 x - 3 - 3 x x - x - - = - O - O O N O C * - *-• 0 0 - - - - - - - " , m - m u - * - C a - C - C a - - N O -- Y a u a w e o - e - b u - C f a C x C r - - a - W C - - a - * - a - - a - x u r n ~ ~ ~ a ~ a ~ x n ~ x ~ a a B b c l L B L I C - w - I W - - m a * m - I - - - m -
\ - a 0 - - 7 t + a 0 \ - 0 0 \ - -8 In In \ 0 - 0 x - a - = a * - - L a 4 e a a a \ L L - 0 - L O U a L - a \ cn - t o - 0 - - 0 a - cn
u - 7 U L * u - - W C * L *13L x a - O m - a - - - 0 - - a m a - 0 z x C * 01.1 L 3 - 0 - h X c - (Y c a
x . - 4 < l tn u s * - - a x r e - - - "1-a - - 2 - x - * x f r n 0 - - - 0
- x 0 - C . 0 -Nu \ C m c c e 4 - - m - u a r a n - a
n u 0 7 a - - + - a * e r n - r c a a S O *L - r n - Q - a o 0 x 0 < W - - a ~ P X C - N -r - 2 x 0 K x - - 1 - 0 h U O * - - h. c Y - - a c - -a x u o r r a m c * c w e x a - * - - c t N - - - r - - 100 3 - t 4 r 0 XI. - 1 - e o a - a - C L * - 9 . a - m a + * - - -un= O L \ lllw Y O - t - - - - a & - - h u a ~ a m - - a
• O u - = - O r 0 . t - x -a. - a x - r w - m - - 9 atnatnn w x r U U L t X * - u o u a - X O C - - z-* W - a - L Y - 0 * a h & ? - c - Z W O - - a x x b - n - - - 0 - 0 - n o - - a -.( w . - t n t n c - L O O = l u x - 1 ~ 1 . - c t r d o - - V O W e t n c m w z m * o t e o 0 . 1 - b - - - - - * - ;a,=-:= 0 - z 2 0
- O m x r x - --I-- 0 0 m n -L x - o m - * = - m u a x m a - \ -a 0 . * - h - - - - -- - z c a 9 a C W rL 3 - N - x x - - t n W - ~ t c * 0 - 0 - - m r c m o - r e - - n c o - m r A - L - 0 0 - - 2 - c N c Z a o - & t n x 0 L W
P I I L L - D - - 2 - 0 a r - o a o u > m n - z t Z - O L -1a - m n r c - a - - t n \ Z ~ - 2 a - - - r c V I L L O X - - - 0 t n t . e n c m x u - - n - t r a - Y ) - c - n o w r - n w u - 0 r r X - r 0 0 O X 2 O d n - - C O O X n C a
C - 0 0 O X - x x b - O K - a > - a d - m = a m k - - - m X m h O - O U - e L O X b . O - - t a - 0h.L- - - r+ - - C - - * R U - Q W U O d - L U L x - 0 - * - - 9 - c- .d - - N \ - c - - - a - x - - N - W O R Q - - x u a n m a x w m ~ ~ - m - u 0 0 - ~ a ~ w c n ~ ~ o c - a a ~ z - c c x N - - 0 - x Z - > . a - a a + - - a - -c L O - . - - , & - , , x - r o o - n - 3 u - - - - - A&-, n- Q O - L N O - c -u - - = m e - & a k a -- w - - \ \ u - ' J - - - - - x x - - X C - - a - c n - x m a - w x \ \ *,.-,#- - \ - 0 c - v w - - - o w -- - , - - w --c\\ .. - L N n C \ x - - . ' r ( r - 0 Y C W O O - c . I I ) I a C u Z - r - - - m -- ~ 1 e 3 a n a r = \ a c m > \ -
a m r - o o \ - c - - n a - o L " c * o \ - w \ o \ m m ... N - m 0 0 - 0 L r., . - i r - .-..-. - e m - - . - C \ 3 t - C - & - N - O m - c - r - r - 0 - C - Q - - Q - I . w ~ - - . - - - a w < o o c w - P - O W - - - - I .L - t X * -. x X X P * O Q c I I - f X - b - S - - ( P O - x - a x r 0 0 c ~ m O V I - ~ O u a - Q m r - a r x c o - x e a o c a - - - - - - a w a ~ o a o ~ - m a o m ~ W L L N - B L - B L I L \ S L N - h n - * N O Y O ~ Q - N * -Hn - N
n m h - - - - - N
W x n n t o r - - 4 h 0 C - .- - *
U U U * U
35
9.
DO
*I K
mI
,*A
M
36
0.
WR
lTC
l1O
UT
~2
0)
36
1.
20
F
OU
MA
T(/
//,Y
~X
,VC
N IN
cn
MA
L
cO
EF
r A~RAY(IP~~,JALPF,KAM~)*/I
36
2'
WR
ITE
(IO
Ul,
lT)
AA
MIK
) 3
b3
* 0
0 *
I I
.l,t
PM
l 3
6**
$
1 WRITfblOUTelO)(CCNII~J~K~~J*I,JLLPT)
36
5.
00
W
Y r
ml,
wrn
3
66
0
WR
ITE
(fO
UT
,2b
J
36
7.
26
rO
RM
A7
1//
/,Y
Olt
*CL
M
(PIT
CH
M
OM
, C
OE
FW
~~
RA
Y~
IP
U~
~J
AL
PH
,K
AI
OI
*/
~
36
0*
W
Rll
E(t
OU
1,I
I)
AA
H(K
) 3
69
. D
O
YY
l.),
IPH
I 3
7 O
* O
* W
RI
TL
(~
OU
T~
~~
J(
CC
LM
(I
,J
~U
~,
J~
I,
JA
~.
PU
)
*O
~A
GN
OS
TIC
~ TH
E
TE
ST
F
OR
C
QU
AL
ITV
BE
TIC
EN
N
ON
-IN
TE
GE
RS
H
AY
NO
T
BL
ML
AN
I~G
Fu
L.
37
1.
lFfO
PS
YM
.fO
*lw
) G
O
70
9
5
37
2.
DO
4
2
K-1
,K
AM
3
73
* .
WR
llC
(lO
UT
,22
I 3
74
. 2
2 F
oR
MA
T(/
//,Y
oX
,*C
Y
(L
AT
C~
AL
CO
EF
* A
RR
Av
(~
P~
I,
JA
L?
C~
KA
~~
~~
/~
3
75
. W
R~
TL
(IO
UT
,I~
) &
AU
(u!
37
40
D
O
42
I.
I,IP
MI
37
7.
$2
W
RI
TL
I~
OU
T~
~O
)(
CC
~(
~~
J~
U~
~J
.I
~J
AL
PI
J
37
8.
OJ
4
3
K.),
KA
M
37
9.
WA
ITC
IIO
Ul~
tY)
3
n0
* 2
4
TO
RM
AT
I///
,~O
X,,
CL
L
(RO
LL
'1
0Mr
CO
EF
* ARH*YI~?HI,JALPM,UA~~)'/J
38
1.
WR
IT
E~
IO
UT
~I
~~
A
AM
tK)
31
2.
DO
4
3
I-I,
IPH
l 3
e3
* $
3 Wnllt(~0UT,~g)lcCLL~I,JDK~,J~I,JALPM)
3 9
Y 0
0
45
K.1
,U
AM
3
85
. ~
RI
TL
~I
OU
~,
~~
)
38
60
2-
F
OR
MA
TI//
/,YO
I,*C
LY
(YA
W M
OM
. co
tr.
AR
RA
Y~
IPM
~,J
ILP
H~
KA
H))
*II
38
7.
*R?
TE
llO
UT
,I7
) A
4H
(K)
30
8.
DO
Y
5
Ir)
,lP
Hl
38
9.
05
W
RI
TC
II
~U
T,
~O
)(
CC
LN
~I
~J
,K
),
JI
L~
~)
3
90
. 9
5
CO
NT
INU
E
~O
~~
GN
OS
TI
C~
T
HE
T
ES
T
FO
R
EQ
UA
LIT
Y
BE
TlE
EN
N
ON
-IN
TE
GE
RS
M
AY
N
OT
B
E
ML
AN
INC
FU
L.
3Q
1*
IF
~O
PD
AO
EO
O~
~)
G
O
TO
3
9
39
2.
GO
TO
3
8
.
.
. .
39
3.
- --
- -
JV
CO
YT
INU
L
39
**
WR
ITL
(IO
UT
,B~
3
95
. 0
F
OR
MA
T(I
HI,
IOX
,~P
PH
~D
E
(RO
LL
A
NG
LE
-DE
GR
EE
S)*
/)
39
6.
WR
IT
E~
IO
VT
~I
O)
~P
PH
lD
~~
l)
~~
~I
~~
PU
~D
)
39
7*
WR
ITE
(IO
UT
,9)
39
0.
t
~O
I?M
AT
I~,I
QX
,~A
AL
PD
E
(AYG
LI
OF
A
TT
AC
K-D
EG
RE
ES
J*~
) 3
99
. -
WR
~T
~(
IO
UT
,~
O~
(A
AI
*D
E(
J)
~J
~~
~J
AL
PO
)
*DO
* *R
ITE
IIO
UT
,6S
J
YO
I*
65
F
oR
MA
~I/
,~o
X~
*A
AM
O (M
AC
H
NU
MB
ER
)*/)
q
02
. W
R1
TE
tlO
UT
,IO
)(A
AH
O~
U~
~I(
*I~
KA
MO
) '1
03.
00
9
7
K~
I,K
AM
D
qO
O*
*R
lTt(
10
UT
c3
0)
YO
S*
- 3
0 F
OR
MA
T(/
//,~
QX
,~C
LL
F
(RO
LL
DA
RP
INC
M
OM
, c
e~
r. AR
RA
Y(I
Pw
~D
,JA
LP
D,U
A
YO
&*
. IM
D)
)@
/I
90
70
IR
llE
(lO
UT
,I7
)
AA
HD
(~
)
$08'
D
O
97
I.I,
IPHI
D Y
OV
* $
7 WRITLI~OUT,~~~(CCLL?~I~JDK)D~~I~JALPD)
'I 10
. 3
4 ~
OR
MA
Tll
XD
8F
IE*
3/)
*I I*
Do
9
8 K
=I,
KA
ND
.
$12.
IR
ITL
lIO
UT
,33
J
x - x - w t a a ~ 0 0 - z 0' a - a a - a w a ~ a a - ~ m o o = n t t o z a n o a * u a a - o w L a s * - 0 C C r. 0 C u u u a G w * L U m m o - C r r a U - a u ~ a a a b ~ 2 x 8 - a a - 0 - - 3- 0 U W L W . 0 0 P 0 u 1 - - -a-=c - . - - - n - z n W O t - U U O C 8 9 - 0 C X t X N O C X t 9 - C a r O U O 0 0 0 n - % ~ ) * n n - a . n a \ * - 3 u s z - - " w n u t z a = r r r r y r r r r r r \ r r t a r X C r . - C Q 0 a Q O \ + - - - b - \ C - C C \ C x a C L C - * 2 * C 0 0 \ 3 - 3 3 3 - 3 3 \ 3 0 x a r u 3 x c - u - - W 3 Y a * * \ - o m o r n o - 0 8 3 0 \ 0 c v u - a ~ ~ 0 - - O Q U + 2 0 F e - w r - - - - - - - - ~ X O - ~ - C C W C ~ ~ - 9 a a m x - 0 - - - C - - - C - C T < Z O z - 2 w . n C 7 0 7 - a m a * ~ - W P W D W ~ - W ~ W W a W ~ r r - - - W - - 7 I D 3 ~ 7 W - 2 8 8 t - C * C t e x . . C t C C r C x r - - - C C u C 7 u - m u 7 D - m u A W W a @ - - - m r - - - ~ - e ~ x x x z - - - - a 9 n - m - 2 2 - w ~ ~ a ~ a ~ a ~ - a ~ a a ~ a ~ ~ ~ ~ ~ n o a r r ~ r - u ~ ~ ~ ( ~ a ~ a m x L ~ ~ O B O B L - B ~ ~ I L B L N N N N U ~ - - - O - ~ O - I ~ u U P - - - * - N O n o N o m * L . a o - n n t n a a t o - 0
0
u u u
ZP
DD
*FF
tS)
CU
IF
~~
TA
N~
~S
NG
L~
YD
II
N
GL
~X
OJ
)*O
PR
G
AM
~~
TA
N~
(S
NG
L~
~O
~,
SO
RT
~~
NG
LI
XD
~*
Z~
YO
*~
Z~
)~
IF
IXO
*L
T~
O~
I bA
Mm
3*1
;15
02
65
-GA
M
GA
MP
~A
TA
N~
~S
NG
L(
Z?
D~
~~
QR
~IS
NG
L(
XP
D~
~Z
*Y
PD
**
Z~
~)
IF
lXP
De
LT
*@
ml
GA
MP
o~
*~
~~
SV
~~
~-
GA
HP
G
4WE
mG
AW
eDP
R
AL
PP
E-A
LP
P~
OP
R
GA
MP
E=
GA
MP
~O
PR
0M
X8E
.0M
X8.
DP
R
On
y~
E-o
w~
e.0
PR
O
MZ
bE
mO
MZ
R*D
PR
C
O4P
.CA
P.S
P
VD
~S
QU
T(
XD
D*
*~
*~
DD
.*
~+
ZD
O*
~~
)
VP
D=
SQ
RT
IXP
DD
**2
*YP
DD
**2
+Z
PD
D**
2l
TP
OR
B-S
QR
TtT
P~
YB
*~
2+
TP
DZ
Bo
.2)
IFlL
Ye
LT
*!T
(rl
GO
T
O
3
Pu
LA
N-
AT
AN
~IT
PD
RB
,TP
DA
~)*
DP
R
GO
TO
5'
4 5
3 P
UL
AN
m
18
0,
5'4
CO
NT
INU
E
~R
IT
EI
~O
UT
~~
~)
T~
~~
XD
,X
DO
~F
X~
CA
~V
~T
EN
S~
X~
~X
PD
~X
PD
O~
~X
P~
CD
A?
~C
M?
~
ITX
R~
Yr
YD
,Yb
D,F
Y,C
Ne
AM
,LT
~Y
P~
YP
O,Y
P~
Dr
F~
P~
C~
P,A
~P
~
2T
YR
,Z
,Z
D~
Z~
D,
FZ
,C
Y,
DY
PR
~T
P0
~Z
P~
ZP
Je
ZP
DD
m~
Z~
~T
YP
B~
0~
~R
P~
~
T~
R,
PS
~E
,P
S!
~~
~P
SI
~D
E,
QP
~~
,C
LN
,A
LP
E,
O~
X~
E,
PS
~P
~,
P~
IP
~E
,P
S~
D~
E,
Y
QP
SIP
eT
ZP
teA
LP
PE
, ~
~A
ME
~~
ME
E~
TH
E~
~~
T~
E~
~~
~O
~H
E~
~L
M~
~~
~~
~~
O~
~~
E~
~~
EP
~~
TM
E~
OE
~T
H~
DO
E~
~
OT
HE
P~
TP
PX
R,G
AM
P~
, ~
CH
IE
~P
H~
E,
?H
~~
E,
P~
ID
D€
~O
~M
~~
~L
L~
~M
JA
E~
O~
~~
E~
K~
~C
LL
P~
C~
MQ
~C
L~
R
~~
TP
DR
B~
PU
LA
N~
MP
AL
~M
PA
S~
O~
D,
QM
AX
PB
~J
P.
PB
,I
~P
~~
SP
~
9,
SP
~S
PR
U~
SP
RL
DT
~N
T~
TN
IN
Y~
TF
~~
XP
BO
I
JJ
m2
SP
EC
IFY
V
AR
IAB
LES
T
OR
L
OT
T
AP
E
TH
E
TE
ST
F
OR
E
OU
AL
llY
B
ET
WE
EN
N
ON
rlN
TL
LE
RS
M
AY
N
OT
B
E
MC
AN
ING
FU
LI
IF
~
OP
PL
OT
~E
Q*
O*
) CO
1
0 1
10
-
Nm
n*
1
X4
fN)-
T
TIt
Nl.
PS
:E
Y2
INl-
TH
EE
---
I3IN
lmP
HIE
~
YI
NI
~A
LP
E
YS
INI-
Z
Yb
(N1
mA
M
V7
tYl-
G4
ME
Y
8IN
lmD
YP
R
Y9
fN)m
TP
D
Y(g
INl=
4L
PP
E
11
0
IF
lC
ON
ST
Rl1
1~
ell
~,1
11
Ill U
RIT
EII
OU
T,S
2)
52
~
OR
MA
TI~
~X
,*R
UN
EN
DE
D
By
C
ON
ST
RA
INT
S*/
) I
F
IOP
PL
OT
) 9
9,
9V
~l
iZ
I
12
If
tNe
LC
.1)
60
T
O
09
IR
E7
l 0
G
O
TO
1
16
a a a a a a a a a a a a a
U L a a a a e ~ e a a ~ a a ~
i 0 C Y N 0 L.
- . o o o d o r - - - - - - - 0 0 0 0 0 0 0 0 0 0 0 0 0 O O O O O O O C O O O O O
W Y & L 0 - x 0
- - a 2 W W * -e L - w 4 0 a m - - - ru- m 2 2 Z a U - A & a - * c - *
.L - - a & - - L a a -1-
a n a - a u - t t L cn o w a ~ c r -- a o a - < - m - - 4 a a m * N A b - ( W m a - 0 - - U r n - L c c L - U W D ~ Y - 9 - * a n - a o 3 3 e o - - = N A P - - e c n - 3 - 0 - n m - a . x o - a - 0 9 - 0 3 - - z -K - m - - - - a ~ a a a t r m 0 3 - e 4 - - a c u d - n - r Z Z m a f - a U O m O * c - - - - * - a - - a W d - - I - u - a * - - 3 2 - q + A < - - = a * x x L - a - - 0 0 - n m o a x e - r x m u a a - W & n - -- 0 s - a - x - a ~ m - a - - a - k w m t~ -a s m ~ k 3 1 - a ~ - - u - -s - c s c n a a ( a - a a - - I L U - ~ U I - = - u r n a - - a - - ~ r * C m l -L O - - M I - a - m -- - -* --el -0 m a u a r d t x r - . ~ , r - - - ~ n o o ~ - a r o t
I - W - - u x O V ) - C - w O x N - - s c g - v - * - Ur) I l m x - O - U - x o - u n - A 3 a w * - u - < - a - x w e x c -t - ~ a e a * u a . u a - m a m m a - - 3 - - 3 l * 3 l * - i - U c U % Y S J X W * c c P C X O W o - z a - - & = a - - - - A - - - t - o
-4 -2 -4 M c u O - I a - a - a - s - a m a - - -c = a x - - - - p a - - - L * O * O W O - C o o f ? Z C ; k o a c n a a r - - - - 0 - - - - - = d o
I
- 0 0 - a - ~ ~ o ~ * ~ . c n s u m a m m - m r c c a r I X Z w a . a u a a n - 0 - u - C - a u U & - m c 0 0 A x w r - 0 - - N W - -m - a r m - - o r 0 I - - N - - - P C x t a - + a & - w - - m x - a n - a -c mtn * c - - - w - - I Q A ~ - a r m a x u - - c n c c a - Z Z - 3 0 - 0 - 0 u a - * - a # & - A 3 0 -no- -*-. Y W - x a o a - . ~ Q O L X A * A & - O X - = ooa x s a u - - a - - 2 - c n - - + a - u r - a a m c c ~ o -- N Y ) - X U t P V ) l l ) d - X P B t O * - U - O x t O A O W & O Q f x - a O u a a u D r X L 0 - O - O - - ( = w e - - m !
- S O - I - - - - H A ( r n - - - a m - - ( ~ z n - * - A * W O - N O x ~ L - 9 - Q Z Q w - - - x - - 0 - m s ~ x a n - a d - - t h a - u - - a v l x o - m a r
Z Z r - - + N - - N U - 0 - x a - u e - z a e - - Z O C ! - - x a - - m - c w - - o u a x o - - a - - - --
- a ~ - Q * I U O * N - a - -I.- - - U X ~ X X V ) -r - - - - trim -r i t n a w - a u r x w m - a m m w - x - o m ~ n z - - . - ~ a - - - a ~ R Q Q - u a * - -+ -- - m x ~ -en- - 0 -a - 0 - - - - - w - - - - y
u ~ t n - - x -z - m w r a ( * z - o ~ o m w c c ~ - w - m x - w - m x Y W - a ~ z a n a o u o a - o a * - ~ x - e x - o - u - m x x o c m x s a c a a a -r - o a - - a - * * o - - - - u a - a ~ a w a c a - - ~ v - a - - w w a a m o - c - x a w x ~ u r u o u - o u a - x 0%-u -z O - - - A A - - - A A l a a z + m z n - - x - u - x - x x o - - - a x 0 - c u m - N ~ ~ O A A ~ ~ S ( ~ A ~ . J O W ( Y ~ ~ L ~ N .
cn-* - - a a r a - - - U f 9- -*--a - o * - - 0 C l O O O ( 1 U Z Z L I O O Z Z . . a m * O c.: - ~ ~ U U O - L I D - ~ O U O C I F C U V . C O ~ L C - l e a t z t n t n z x z m m ~ . w e m u = ~ . w r n ~ l i ~ w x ~ - x a ~ m - - - - u - - U Y X - c A u ~ N o o o - I I ) ~ ~ - - ~ ~ ~ - - ~ - s ~ - s t x z r - - a + - t n ~ ~ - u m m - u - e m - x a * a ~ a t - - U - t n - - - z z - - - c n n x x c s x c N X - . - A U I ~ - -a - e - ~ t a c n u - - -- - - -I - a v m t x a z ~ z - - ~ ) c n ~ ) ~ ~ ~ a c n ~ a u - - c - - A P Q m * u o a a - a m o -- - - m o a c a e - - - m m o o o u u m m m u u m
ass - w - - 2 - z * - c - a t - * - - a r m x o - - 8 8 m w m m c n m n u u u m l l m w m mcu u w o t ~ ~ r n a o m c - - - w - 0 0 - o r m m - u c z r m w r v 0 m l m a a l l m a a m n a u ~ a ~ f a a - u ~ a - f a ~ - m ~ - r n ~ f a t x - - ~ C U D ~ O O O ~ - W - - W - W - - W W - W W - W . ~ . ~ r n ~ a c n > - c t - x m a a - a ~ * - r ~ m n r m ~ a a ~ , - W - - y l x ~ m x t n x x m x r r x x ~ x Z Z ~ W L - I O ~ ~ O ~ O X ~ - X - ~ O ~ U U ~ ~ ~ a - - m n t n x r c n a c o a c ~ ~ r a r ~ a c ~ a ~
m a - - O N - U N U - - U P C V - U W - - w a + a a c n c n m m m u u u u u m w m u u u Y W - u - N - W m P a h - r v n t m a r . a * ~ m a x C C
I
I
e9
e
so*
5 I
st*
5
3*
I*.
55.
Sb
* 57
. s
8*
s
9
bO
* &
I*
b2
*
43
* b
**
as*
bb
* 67
. b
e*
b9
* 70
.
CC
II,I
~.C
PS
I*C
TH
E
cc
tl,z
I.-s
Ps
1*c
Tn
f c
Cll
,3~
=s
Tw
E
CCt2,Il=-ClS2s3*S?sl*C?nI
ccl2.zl=sls2S3~C~sI.CCnl
CC
l2,3
l~C
T~
E*
SP
nl
CCl3,I~.-c1s2c3-sPsl.s?nl
CC
l3,z
~=
slS
2C
3*C
?s
I.
PH
I C
C0
93
1=
CT
YE
*C
p~
1
CC
pfl
,l~
=C
~S
1P
.CT
~c
?
CC
Pli
tI
=-S
PS
lP.C
Tn
E?
C
cP
l1,3
).s
ln~
P
CC
Pl2
~IJ
=S
CS
lP
CC
Pl2
,2J
=c
?S
1P
C
CP
I~,I
I--C
PS
JP
~S
T~
LC
c
C?
l3,2
~=
SC
slP
*s
T~
fC
CC
PlJ
eJ
)=C
Tn
EP
C
WG
O*
IR
E/
(Z
~R
L)
)~
*~
T
HC
T
ES
T
tOR
E
QU
AL
ITT
B
LT
lLE
N N
ON
IIN
TC
CC
RS
M
AV
"Q
T
BE
ML
ANINC
IUL*
IffO
~L
tRC
.CQ
.I,~
6
0
0
10
CA
LL
D
EN
SIZ
e?
R,R
~O
eV
S)
60
T
O )I
10
CA
LL
D
CN
SM
IZ,?
R,R
NO
, $
1
11
CO
NT
1N
UC
V
=S
O~
Tf
SN
GL
tI0
*~
2*
VO
**
2*
10
~~
2J
~
V?=SO~TlSNCLl~CD*.2~V?O**2~ZPO~e2jJ
rnw
v/v
s
AM
t.V
P/V
S
OT
CR
=~
S*R
HO
eV
*@2
O
V?
R?
*.~
.RM
O~
VC
**~
IOO.X~*CC~I,I)*VD*CC~I~~~*ZD~CCI~~~~
V@O~lO.CCtZ,l?*V0*CCt2~2I*ZO*CCI2e3~
ZO
~~
IO
~C
C~
~~
IJ
*T
O*
CC
~~
~~
I*
~D
*C
C~
~~
~~
. ' A
3 . L r. 2 L
2 < W t
u) b *
C 0 2 - I * - * - x
I : - . cn - a - * w - -u (* - C O Y 0 8 C . C " A D o r 0 -cn- C " t a* - e - 0 - ' 0 2 0 *-a : 0 0 U . *-I.. t . w
D -C - - w \ - - ) a - L Z P \ t = --- z w t - C O - - \ - Y o o n a & - 0 0 C I B- 8
0 - m - d 0 . * O . . t cn- - t w c c o C I I 0 0 3 0 3 0
0 --- C C . 0 . - a \-- - * w e - - 0
- - a 0 .+ 0 . u - -L(n w Ou 0 - - o m . = 72: - z a m t w - **.mu) a - "-- 0 . 2 0 .a . & - * - ? a C O . L O - - a--00 c Z W O Z - t C - - a m . 0 0 - W ' C Q t O * D m a C r . O U C * P . C O \
o r - m x a m o u o l * t 0 0 - - \ t r r W - D - - 3 h Z F O t L O W - Q - C
C P D O D Z -- u C V L 0 - 0 - c n P L O P - O C - 0 ~ C U * c O c O
t + C D H . O H C O t r Z r C C ~ P t - * - C h - - > . o m - P Z 0 - - w - - - - . L . O m a a a L O O O O V L O ~ C ~ L L O ~ O - - - e c m m m - ~ x u ~ -ID- - - - a - Y
W 0 - 0 - .) o x - Cr
,O 5! - 2 2 o * 0 0
t
U L.
t n 0 ....*. e e ~ e . o ~ e . ~ ~ Z ~ ~ ~ . ~ .
~ ~ o - ~ n t m ~ n a ~ o - ~ m t u m ~ ~ ~ m ~ o - - ~ m ~ ~ m ~ h ~ f i w n m n n n ~ n n n m n t -.-~o--"oo--e-----.--o.-
0 . \ .
- h r n a e ~ - ~ t m e ~ ~ - n * a e c ~ h w ~ c t ~ ~ ~ t ~ a ~ n n m n n ~ a a a ~ a ~ c h c c ~ f i ~ n n h f i f i ~ ~ h h f i n h h n f i m n n n h f i 0 0 ~ 0 C 1 0 U 0 0 0 0 0 0 0 0 0 0 0 0 0 0 C 1 0 Q c ~ e c r c 'c c . o c c c c . r o c c c o o c . c . r c
---- - a k a . n UUU. a U U V O a . . o m - N Q O O I O P P I Y r n - -**,a t b
t l - W 3 a 1 u 1 C C C C 0 0 0 0 C C C C .... W W W W O a a a a a
r g C C C C Y I - - - - a L L L L L n - - n c n
- ,
. . . .
. , ,
-I---
. -.
* '
, I,.
. ,
.I
.
, .,
.,
5.
ti*.
1
60
. 16
8.
16
20
I 6
3.
Ib*.
'
16%
. 16
6.
1 67
. 16
8.
1 b9
* b
70
* J
7 1
872.
1
7~
1
79
17
5.
17
6'
177.
1
76
. 1
70
. 18
0.
I at
10
2.
1@
3*
l09
. 1
85
. 1
86
. l8
7*
*D
~A
CN
OS
TIC
O
t 88
s
18
9.
I QO
L
9I*
1
92
. 49
3.
)9*.
19
5.
19
b.
497.
8
*8*
)99
. -
200.
2
0t*
2
02
. 2
03
. 20
'6.
20
5.
266.
20
7.
20
8.
209.
2
10.
21
I.
212.
I-. B b4
: g
.-
- S
PR
L
SP
RO
S
PR
U
w
SP
II
PC
TO
I*
PC
TO
I tR
U=
T
OT
Rl
~~
~~
~L
c~
~~
~R
~~
~N
D~
~~
CT
~~
IO
TR
JW
O'
I~
~
x?
~U
J.X
YIP
G
O
TO
19
0
2 T
INT
* T
OT
Ql
SP
RC
. S
PR
I S
PR
U-
SP
RI
PC
TO
I.
PC
10
2
TRU
. T
OT
R2
I
F~
~~
LL
~T
QT
R~
~~
MD
~~
DC
T~
~I
OT
R~
WO
T)
)
XP
8D
l. X
PD
O
GO
TO
1
90
I
TIN
T=
T
OT
Q2
SP
RL
=
SP
R2
sr
nu-
5~
~3
P
ClO
l.
PC
70
3
TRU
. 7
0T
R3
6
0
10
19
0
I9
0
SP
OE
L S
PR
U-S
PR
L T
FI.
~
~*
.~
,~
~~
~~
I~
YR
T~
SP
OL
L~
~~
~O
R~
~/
X~
BO
~
TtN
Fr
TlN
l*T
Fl
TF
UT
I*
tl
~r
-t
l~
T
TO
IF.
*S
*P
CT
~l*
TF
MT
l IO!.
TtN
f*
TO
IT
TO
IR.
TlN
f*
2**
10
1F
S
PO
I-
PC
1OI.
SP
0CL
SP
OIN
F-
~..
SP
OE
L/T
C~
T~
lF
(tT
RU
-lf
iIR
lr~
T*
0r
Ol 6
0
TO
LO*
IF
~T
~L
L~
T~
N~
~~
NO
*T
~~
~*
T~
N~
~
60
10
19
1
tW
t T
ES
T
TOR
E
QU
AL
ITY
B
ET
WE
EN
N
O&
-IN
TE
GE
RS
M
AY
N
OT
8
E
n€
4w
lNc
fUL
* ~
FI
PC
~O
~~
EQ
~O
D~
~A
NO
~~
~~
T*
T~
N~
*A
N~
~T
~L
~~
~U
~
60
TO
8
96
IF
~T
~L
E.T
OI.
~N
O.T
.GT
~T
INC
~
GO
TO
rvr
IF~T*LL.T~~R*ANO*T*GT~T~I)
60
TO
lV
7
~F~T.LL~TRU*AND~T~GT~TOIR)
GO
TO
19
8
I
DE
LS
P=
(S
PD
LL
*(
T-
~IN
~)
**
~)
/T~
MT
~*
*~
S
P*
SP
RL
+ O
EL
S?
- -
SPD-(24lSPDELJ.(T-lIWT)~/~TfATl**2)
XP
I)OR
. X
PB
D
GO
T
O
20
0
I
OE
LS
PO
S
PO
INI*
IT-T
IN~I
.
SP
=
SP
IU*
DL
LS
CO
S
PD
=
SP
DIW
F
- .
- .
- -
.
-.
-.
XP
RO
R.
XP
BO
G
O
TO
20
0
1V
7
OLL
SP
O.
SP
o1
- S
~D
INF
eIT
-70
1)
SP
* S
PR
U
OL
LS
PO
S
CO
* -S
PD
INI
XP
RO
R.
0.0
G
O
TO
20
0
11
8
SP
* S
PR
U
sp
a-
O*O
W
l)OR
. 0.
0
- - - W C) rr U (, r A m - 0 - m - 2 " 8 Cn ) - - 0 1 . 8 - e , * Z ry Y) - (Y * - - -- -- t -. P 0 - A - - - - a --a - 7 C L a t "o, t - - 0 a - t - e -
tn cn - a C - 7 - 7- >- - 0 - 2' % I - - 2 - - 8 - a - a r W l l cn • I N - 0 "-- -- -."
m - w o o m L O o q u e a - e - 8 - L cn
- = - B - an-• c - >-& - a - a Y 0 x * . C
Cu L & o u t - - N - - - * a - z a - S A
I e - 0 Y O Y ) o . & o r - w a n a u n n u m . O * o w
C - m . . a m c 0 A . t * * - 3 a * u a z u a x u a - 0 - 0 1* . a m 2: : utm..,,
. & A 8 - - O I I + O U + O U + L I 0 0 - x mu)*.. - a W a u a A u A A u a . a O *I' - . O ' O x o t u u O O I U - - a ~ n m * @ c n I a t n e \ - - + I ~
3 - * tn m o m m - c a * + - m u - a - - n . - a ~ - a ~ a - - u r - a w w A O I * ~ ~ A - ~ O Q X O ~ ~ O O - a A 8 - 2 - - s - - z - * a a o m ~ t x m . a o m n u x * 0 l n . t L , ~ x . U 1 * L t - u I - - I * - e - - * - - . = -
C - t n a x o . - o > r P C U - C ~ ~ J - x a ~ ~ - ~ o - ~ . - n x s o ~ a - c n w m ~ ~ ~ o 0 . 0 ~ 0 t e - . - m - a - b a s a - - - - - - - - - z r a m w - 9 - r C - d \ O O O L W r L P 8 r . b * * A + O I @ W a w - # - r @ - r @ U ' - - t + N o m - ~ c u o o o o o u s s ~ ~ m x m o + o o I a a * - - - - - - - - - - u I z - N - U O - ~ a t n x x x o a w t o o + 4 a o ~ ~ n * o r x w O X o - a - - a - - a - - + o - C Z - r - a - ~ n a ~ a a . r x mo - 0 - u ~ y l r - s c - u u , u - 4 - w - ~ - - z - - s - - ~ u n c r ~ m o n t a ~ e t n P ~ t - o o > * d m r a s . - 0 - 0 -- P * m u @ n u ) @ U I a s m m o c r - 1 - 3 0 m m m o m * z a a a m L I P a u - - u - - u - . - u - - m a & m Y t - t n ma0 m a c n o m m - s s w m m w ~ m m m a x - 0 2 2 - o n - - - - - - - - - a a w m ) C I Z - m a d o a m r ~ o r ~ n r m m s 0 0 0 0 m ~ ~ ~ c + c m * + c a m a & m a a m % & m - - ~ w - - - G o - I a u - o - - O * W - ~ e z ~ m r t n * x a .LLN-.- a h - ? a n a e a t z ~ t L I T s = O - a ~ ~ ~ ~ ~ ~ ~ x a a t x ~ ~ ~ a ~ o n a ~ U U ~ - I * - X L ~ L L W O I m~ m ~ a ~ ~ u z u u ~ u u s ~ a ~ ~ ~ ~ r n - o a a u a x x o a - u r r u r t r - - - - - - - - - - u ~ ~ - ~ o ~ u u u u u u u u u u u m ~ ~ - ~
g- *'
,
i 2 y. c
0 2 .
2 J . a (II o s a -
I L C * w
Y W f 2 - a I- -. W 2 t z - x < t 0
4 u
4 n Y - L x
L - w r (L L LI - z t - - - - r U U W
, - t r r O W L 0 - - 0 0 7 + u u u 'I - w x
t i , a 7 = * - - o t - 5 I - C * - 0 3 % f w - - - ,,, " , ; ; 0
I t n x r - -
I- w e a - - . o r - - U I o r - ( - o r *
2 (L - * Y a * C - - e x - 1 0 * - ? L
p s a W * L I I I ) O O - r a - I O K - - n u b - 9 - m o a t - - * o - * - -
I w r - r w r r a
w & n
Y . - - a -J I- O W - m u e 3 e w r u m - ? - - - Y
W W X J C - - # A B rn
C Y U L X t a J ( n 7 A l n m t X I t. - - Q u t ~ ) - r m - r n -. t 0 - a - w o A U L - - = J I - X X e U t o 0 - * x o t t a z t 6, % 2 W 8- m e * h Q X '5 a
0 I. u , P I % % J ~ 0 Y B E . a a I - c e o - 0 a c \ X Z A J H a n m - OW -. C C - - C O u , u , t z ~ x x - I x a r o - 0 - 2 - L L O U L L L w L > X 5 0 - 0 a * r W m & a 4 W c a A A - Y
o* L U u - e s t a A * * A * a m r~ w r m o .I B - u w - C 0 0 t 0 Y U I k ( n X - - W O . o A A - - - - .- 3 - W II - L - J a cn m c a u, m m e a t o . - w a n m m r (I) L X 0 X z s 2 3 z W - t t = * * * o x * * * 0 a U Y - u a u W * W - C V - - 0 t - m s t -I. - x r a x - - r - z - * - u * - e n . x o - C z a a - * a a A e .I - C a z a a x u u u * - V U U 5 : . & - a 0 - C t u U U 3 - U w 3 m e n r r w w w m c - - o * - - = w z O W * m u a . O W c - a - a a d o a & & Y ( - - A - = a u a - w a w t a a
c n c a-• w z x • 0 0 4 o a • o r a o* w w w a c w w w a x w a w O X • O - w 00 O - -u t c I- a = + c e * Y O a, - & I - 0 2 t 4 nc w a m * - w - o m x w r z z L * + z z * Z ~ L
m z n m - a * * • - a - rn n + a c P - - -r m - r r - r .Z o - W W X 0 o r - o m A - o n 1-02 a 9 0 m r e r L u c O r u U O -00 -c . t - C - I 0 0 0 - + + A * + t r P + w U l N ( 9 C t - * A U A U A . t t l O A U J O 4 O C O I 8 h W
- 2 - m - Y E a U-I c n m - 3 a m - - -a-• a L ~ I A O A L W - A O ~ A ~ Z I ~ Z Z W O X O t t ~ XW U 8 L . O S . L * O A * C . O x 7 L 7 O a C a * ~ z a a C L - * * A I a O * a a O C U I U L L C * * n . . Y U l ~ - ~ O * - - - V I ~ - 9 V * u r n u u u u u u - U u U u U U O u U u U
xzl W Y + n o - 0 - 0 o t t 1 a e x - .. 0 0 0 0 0 C * I -
n r r c h 0 , O 2 I I m
WUWWWW U U V W U .. w C
- , C
m m 0 0 ..*.. e . . . . . . . . . . . . . . * . * * . . . . . . . . * * * * . * . . . Z . . . . O . * . . . * . ~
~ e o - h n t m ~ h ~ e ~ - m n ~ a ~ h a e o - ( ~ n t ~ o h m e o - a n f m ~ h ~ e ~ O - m ~ ~ m ~ ~ m o ~ - ~ m m ~ ~ m a m m m m m m m ~ m m m ~ m w ~ m m ( ~ m m m ~ ~ w m m n n n ~ ) n m m r ) r ) m - n m a ~ ~ n m ~ m ~ n m m
0 0
t ~ + ~ e - - - - - - - n w ~ ~ o - h r m a ( ~ o k o ~ ~ ~ , e o - m m m r n c * n * m * n o - - n w m * ~ 0 - m n r r * - r - ~ h ~ m h ~ ~ ~ ~ a ~ h h m n n n n n n ~ r t w t t n @ n n m m m r n d w m a ~ e ~ ~ ~ ~ a h ~ ~ h h ~ h ~ ~ ~ m m m m w c n n m m r c n m m m ~ m ~ m ~ w c w m ~ m ~ ~ m n ~ m ~ r ~ r n r ~ v m ~ m ~ n w ~ @ m ~ m m n m r e 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 o O O O O O O o O O O U O O O Q 0 0 U 0 0 0 0 C , 0 0 0 0 C , O o U o O O O o 0 0 0
C . C C O C * C C C , C ~ t O O C ~ O G C C C . t t C CaC'C Ct C i t C C-C C C C I C . C C C C * C * O C 0 0 C O L . C O m C O C , O Q
-77-
e - i 0 - -
N X - - a x \ a 0 0 . U 0 2 0 m a x U U A Q r 0 U L 0- 0 u .--a a t - x 4 a x A U a . t U. m a + = * a 2 A T U 3 U A l U B u a ---Y)
0 0 . 0 a a n a . * . a II)m m a . a * * a a w o a a a B - * > D O 0 I ... m m 6 6 O u - r l u C C C l L
0 0
Y I - -
O r r r W - a 0 W N X A X Z - a a A J - - A n u U u -- I U u U 0 - I
w - - - x x W a L a 0 2
- a a I -- W o x w W W zwl w o a c C c -0 n-zkzh.zLtnu W 8 - U - U r U D . C
0 a I . 8-- * C m a a a u a a - - a - d a ~ x A 2 x x 3
O L * A * A * A & L U W ~ U U U U U U ~ U A
.I 0 9 w 0 * n
UUU
9 - h w m * h O k n t m * ~ - h n t m a ~ O - ~ N H N ~ I ~ ~ C O - ~ ~ ~ ~ ~ C O - N ~ ~ W w m # t m 9 e C e - c ~ e 0 ~ ~ 0 0 o - - - - - - ~ m - ~ m k ~ ~ ~ n n n w n n n n ~ 1 n 0 t t 8 ~ w w t ~ m ~ m m m m m m r n ~ m a m m e 9 ~ Q ~ ~ Q ~ ~ ~ ~ Q ~ Q ~ ~ Q ~ C ~ Q Q Q - ~ Q ~ ~ C Q ~ ~ ~ ~ ~ ~ Q ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ Q Q ~ ~ ~ * ~ ~ ~ ~ ~ O O O O O O O O O O O U o O O o O O O O 0 0 0 0 0 0 v O O O O U O O t r O O O O O O e e O o o o o e o o v m o o ~ ~ o ~ o O O O C C , C * O C C O t G O e c . O C I C C C . O c 0 C . C ~ C C . C O C i L C O C C < C c C * C . C O O C O O ~ C . G t ~ C . C ' C ~ O L . 0
378.
3
79
. 38
0.
38
1.
39
2*
38
3.
38
""
3n
s.
J@
b*
387.
3
n8
. J
aV
* 3
90
. 3
9 1
3
92
. 3
93
. 3*
(1*
39
5.
396.
3
97
. 39
8.
31
9*
'100
. 4
0]*
q
02
. q
03
* YO
Y.
'40
5.
'143
6.
'107
. '101.
00
9.
'110
. 11
1 1.
'1
12.
1113
0 '1
1'1.
(115
. *
I&*
1
17
. '1
18.
'1 1
9.
9?
0*
(171
. 9
72
* 0
23
. 'I 2
'1
'12s.
02
6.
(177
. 9
20
* 0
?T
o
(130
. s3
1.
(132.
- lr
fLT
.L?
.LT
O)
60
TO
0
C C S
PR
IN6
C
ON
SlA
tiT
A
S
A
FU
WC
TIO
W
OF
ELO
WG
IT]UN
C
00
60
1
mle
11
. .
- - -
.. -
PS
lt(1
Im
PS
(!)*
NS
*CS
I6C
P
TX
tt).
P
TfI
).N
T
OE
~s
Xf
ll=
EP
L~
~I.
LS
*CS
ICP
O
EL
TX
tl)=
L
PT
tll*
LT
O
60
CO
Nll
NU
E
00
6
5
1m
l~
B
PX
tII.
PS
X~
I)
J=
2
&I
IF
~P
%~
~)
WL
E.P
TX
(J
II
CO
TO
6
2
Ja
J+
I 6
0
TO
61
6
2
PT
XS
L-
~PX~~I~~T~IJ-~))/~PTX(JI*~TI~J-~~~
-- -
..
DL
TX
m
D
~L
Tx
~J
-I)*
~D
EL
Tx
~J
)-D
EL
TX
~J
-I~
)-~
TX
SL
D
Lx
tl)
= D
EL
SX
(1)+
O
LT
X
- 6
5 C
ON
TIN
UE
1.7
6
6
IF~
OL
T.L
L.D
LX
II~
I SO
T
O
67
I.l*l
60
TO
6
6
47
OL
xSL
. (DLT-OLX~~-IJI/~DLX(~~-OLX~I-II)~
TEN
S.
PX~~-II+(PX~~~-PX(~-I)~*OLX~L
US
=
TE
NS
/DL
T
CS
m
~.*
CC
RIT
*~
~R
T(K
S*
MP
L)
DA
MP
=C
S*L
TD
6
0
TO
31
3
0 T
EN
S-O
. 0A
MP.
iJ.
31
FX=FXB.CCtl~ll*FYB*CC(2DI~+FZB~CCt3~~j
FV
~F
~#
.C
C(
~D
~)
~F
~~
*C
C~
~*
~~
+F
~L
)*
CC
(~
D~
)
F?rFXB.CC(lr3)*TYb*CCt2e3~+FZ#~CCt3~Jl-k*G
-..
~XP=FXPB*CCP~IDI)+F~PB*CCP~~~I)+~~PB~CC~~~~~~
~~
P~
F~
PB
~C
CP
~~
D~
)+
~Y
~B
*C
C~
~~
D~
~+
~L
P~
OC
C~
~~
~~
~
~ZP.F~PR*CCP(~~~)+~~PB*CCP(~D~)-~P*~
QP
SI
=-
IT
XB
~C
~~
~,
~~
*~
YB
~C
C~
~~
J~
+T
~B
~C
C~
~.
~~
~
OT
ME
W-T
YB
. P
HI+
TZ
B*S
Pn
l O
PM
lmT
XB
Q
PS
~P
=-
TZ
PR
-C
CP
~~
D~
)
oT
ntP
=-r
yp
g
~~
OP
S~
~-
A.
CC
(~
D~
)-
~*
C~
~~
D~
)-
~*
CC
~~
~~
~
DADTHE=A*CIS~*B*CIC~S~+C*C~CIC~
DADPHI=-~.CC~~D~)+C*C~~~DI)
bA
OP
SP
=A
Pm
CC
P(I
,2)
-
DA
DT
WP
=A
P.C
CP
~~
,I)
OB
OP
S~
~~
.C
C(
ID
~)
~~
*C
~(
~,
~)
*~
*C
C~
~#
~~
oBDTWC.-A.5ls2-B.SIC2S3-C.s1c2C3
O~
~P
MI
~-
B*
CC
~J
~~
)+
C*
C~
~~
~~
)
OB
OP
SP
=-A
P*C
CP
(~~
~)
0lb
lMP
~~
?.C
C?
t~,~
)
. .
..
~ O
CO
TH
C~
-A*
C~
M~
*B
.S~
S~
+C
*S
ZC
~
'OC
OP
M~
~-B
.CC
(~e
3)*
C~
Cc
t2e
3j .
DC
OT
M?r
A?.
CC
Ct3
,3)
U N N W O O Q W & L a c C C C L - + + + - I-.. 1 u r u N z r 0 -0 -u n z - - - u u u z U U U W e e O W * * J - S O P O O C t a a a w C t t m O . + + D- - - -* - - - c o
C C C - 0 - r e A A A - u * - I \ \ \ - - - a -
a a a ~ u u a - L * * ( U U U O - ~ a a ' 6 e e . W . O ~ m u u x I 0 + e * . O O O I Y H O O O I a a O a Z P A a C C t L m W C C C I I . 0 C I I ~ ( D Q ~ C L . K * N n * W c n o O O O C C O P W - & A a a A a a C L C C r C C C C -
Y
e . rya a 3 • i m w o r r u o cn M I - C l n o m I - r F cnn i -a w .-ln o n m c -. i - 0 t o .m 0 ' e J I - 0 .u 0
0 0 x a m + x I* P * u- 0 P U N c O X L (Y 0 0 cno * P \ cnw . t -0 \ A X m w e m * + C N X N O b- - 4 r c u l e a
(1 w - * a x a o I
W L - - * C * m x C - w e 0 u m W c r u r c n o o - m x u c n & D f Q : .+ 0 .
I w I % ( . N-. - - v + o m O X I. 8
N O ~ C C * - a z ' L. -=-I+- cno x ruw\-• a * 0 ,
0 a V) 7 U F Z Z Z 2 2 : - z ,
a a - t n ~ t - -a .I 0 m U l P r u I f a cn = C N e m n u a m a t . L. m x s . r y ru * W *
0 6 N - O m - U - - 3 un --.ryx e. - a
- 0 N I e P L A Q O Ul W - X X - I * W ( * U M a (L - L a x - L O - * X C U C. W H ( 1 A X C ~ D I ~ * I - XSF A e H W Q u t r n N r- - * b. \ \ . W O O l m - - x - t m - + c - - - * * a Q * O m - I 9 Ul C C C .
C - 0 - W ~ O ~ I N W M ~ - 0 l - 2-1, - I - e n - + x c a - - x - * n ~ w m - \ \ .+a - . - r a * e ~ ~ - - - a t a a a nnn- ru 1 0 0 w a ~ e e w o w o a * - - I d = - - ~ c * s ~ * w - w C - u m a - W L * U - I - X + ~ C ~ C X X - z ( m u X C - U W A n - - * v a ~ t a 0 - - - 0 0 . tU9 .e -u - - . m e * - U U O D - - 0 0 0 I I ) . W ( D r , * - * x a o - - e U l - ohm a a a - m r n - I E - ( D @ P - x l W O e t w - x C t t m u ~ - a * - a . ~ t - ~ e n m x t m C W C ~ + + + O - U I O - 0 0 0 - 0 m C O X + O W - 2 - U * N ~ I I I I I I ~ I D I ~ ~ ~ U ~ - ~ U I I > S O \ LLI-w ~ - - - - - - n - ~ a w + ~ - w a - - - w - - - - - w z e r N c ) - r u ~ - N O I W I I * O . . X A ~ - C .#I I r r - - r r - r r r - Z - - W . - N a W W e - - n o - - - - - ~ r u u n a ~ - t x ~ x ~ ~ o m i l - c ~ ~ w m a ~ ~ - - - - u - - - u - - u A - C ~ - u + - I - - e o - - - 2 L ~ ~ Q O O O ~ ~ ~ W I ~ W N ~ W ~ C ~ L ~ O C ~ W W O O - o ~ O O O O O O O W - ~ ~ W U ~ U W C A U - ~ L L ~ W W W V U -- " N - w
* cn N 11) 1 - l - m
r Y 0 u (Y -. cn + -
0 - - m U - - w u u - * cn r n - 0 - c ) - * I o n - u n
i w * ' U 3 0 -3 X O U . ry x - - * a u a D U X < c. u - w a r n U l - * O N f m r m r y : - P N A N + f * - m * - - a m + W I N - u N
t r u e * - c y cn I - m - N - - . e9 O X *cl-u m cn C a t -uu N N '
Ln -NU. * U 0 . # - . a -. -
I - 0 - - u a a + u
* - x u - - u-. - - . x ~ e c a * - x a - - A N r n % # T O * H e u r n * + - a m m a " l e e u - - m e * W W 9 h m - ( 1 - e m # a X N .+ - - Q N N I-**.--- N=(n A U r @ n - u - - -
I . + + w - u w + mu - w - t - U * ru-• - u x * - e ~ : e r a w - a q u a * X O m - U ( D W L U O C N . . * m x r w e u U l W r n W D U I - e m # - x * x x - - - m * - c n c * c - 0 . * o * r * m m - u I s f 00- . -x t l-.....I.II 1 - e n - - - - - - - - . - w r u n t t n a - m a t
I - r x r r r r r r r - 0 - - + - - - - e u ( 1 m ( 1 - - w - - - - - - - - - L O H O O O O O O O O O - 0 l 0 0 O O O O O O O
4 1 u w u U V U
Y ' I) I t m I
0 . * e * . e . z e . . * e * e e e . . . e . e . . e e e e e . . * * e . . . e e e . e e e . * e e . . . e . * . t ~ ~ ~ a ~ ~ ~ - ~ ~ w m ~ ~ a ~ o - ~ ~ ~ m ~ ~ ~ ~ ~ ) - r u a ~ m a ~ ~ ~ o . ~ - ~ ~ c ) ~ m ~ n ~ ~ o - r u ~ t m ~ n m a 0 a n n n w ~ . t * w . t * t t n n m m w 1 1 1 d y l w w * S O 9 9 9 9 9 Q 9 . ~ ~ C L * L L * * C a a ~ e a a a a ~ * w * * * r t - C t * * t t t t * * * Y t t t t t Z t t t * Y t * * t t t t t t * t ~ t t * * ~ * * t 7 * * w w
0
~~
~~
)-
(-
P~
~P
D*
~~
EP
~.
S~
TH
EP
~~
XI
IY
PB
-T
PO
~~
A~
A~
*O
AO
PS
P*
~~
~R
~D
~O
PS
P~
/
~L
T~
~P
SI
P~
/~
~X
?~
~S
TH
E?
~*
~~
PB
~C
TM
E~
~)
I
~(
~I
~I
PS
~P
O~
.~
S*
S~
T
EP
~~
xM
Y?
B-
~P
O~
~A
BA
R~
DA
OT
~P
*~
BA
R~
~~
O~
~?
*
ICBAR.OCOTMCI/L~~~~HEP)/~~PB
F~
(~
J~
~-
TP
D.
A~
AR
/L
T*
~X
~)
/)
~P
~
C~
~I-
I-T
PD
.BB
AR
/LT
*FY
P)/
MP
F
FIS
)-I-
TP
D.C
IAR
/LT
*FZ
P~
~II
~
RE
TU
RN
r)
2
CO
NT
lWU
E
I n *) x
0 u - - C * .)- 3 a m 0 u r n Y 9 9 - * b - 3 t
n o Y 9- C
C) - 3 2 & * & n I d - 0 0 2 Y I C O - m m w
L 1L o--- 0 0 0 - - 0 . 9 n-- C L . Y L I -ccn e z z r --- u - 0 - a & & > I - \ Y Y Y - n m - r u u P I P * m m m t
U Y - - - m e a--r n t m a a r - e - - r - * o * L C O C l b L I Y Z w - - I - L I b e Y -
m L
a - u c a * n & m o w - a 0
W - Y E r r r r
4
r n- Inu - w v m n m u ( S
n
W C x C
n L - 0
W W A A?- *-. - C C
w u s x C C
r
In W 0 n- L 8 -. .t u- W LL n- u -
- - 0 - r -- 0 - n u * - L- 0
r t ..- - t
" A - o r 4 r -. - 0 . r n r r r r *- = ; - w ,
- I 0 - a I
- m a - 0
r r r r r~ r r r \ r I - - r r r r -- x - N - U . - N N Z : * - - " - * - - t * * - a - - o r 0 .
- 0 \ - - - \ n - n r r r
A * Z Z Y I u ' O O O I - - .a0
0 , A
HII) A ( 2 1
k , . . W w A * r - - a x t t * t r - w a a 0 ° = 0 n-r)
0 ! I
C o n 6 - P a C O O W
s a t a n z X I Y -a 0 1 W - O t
t b 0 0 0 0 00 0 0 0
- I.. k
I , -- \ I
E
I k 8.
k - - 0
Z 0 -
r 0 - 1 f.
I t b r r r r r r r r r r r r r \
- 0 r - r r - - - . - r r 'i: - 0
9 U V. w u. a i : . u a t
U W L o r - - r C I I - r -
I I
0 - 0 - - - - - o o o - r r I
- - c * r - - r - r r - - r . - - - a uun \ • u ! 2 V O W w m C Z I 0 W 0 A b - 0 - - n X U * a 0 w a *
U U Z t * * t 0 O w e t o L O r r r 0 - - r r r r r r r o r O I X w ..
r c r r C r r r r r r r r r r r ~ - o r x - 0 - r - - r r r - - r - - w - r n w a o c n o n ? f 5 5 + - x n e - 3 - a ~ m m - * * t - - a - c * u + w a • o P o o
C b - t . r m a = W e b - A A a 2 A z w - c C L t L I W Z Z L L L
r r L * Y O C * C O 3 Z H 2 W 0 0
y L r r c - - r r r r - s 0 x - - - - n - M r W - r w r - - o - o - - - r - w * ---- - N U N N ~ a r t t t ~ -2.0 r c r . - C w x x u * * * * a **** . . * * * * A * * * * A - > - - * a + Z C - 0 0 0 "--- - - " -b - - - -** n -0 -L on--& - C t m S W S t - W t C 1 - r - 0 - o r - - - - - - I - - o - o - r - - - o
n - U A A X ~ ~ R W C I - * - w * = o - u m t s o - n m t x o - w a r s o - m * r ~ Z . - m C A y y - Z z t > 0 0 5 w - w - - - - -* (1 - - - - . *n - " - - * f ----* a U ~ * - O Z * - ~ ~ W - X W u O * w m n n m n n n n n m n n n n m a nnnnna n n n m r
~ * * a ~ a L l r a t C U w U O * A A A A - A A A A - - A A A A - - A A A A - - a A a a - - n O r r r r r r r r r - - 0 - - * W U i J U U Y W Y J - b I W W d U O L D W W J d - u Y U W Y
r t * - m \ m \ o - - w m m m m a c m m m - n m c m m - m m m e m - c n m m m r - - * * O O * L - * * a m - a * * - m e . e m - m *-**-a a * * * - * n \ - g - r J A - L A r A A ~ A W - a a a a w ~ A A A W - A A A A W c u a W W W - c c o n o n o n o a .. t * X Y XLA- t - + ~ r m m m ~ r u ~ r m B B D B X ~ m a ~ r u a 0 0 c m t - r w w n m w - a u r a w -(nu
. B y
U A * c t o . - - - - - 3 - - - - 3 - 3 - r - - 7 - 3 - = , = , - 3 - 5 ) n - c . - 3
C' f a z - ~ - - ~ - ~ ~ z - z - ~ ~ t ~ ~ t ~ u ~ t ~ n z - ( ~ ~ r ) t ~ t t - m c ) ~ a - r * y r - - -* - - r r r * - - - - - * - - - - - * - - - - - #I( a < - - A X X X X X X - ~ * * ~ D A - ~ * ~ * * A - F * ~ * * A - ~ * * * * A . - C b * C C - a - * C O C O - Z 0 0 0 0 ~ ~ Z 0 0 C O A - ~ O C O O A - I ~ ~ ~ O O ~ a - n e n a r r t c ~ A A ~ ~ ~ o A - L ~ A ~ L O A ~ A A ~ L O A ~ ~ A * L O ~ ~ - ~ ~ *
0 - 0 - Q - - U - - L L l b r U L L I L U - U L L L L b - U L L L L U - U L b L L U - W -a - n w t a - h m e a m ~
x 00 0 . 0 0 0 C a a OJI - m a C - - m II, m , m
u u U U - L. c.
h I I I I"'
I w 0 . x L. * 2 t 3 A 0 u
0 Y C + * - - T
C 011 L. 6 a w % -- Z c -0 .0. - w A d h a 0 0 3 - - 2 U N O A a - 2 - m o \ 0 - 0 - C W 8 n n o - r 0
w n 8 O C Y Y ~ 2 - - u 8 0 o w a w . - f u n
m y - . --I.
* a - I - m - - - a w i m n W I L L W - + + - - l l ( L - T r
m e - * - 0 0
- - - a"\-
* x m u - - - 2 0 c u o n
w O N
z a o - OIQO
z no z 0 0 0 0 0 Q
0 * * a n *
Y Y Y Y Y Y Y Y Y Y U Y Y Y U V U U WY Y U Y Y U Y Y
Y U
a r m * a 0 0 0 0 . 0 0 0 0 0 0 0 0 0 t
n
0 - r 0
t II) I * cn - a Y
nh 6 h - 0 . 0 0 * h a t 0 - 0 0 . - a .n 4 0 0 C ) 0 0 I C) - r t a h 8 - - W a * & - * O W N - a n - * * O t r a n - = ~ ) s r - \ - a n - s o o r * O D o - c n c s m * N I ~ - n n a m 0" 0 0 - * 1 4 . a . N m - h e - N ( Y ~ O W . O - m w m # - - a -Lno 4 a h r Z a r rcr.64 - 4 a m - 0 . - 0 0 a s m c n O t O h e r r 0 a n O o o N - - - h u b * h a t . * a = eooa-• - UI 0 . r - - I u n 0 ) s - W . O . e N U l I - - 4 - 0 - 0
D f ( D - - 0 - c g - a u n a - 0 a B . B - D ~ N n r u 0 - t a o r s - r r n - - N - n r m ~ - - m a o * W N - - m e n n't - - \ O = W - B - n 0 . o o ~ r u a r n - c n a 9 9 c g - m c n - w a 0 0 - a * - m a * m - - r * w - m t c g m - N t 0 * I I O t C
r r . 0 - a - - # N - C I - I I W r r r * . - a n r h * r h - d - n h t ~ u ) u b w - . r - 0 r . ~ ) c a r ) u ; ) * m e n - h z - m - a - r* . 9 n o - o o o n 4 . . C I O h O r * # 8 rn NQOdh.9 - 0 - O o 0 - * o h - * - r r - 0 r U 0 a m . o a v o - \ a . - t h 4 r * * * ~ n * w o - 0 . m ~ m c n u - v - c n - u a n t - - - J ~ o - 3 n - 4 r l * D n - . w t O
- o s m t t - ~ ~ r * o - r m u ~ s o r * + - - - = = ~ncg - ~ ~ h ~ , n a ) s - r m o ~ n 3- N -. . Q .u --. - 0 0 - 0 r n u u - L O - ~ m n a ~ t h r r * r m o h m I r - - - -C I
11) - O , ~ - - - * O D V . Q o m a w + r a t t 0 0
V . n h - N a - - n - 8 .- - 3 - 8 u m O * N * F . - a * - u * n w m - m t c w w w w - u * - ~ ~ Z ~ N O * O P W - * - Z - N O . Y ~ - -
o - - N ~ - N . o o ~ N - 0 . ~ m r c g ~ - m t x - * - - - o - N O - - I n 0 - I n - - m - * m a a w m c n s r r q m o r m a o O N P P - Q O 0 0 -
C r a m - n * h J O s * . m o a w * * * t o 3 - - - - r * 0 -00 - . O N - t r O N O o * n L W ~ 4 O * N r * O * N Q 4 -(n * h r * I n o Q O h
~ a m t h h ~ o - N - o m I I - - w t o a a - N - i ) 4 t i ) . N r r U - 4 W - 0 - 0 .t m r N - . D ( Y U - Q -L - . * . a n - - - t c n - n \ - ~ - 0 m - c n ~ n ~ a ~ r u 9 0 0 u w c n a - - I e r , v t o m ~ n t n - - r e t u - - r * m u u I I * a - o m z - w t : - - - m - m n w m c r t m r w d w m - m a - w n - . o n t t . ~ o n n - ~ n \ n ~ , o * o e w * a m O - s n u o - e - c n ( ~ 9 0 0 - h m * = - c n n - o a
u x - a - o - m n m s o - r s r , r c r u , c - o n W W * - - - - N m O Q Q m 9 9 ~ 0 1 ~ N ~ 4 a ~ a z o n c ~ a ~ r u - o ~ w ~ r . c n w o + . c * o ~ n - a o \ - m * \ m o o - x 0 0 I ~ - a o m o n t m C - m 0 0 n m N 0 0 ~ Q O ~ x h ~ Q h \ . l l 8 ~ w m r o m r ~ a o o o m c e - a - - - t a h ~ * ~ 0 N O - - * . . N N * r . * - # I - a m w * - - - m = x + C C b C
- . - eo I. - - W W
- - L - m k t N r - x n - . L - N - w - m . w t n - + * a - a c e 0 t . -a x . n n - r - n x W I I l h
W e - - - w t - n 3 B - r r r m z n x m m - m m m - 2 -(I- * - 0 L a .: t ma ma 2 . 1 1 8 - r x m o 0 t W L W O W X W W t U - - O t t - * O C L a u I w
o o o o m o o u o o o o o o o C O O O O U C C O Q O Q 0 0 0
A m - * X C 0
0
0 *
r. oln
I - 0 Y . nl)lr - 0 -
m b Q O - *
h Q t - - - .* d r. 0 0 0 n -4 o m B a = - .*- L n h B - - W h - D 0 0 0 - P O W N - * r o t - 3nn * N r . h \ - n n - * * O m -or)* 0 - m o - m a u m - nlr-0 0 - - e n - + # Q . r . 1 9 . N * - - 0 - U N r . r ) W . O - m u # .--.a -n* a e r . t * o J I . h h W * * * N - 0 . - 0 0 0 9 Q m O l a h 0 - 00 a n - * I C B Q P N - - - h m t r . * t . t m t -00n-• - m - a n - 0 01. 0r.N N - 0 .O.N,n a - s * - - 0 -
**a - - 0 - a - 0 W a c - o 0 m B - e n - U N o - t r l ) - 8 - N W h - u-n)r m m - r c u ~ o Q W U ~ - w - h . u t - - \ v f r y q r o n - o r --e41nan - m m Q * Q - m d l - N 0 0 0 - - * - * a * n - w t r r y n y l r r g m - u t m . a @ 8 o + +
- . r 4 - a - - 8 u - a - a @ w 0 - N r - ~ N v - -rhi3*1 -n h t d r , m c u - - 0
- I D L I Q ~ ~ C I J I ~ I C W * a - r . t - a - a - N * 3 W 0 - 3 0 > 0 0 0 N D h r ) N B 1 8 N Q C I U I Q - 0 - 0 0 - 0 - m o h w n - - 0
- a N 0 3 4 r 0 3 0 0 - \ m . - t r . * *** ( ~ ~ Q N O - 0 - m h ~ n o - t - m - @ . r r , t - - - n o 0 - 3 n J - e w - 0 0 - . W t O
r a o m t r - : > . + N O - - r o o m - + v r m - - - r t - n o .ar.Nm a h - t o * n o a o cn b - . a * O N - - . - 4 0 r . m . . . w B - h O u W O ~ ~ N ~ C C Q - n - h i u B - - - - - w e
- 0 m - h - Q O D h Q 1 U Q W t m Q t * l l - n r . - ~ m - - = - ~ * - - a - o 8 B @ * c u e 0 c 0 - n O - U Q m W Q - m t W W W W W h I B x * 0 t h N 0 9 n t W - Q - t - W 0 N 1 1 ) - - a - - N ~ - N . Q c , L . N - - ~ o r . b o * m r
- 9 - - - 0 - N O - - a n 0 r U i N - W - 8 ~ I - ~ m m a t r . n ~ ) o t m ~ o - ~ Q - ~ D C ) - m - a N a n - . * r . O U 0 9 - m o m W r * * * o o -- - r e e - 0 0 0 . S O U - t 0 ~ 0 0 o e - W ~ O Q O P N N O . N O * - a r r . n d l o w a r - - 0 n 0 ~ r . h n o - - - o m I B - - w t o a m u u N - O J ~ ~ - N ~ U - Q W - - - 0 r r r r - 0 - C N U - Q - C - e m - a m - - - ~ n - h x x z - 0 m - m h a ~ ~ n u ~ o m 8 t m * - 0 ) # a c n o u w ~ m - - - m t e - O N O W B 8 B * w 0 - C )
Z - N . 1 - - - N - O a N W L n L l t L n t J W W 0 - O m - w m - ~ m t f r u e n r . - a * u t m o . o o w h a 0 r u O - D O N O r ~ - m m a 0 0 - h m ~ t - m ~ h o *
v x - 0 m e - @ n u d e - u t ~ n * r r m . o - n m w w 0 - - - - N U I o ~ ; o r ~ ~ u o a o n - * n z a z o n m ~ u r u - ~ r . r u ~ l r r n ~ n * - Q - o m - a o \ - r c t - n o o - N - 0 B N - ~ o m o u + n * - m 0 0 n E h O o t C O Q ~ N h a ~ \ . ..) 3 1 4 m x o a r u a - 0 > m c ~ - ( r - - ~ s h n 0 ~ C A Z r u t - -..- N N B ... - 1 B N a m w t a e * m 5)s C C C C C a 0 - a * * * c n o o o a n 0 0
- ~ n - ~ n t - w - n o * - w n w
n w
(I) - W C o - n m b * m
r c. L;; a x -a*
0 ah.= - - • r r n O 0 - - - . e m
N - --L - O h * . m - w - r n - e t a *
w - W - m . ~ a * -a E + C * - O t t *
w - x • -a r o n r - - 0 - n r w - 0 0 I -hl 8 - r - w C - O *
103 # a * urn0 - x z m r a m - 8 D m - 2
) - m u C - 0 a . a - 2 ~ * L C & I a
- 2 I. r 8 - L . I 0 C O o k o I U L W O W X H w Z O " V - - O C C - C U C I e w a Y
(* m *
C .......... r.." ....................a. r.... .......... 0 3 0 U 3 a 0 0 0 0 0 a a 0 0 a 0 0 0 0 0 0 0 0 0 0 0 0 3 0 0 0 a 0 0 0 0 0 3 a 0 a 0 0 a 0 J 0 0 0 0 0 0 0
i 0 0 0
I - 0 0 0
1 - N O $ 0 0 . .
C O O U C C I I uu i
0 0 0 G O . - t r - K O
C C Z D O 0 O C X . . O
0 Q I (Y
' G O O 0 O O C 0
I * 0 C N
i I 1 , ! " i I-
i : ! ; f : * I , ! "o I C : II)
10000 0 0 0 0 ln a n a
i ~ 0 0 0 0
O O C . 0
N F ? O O . . * n m a z ' O Q t o A t b b O O d O
3 C' C a a a 5 -*t -.La , - C O m m . .. U * N t*
I
0 0 0 0 0 0 c :o
O * O . . I - o n - o a * 2 m a V L W *
P o n o - o u o a o n 0
a w a n . 0 .
a a w a n rt - a m x cn I a
0 O O C 0 06 0 m OON
N N C N . .. . * - - p - N o - m a & -J I 0 L D W I I ) N * V N
r tc W O W 0 . -
L & Z 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
- - w a . . . r x c z O O P O a m . * * a - 0 z 0 4 -Nr.
It' ' I S 0 - 0 f C1 4 0 0 r n-• 8
4( 0 . O
a n i 1 I
O C O a d o o o o l a a a m m a ? ? t i X * N N I
, t I
0
i
. * * 0 0 s e o o * * U t o o r - r - 0 0 - * *
an -end
I 1 , .a
- - - 5. - - . - s x x = r a 0 0 a 0 0 a 0 0 - 0 0 0 0
I Y 0 0 Y 0 0 Y 0 0 Y 0 0 Y 0 0
I - * 0 ... 4 - - C) * - L L' . 0 t 9 . I . . s ? 2 x . . L L r) n L 0 n L a a 0 7 ?, I I a . I a a I I A a 1
1 3 . a - ' - ' 'I : 3 t - 1
7 3 1 . 1 - 1 i ? o o - o o o o o - o o - o c i ' I . - .L 0 0 - 0 0 - 0 0 0 0 c. 0 Q
T h W t rY tu f em n x n a 'L * e 6 (L . 9 . . L - 0 8. L 9 . - - - - m - 0, a 0 - m a - , * * - ~ C I 1 :
. . - I I - I I - 8 I, 8 I - ! I ei - a - L * * L .
1 = a e a . I I I o o a ' 0 0 a o o .c o o I .
I r 0 0 a R Q 0
i e o o = o o e o o o o e
i - - a - - * a - a * i t 1 , ' ? : I :
I . . . . . - • ry w . w N . n w I L J a L I ' L L' I Y Y Y ,
0 : I c : ! ? g ! i 0 0 U 0 0 U 0 0 U 0 0 U 0 0 0 0
.: I 0 0 0 0 0 0 " . I 0 0 • n P Z - 0 e * - 0 m a . o
i i 0 C I t u
9 . . . r o . . I ! f m . . t 0 s 0 i l
$.
[ b r , .. i:
L " ? 0 * 0 t t l o * I * O W r . L r 0
x x ! I t x i u 0 0 u o o u r n d r ~ u 0 0 - 8 0 0 I .
C 9 0 0 C . 0 0 D- . 0 + - 0 0 c .. 0 P - 0 c A ~ n - 0 n a - 0 2 - - 0 I y h - 0 * w . 0 L . L 9 % * 9
mi Q, N W - L L'
r r ir - I j /
t U . X U J * 2 4 U I P O . ~ t O Q
0 0 0 a - '? a CL I . . - :* - .
w r e * 1 I
I ! 0 0 0 0 , . u lbl 0
j 9 - 9
I n , ,
u N I
-. I ! : 1 1
1 0 0 0 0
, ! ? 9 * 8 - - I
0 U 0 f! -, h
- ,. . , ,
c 1 1
i ' I , . ,
I
C i i . . .. . . . . . . . . . .
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 O O O O O O O ........
o o o o d o o O O O U 0 0 0 0 0 0 0 0 0 0
'C a 0 Y O
Q Q O ' L O 0 1 0 7
0 L O
O h n n - 0 n o o o o n
- 0 0 0 . 2.. 0-N o m m r a m
I
n o . . 0 . 4
Q O ~ ~ O O - N O I m • • n n 4 U,
0 - C - 0 0 O t O O O O . - 0 O O C I
m 0 0 - 0 . t O O O O O O - 0 N V O U o . . . b o o .
OD O O O C O C O O O O . . 0 0 0 . 0 . . ..
. . 0 0 0 .O 0 0 0 o O N - m m o o t .
C 1 . .ma
W W W m m m
a m - N
Q 0 0 0 O h 9 ~ 0 0 0 . 0. o o m 0- 0 0 .. - 0 a c n o
0 (Y ..... m
0 0 0 0 0 . 0 0 - m e - t 4 t m
, ? 0 O t b t 0 0 N Uh S T 4 O . a O N O O O ...... w o n I -
1 W W W 0 0 0 - w - n
....... .- .. i' )-...a- t - 0 . - . O m . . . 0 - l P - . r . W - e e ~ e r m a m o ~ o r o - <', 0 - O a a - C O t O n N - + O r O - J ) - + 0 - 0 - 0 - n o n o n n -
k C O W = m o + @ S O W .* C 1 1 - 9 0 t C ' O @ N C . C * O h I ~ m m o 0 8 a h a m Q * * & a a n n n S O o m n . t n r n 0 0 0 - 0 0 - m m n .- or * O ~ C B '
4 . . C Z * N & & ( . - r a m c h o o * z n 0 - m m m N o * 9 n o m 2 ~ g ? , I . U U C C C C * C) - - m e l. - - , r - - ( Y C ' - - - 0 . h - I.,,.
- , 0 - 0 0 - -. - - 0 - 0 - - - - 0 " . . . . . . . . o o . . - - 0 0 . . - - * O O - . * - O O . r e r e O C ) . e # - e O
a 1 n 0 . 0 ~ 0 0 n 0 - 0 - 0 0 ~ o - O r 0 0 O O n b * 0 0 rn o n o r c r o * 0 r o a o o i
- u a p r mm t o o 9 1 m n o o N @ * 0 0 0 r I + w o o o ~n P ~ L I I , f t ~ t Q O a t h 0 W e 0 C N O * c l * C - 0 0 b O z ? g : @ : z ? g
i ' - O N O Q 0 0 * o n o o o o t o o o m o o * o m o a 0 0
. O D O . . . D N O . 1 1 ; . . 0 h - . . Q 0 r e . . . ,
1 '
z a @ 9 r o y r t 8 8 I @ @
i . . . O C e . . . . . . . .o m 0 , O n N - o t o a - - o r r c o r n ~ . 0 r h o t o * * - h . . 0 0 - h . . o m (
m t n r a t n ~ I ! - c u m - h e 4 0
n - n - 1 ' 1
8 , I i . o . 0 * 0 n . t . O t O h
W W * ; C I - O C O N 0 Q - C C e l . cn D Q C 0 ~ 4 0 0 0 0 r P ~ O O C * . . . . . . . . . . . . . . . . . . . W C O L r L a * 8 h C O O * - o 9 r a c o o o n a - m w o o.-o w 0
- 4 - r -
.> m o m o o n o a0 0001.1 w u N - - a o o . h - = m o o -
mu-( e r n - n o o t . J ) - Q 3 o . r &&-...a 0 .* .. e m o . . . . . mrn
S * A X X O t o m @ o n o m @ o > * O ( L I H m - o 9 m - o 9
n * M l o o o n m * t o o o n n - 0 0 0 0 . m - * t o o .
m - m o o t m . e o o . r , O . . . . . . . . . . . , v 0 0 9 8 o * * I ) , m - o 4 m -n 9
W W Y O m t n o r o n a o r n o s o r 0 0 0 0 O O W V h O . , C I Q O O O . O D O R n o - o m o n r o w e + on
* - r r r 0 -
h on o n c - * O t a a - 0 r r e - 0 0 0 0 I.-n n - N * r C 0
0
r r r o r 0 0 - 0 0 0 n 0 0 o a t t o e - N * r m
r r r r r O O 0 0 * 0 * 0 0 m a - t o o * - 0 - 0 n t o o - C m - a - a
. r r r 0 0 O O O I G O B - - 0 0
r r r 0 0 0 C oa O L O O * 1 - - 0 0
r . . O h - 0 r
II) O t 0- t - n 0 I m o n * r 0 - C . . O * -
a t n r - c u m w -
Y Y U 'V. c - 0 x o m * - u c * x t x L C A * o o o m
= a * a 0 0 0 - - * n 0 0 r . U - C I O O * r * r r r n - * a II) -*
. - 0 o r .- - o* on0 r -ma a ; * m * * * h . ' n u
C I . ' I
N 0-0-00 ' t on oh o r , * o m o - 0 . 0 . o r . - I * - - n , * ' i , i
8
* or on o o 0 o* 00 00 a o t on o r , . - 0 . . . n 0 1 0 1 t n t
II) I . . l l o m O t
0 0.. 0. 0 . 0 * 0. Odl * .. r
e t n a - -.
U Y U I O O O I 030. .
0 o* 0 -00 t O* OC 0 0 O 0-0111 0 0 * o m o r o o r o* 0 - 0 0 m 0 9 0 - 0 0 ~ o r , o o o . * o o o m o r a o * o o o - r r r r . .. 0 . . . r r r r - t J 8 l 1 1 I C) * # a @ * + m a t a n t a n
a a f . . m
a
r . r O* or) m on o* 0 - n w e e o r ,
1 .. r m a - a -. -
I
l - o r on a 0- olm 0 . n 9 o o o m
a 0 . 0
0 r r 0 - 0 4 O O h 0 0 0 . - * a n o n
O r r - r m n
'dl I.
REFERENCES
1 ) Goldstein, H . , Classical Mechanics, Addison-Wesley Publishing Company
Inc . , Seventh Printing, 1960.
2) Wells,DaveA.,LagrangianDynamics,SchaumPublishingCompany,
N.Y., N.Y., 1967.
3 Korn, G . A., Mathematical Handbook for Sc i en t i s t s and Engineers,
McGraw-Hill Inc . , 1961.
4 Doyle, G. R . , Jr . , Three Dimensional Dynamics of Two Bodies Connected bv
an E las t i c Tether - Six Degrees of Freedom Forebody and Five Degrees o f
P-;ecZz= !?n==lnr=r21r~ CER- 15957, July 1973.