7/25/2019 Gogu Catalina
1/24
DATE
Parametrii de frunctionare a cuptorului
-capacitate, t/an-temp mp la intrare in zona de convectie-presiunea mp la intrare in zona deconvectie,bar-raportul abur/mp (la intrare in zona deradiatie)-tensiuni termice,kJ/m^2*h - in ZC:
- in ZR :
-sarcina maxima admisa pe oserpentina,kg/h:
-conversia finala (la iesire din ZR):
Caracteristici Mp (benzina)
-densitatea (relativa)
Fractia "combustibil lichid greu :
-densitatea:
-factorul de carcterizare:
Caracteristici constructive:
-diametrul exterior :
-diametrul interior :
-lungime tub (tronson):
cap 3696 ti c 70 C
pc0 7 bar
rs 0.5
tc 6000
tr 25000
Smax 250
xf 0.8
d420 0.730
d420comb 1.10232
Kwcom b 9.61550
de 88.910 3
di 73.6610 3 m
L u 8 m
bucDurata anuala de functionare cuptor
Timpfunct 840 h
Constanta gazelor ideale: R 8.3145 J
mol k
Fractie masa in gaz cracat f r.masa.g.c 0.27 C2H 4
7/25/2019 Gogu Catalina
2/24
Coeficienti pentru curva VE
A
2.9747
1.4459
0.8506
3.2680
8.2873
10.6266
7.9952
C
0.4209
0.1287
0.0817
0.6214
0.9340
1.1025
1.0737
B
0.8466
0.9511
1.0315
0.8274
0.6871
0.6529
0.6949
t
55
74
90
102
114
130
152
T t 273.15( )
T
328.15
347.15
363.15
375.15
387.15
403.15
425.15
K%v
0
10
30
50
70
90
100
Caracteristici MP(benzina)
d420 0.730
SG 0.9952 d420 0.0080 K
SG 0.734556
Curba VE la 1 atm
Tve A TB SGC
tve Tve 273.15
Tve
352.50068
362.394455
362.667278
363.853958
372.745294
380.010837
385.082907
tve
79.35068
89.244455
89.517278
90.703958
99.595294
106.860837
111.932907
C
7/25/2019 Gogu Catalina
3/24
Curba VE la p=6.75 bar
pve 6.7
Tf
1393
41 1log pve
3.2041
Tve 1393 Tve 41
0.998
1
1 1log pve
3.2041
Tve 1393 Tve 41
0.998
Tf
429.428698
440.826936
441.140489
442.503855
452.694716
460.990439
466.764847
tf Tf 273.15
tf
156.278698
167.676936
167.990489
169.353855
179.544716
187.840439
193.614847
0 50 10060
90
120
150
160
170
180
190
200
tve tf
%v
7/25/2019 Gogu Catalina
4/24
Alegerea numarului de serpentine/circuite in paralel
Capacitate materie prima , in tone de materie prima pe an:
CapMPcap
f r.masa.g.c
tM P
an
CapMP 1.344 105 tM P
an
D M PCapMP 10
3
Timpfunct
k g
h
D M P 1.6 104
k g
h
Ns24000
Smax
Ns 9.6 se rotunjeste superior si la un numar par Ns 10
Gm p D MPN s
Gm p 1600 k g
h
Densitate materie prima :
SG 0.9952d420 0.0080
SG 0.734556
Factor de caracterizare (Watson)
tmVf74 90 102 114 130
5 Temperatura medie volumica
tmVf 102 T mVf tmVf 273.1 TmVf 375.15 K
K w1.8 tmVf 273.15( )[ ]
1
3
SG Kw 11.943572
7/25/2019 Gogu Catalina
5/24
Masa moleculara medie:
Mm p 42.9654T m V f 1.26007 SG4.98308 exp 2.09710 4 T m V f 7.78712SG 2.0847610 3 T m V f SG
Mm p 102.013276 k gkm o l
2.1 Dimensionarea sectorului 1 din zona de convectie
In aceasta zona benzina (materia prima ) este in faza lichida .Dimensionarea acestei zone presupune determinarea diametrului interior al
tuburilor , diametrul exterior si lungimea.
Grosimea peretului (tubului) :
gde di ( )
2
g 7.62 10 3 m
Lungimea totala a serpentinei si numarul de coturi:
Determinarea lungimii se face cu ajutorul bilantului termic in zona deconvectie .
-temperatura de intrare a materiei prime in zona de convectie: tc0 70 C-temperatura de iesire a materiei prime din zona de convectie: tc1 156.27 C
Entalpiile la temperatura tc0 si tc1:
htc0 2.964 1.332SG( ) tc0 0.003074 0.001154SG( ) tc02 0.0538K w 0.3544( )
htc0 149.443913 k J
h g
htc1 2.964 1.332SG( ) tc1 0.003074 0.001154SG( ) tc12 0.0538K w 0.3544( )
htc1 363.569775 k J
h g
Lungimea serpentinei
L c1Gm p h tc1 h tc0 ( )
de tc L c1 20.444958 m
7/25/2019 Gogu Catalina
6/24
Numarul de coturi :
n c1 2
Lungimea echivalenta:
L ech 1 L c1 n c1 di 7
L ech1 31.493958 m
Coeficientul de frecare:
k este coeficientul de rugozitate k 5 m
0.01k
di
0.314
0.037598
Densitatea medie:
c0 d420 1000tc0 20
2.29 6.34d420 5.965d4202
c0 686.575974 k g
m3
k g
m3
c1 d420 1000
tc1 20
2.29 6.34d420 5.965d4202
c1 611.644343
mc1c0 c1
2 mc1 649.110158
k g
m3
wc0Gm p
c0 d i2
4
wc0 546.862408 m
h
Viteza medie de curgere:
wtc0wc0
3600 wtc0 0.151906
m
s
wc1Gm p
c1 d i2
4
wc1 613.857702 m
h
7/25/2019 Gogu Catalina
7/24
wtc1wc1
3600 wtc1 0.170516 m
s
wmc1wtc0 wtc1
2 wmc1 0.161211
m
s
Caderea de presiune pe sectorul 1 in ZC
p1 L ech1
di
mc1 wmc12
2 p1 135.593786 Pa
pc1 p1 10 5 pc1 1.355938 10 3 bar
Vascozitatea cinematica:
x v exp 4.717 0.0029254 SG Kw( )3 xv 806.283049
c0101.766
SG Kw( )0.52529.263
10 6
K w
x ve
x v
tc 0 273 c0 5.096527 10 7
m2
s
c1101.766
SG Kw( )0.52529.263
10 6
K w
x ve
x v
tc 1 273 c1 3.177536 10 7 m
2
s
mc1 c0 c1
2 m c1 4.137032 10 7 m
2
s
Vascozitatea dinamica :
c0 c0c0
c0 3.499153 10 4 k g
m s
7/25/2019 Gogu Catalina
8/24
c1 c1c1
c1 1.943522 10 4 k g
m s
Am calculat Reynolds pe fiecare
sector(mi-a spus Corina)mc1
c0 c1
2
kg
m s Re0
wtc0c0 di
c0
mc1 2.721337 10 4 Re0 2.195498 104
Re1
wtc1c1 di
c1
Re1 3.952815 104
Reynolds
Rem s1
Re0 Re1 2
Re
wmc1mc1 di
mc1
Rems1 3.074156 104
Re 2.832453 10
4
Presiunea de iesire din sectorul 1 in ZC
pc1 pc0 pc1
pc1 6.998644 bar
Presiunea medie pe sectorul 1
Pm 1
pc0 pc1 ( )
2
Pm 1 6.999322 bar
Temperatura medie pe sectorul 1 in ZC
T m1tc0 tc1
2
T m1 113.1395 C
2.2 Dimensionarea sectorului 2 din zona de convectie
tc2 193.61 C
Htc2 532.17 210.61SG 1.8213 0.45638SG( ) tc2 0.0023447 0.00059037SG( ) tc22 4.187tc2 837( ) 0.07 K w 0.84( )
Htc2 736.932051 K J
kg
7/25/2019 Gogu Catalina
9/24
Lungimea serpentinei 2
L c2G mp H tc2 h tc1 ( )
de tc L c2 35.649015 m
Numarul de coturi :
nc2 5 coturi
Lungimea echivalenta:
L ech 2 L c2 n c2 di 7 L ech2 63.271515 m
Caderea de presiune pe sectorul 2 :
dpc2 1 1 05
Given
dpc2 Gm p
2L ech2
36002 2 di5
c1100 Mm p pc1 10 5 dpc2
R t c2 273.15( )
1
c1
R t c2 273.15( )
100 Mm p pc1 10 5 dpc2
2
pc2 F in d dpc2 ( ) 10 5
pc2 0.499837
Presiunea de iesire din sectorul 2 in ZC
pc2 pc1 pc2 pc2 6.498807 bar
Densitatea medie:
k g
m3
c2
100 pc2 Mm p
R 273.15 tc2( ) c2 17.082656
mc2c1 c2
2 mc2 314.3635
k g
m3
Viteza medie de curgere:
wtc2Gm p
3600c2 d i2
4
7/25/2019 Gogu Catalina
10/24
wtc2 6.105325 m
h
wmc2wtc1 wtc2
2
wmc2 3.13792 m
s
Vascozitatea dinamica :
c2 0.0092696 0.001383 5.97121010 5 M m p
1
2
273.15 tc2( )
1
2 1.124910 5 M m p
c2 8.727447 10 3 tc2 c2 10 3 tc2 8.727447 10 6
mc2 c1 c2
2 m c2 4.4609 10 3 kg
m s
Reynolds
Rewmc2mc2 di
mc2
Re 1.628858 104
Reynolds :
Re2
wtc2c2 di
tc2 Re2 8.802553 10
5
Rem s2
Re1 Re2 2
Rems2 4.598917 10
5
Presiunea de iesire din sectorul 2 in ZC
pc2 pc1 pc2 pc2 6.498807 bar
Presiunea medie pe sectorul 2
Pm 2pc1 pc2 ( )
2 Pm 2 6.748726 bar
Temperatura medie pe sectorul 2 in ZC
T m2tc1 tc2
2 T m2 174.947 C
7/25/2019 Gogu Catalina
11/24
2.3 Dimensionarea sectorului 3 din zona de convectie
In aceasta zona benzina (materia prima ) este in faza vapori .
tc3 60 C
tc2 532.17 210.61SG 1.8213 0.45638SG( ) tc2 0.0023447 0.00059037SG( ) tc22 4.187tc2 837( ) 0.07 K w 0.84( )
Htc2 736.932051 k J
k g
tc3 532.17 210.61SG 1.8213 0.45638SG( ) tc3 0.0023447 0.00059037SG( ) tc32
4.187tc3 837( ) 0.07 K w 0.84( ) Htc3 1.950461 10
3 k J
k g
Lungimea serpentinei
L c3Gm p H tc3 H tc2 ( )
de tc L c3 115.868964 m
Numarul de coturi :
n c3 1 coturi
Lungimea echivalenta:
L ech 3 L c3 n c3 di 7 L ech3 198.736464
Caderea de presiune pe sectorul 3 :
dpc3 1 105
Give
dpc3 R Gm p2 L ech3
100 36002 2 di5 Mm p
pc2
273.15 tc2
pc2 10 5
dpc3
273.15 tc3
273.15 tc2
pc2
273.15 tc3
pc2 10 5
dpc3
2
7/25/2019 Gogu Catalina
12/24
pc3 F in d dpc3 ( ) 10 5 pc3 0.5577 bar
Presiunea de iesire din sectorul 3 in ZC
pc3 pc2 pc3 pc3 5.941107 bar
Presiunea medie pe sectorul 3
Pm 3pc2 pc3 ( )
2 Pm 3 6.219957 bar
Temperatura medie pe sectorul 3 in ZC
T m3tc2 tc3
2
T m3 396.8075
Densitatea medie:
k g
m3
c2
100 pc2 Mm p
R 273.15 tc2( ) c2 17.082656
c3
100 pc3 Mm p
R 273.15 tc3( ) c3 8.34831 k g
m3
k g
m3
mc3
c2 c3
2 mc3 12.715483
Viteza medie de curgere:
wc2Gm p
c2 d i2
4
wc2 2.197917 104 m
h
wtc2wc2
3600 wtc2 6.105325
m
s
wc34 Gm p
3600c3 di2 wc3 12.492967
m
h
wtc3 wc3 wtc3 12.492967 m
s
7/25/2019 Gogu Catalina
13/24
wmc3wtc2 wtc3
2 wmc3 9.299146
m
s
Vascozitatea dinamica :
c2 0.0092696 0.001383 5.97121010 5 Mm p
1
2
273.15 tc2( )
1
2 1.124910 5 M m p
c2 8.727447 10 3 tc2 c2 10 3 tc2 8.727447 10 6 kg
m s
c3 0.0092696 0.001383 5.97121010 5 Mm p
1
2
273.15 tc3( )
1
2 1.124910 5 M m p
c3 0.014923 cP tc3 c3 10 3 tc3 1.492324 10 5
kg
m s
mc3tc2 tc3
2 m c3 1.182535 10 5
k g
m s
Reynolds Re3c 6.97 105
Re3 wmc3mc3 dimc3
Re3 7.365357 10
5
Reynolds :
Re3
wtc3c3 di
tc3
Re3 5.14793 105
Rems3
Re2 Re3 2
Rems3 6.975242 105
7/25/2019 Gogu Catalina
14/24
Cap3 Compozitie gaz cracat si entalpia acestuia in functie de
conversia x si temperatura T
SGcomb 0.9952d420comb 0.0080
SGcom b 1.105097
tmM ef 10
Give
Kwcomb1.8 tmM ef 273.15( )[ ]
1
3
SGcomb
t mM ef F i n d t mM ef ( ) tmM ef 393.417455 C
T m Mef t mM ef 273.1 K TmMef 666.567455
Mcomb 42.9654TmMef 1.26007 SGcom b4.98308 exp 2.09710 4 TmMef 7.78712SGcom b 2.0847610 3 TmMef SGcom b
Mcomb 249.999546 k g
km o l
Masa molara fara abur
M g u x ( )1
0.0248769874x 0.0098026456
g
m ol
Masa molara cu aburg
m ol
Mg x( )1 rs
0.0248769874x 0.040332759
7/25/2019 Gogu Catalina
15/24
Cap4 Calculul serpentinei in zona de reactie
Calculul serpentinei in zona de reactie se face prin discretizare; (discretizare -simplificare pe intervale)
4.1 Calculul sectorului 1 din zona de radiatieAici are loc reactia de cracare termina - piroliza
tr 0 600 C tr 1 736.8 CC
T r0 tr0 273.1 T r1 273.15 tr 1
T r0 873.15 K
T r1 1.01 103
K
pr0 pc3 5.941107 bar pr1
x1 0.11436 x0 0
Principiul de calculPresupunem tr1 si x1 la iesire din sectorul 1.Prin bilant termic determinamlungimea serpentinei.Calculam caderea de presiune, timpul de stationare (timpul de reactie) si in finalcu ajutorul ecuatiei cinetice verificam conversia presupusa.
Intrari : Hsp0 si QfurR1
Entalpia specifica:
0 1031.9472 1033.6167x 0 0.42552368 0.14655111x0( ) T r0 1.9110402 0.19126293x 0( ) 10 3 T r02 0.3733196810 6 x 0 T r
Hsp0 796.558502 k J
KgMp alim
H h cr 0 G mp H sp0
k j
h
Hhcr0 1.274494 106 Entalpia gazului cracat la conversie x0
Iesiri
1031.9472 1033.6167x 1 0.42552368 0.14655111x1( ) T r1 1.9110402 0.19126293x 1( ) 10 3 T r12 0.3733196810 6 x 1 T r
Hsp1 1.460749 103
k J
KgMp alim
H h cr 1 G mp H sp1 Hhcr1 2.337199 106
I1 1.909490T r1 T r0 ( ) 0.0174227110 3 T r12 T r02 0.103761910 6 T r13 T r03
7/25/2019 Gogu Catalina
16/24
I1 303.637295 k J
k g
Q s1 Gm p r s I1
Qs1 2.672008 105
Q fu rR1 H h cr 1 H h cr 0 Qs1
QfurR1 1.329906 106
n R1 2 cot
L R1QfurR1
de tr
L R1 19.047122 m
L ech R1 L R1 n R1 75 d L echR1 30.096122 m
dpR1 1 1 05
Given
1 R L echR1 Gm p2
1 rs( )2
3600
2 2 di5 p r0 M g x 0 ( ) 100
T r0
pr0 105
dpR1 Mg x 1 ( )T r11000 T r0
pr0 100 M g x 0 ( )T r11000
pr0 105 dpR1 M g x 1 ( )
p F in d dpR1 ( )
p 8.011124 104
pR1 p
pr1 pr0 pR1 10 5 pr1 5.139995 bar
pr1mpr0 pr1
2 pr1m 5.540551 bar
Densitatea gaz cracat sector 1 din ZR
r0p r 0 M g x 0 ( ) 100
R Tr0 r0 3.144964
k g
m3
7/25/2019 Gogu Catalina
17/24
r1p r 1 M g x 1 ( ) 100
R Tr1 r1 2.197227
k g
m3
r1m r0 r12
r1m 2.671095 k g
m3
wr0 Gm p 1 rs( )
r 0 di2
4
1
3600 wr0 51.402027
m
s
wr1 Gm p 1 r s( )
r 1 di2 4
1
3600
wr1 73.573421 m
s
wr1mwr0 wr1
2 wr1m 62.487724 m
s
1L R1
wr1m 1 0.304814 TIMP DE SEDERE
T m r 1 T r0 T r1
2 Tmr1 941.575 C
k1 1011.25
exp210000
R T m r 1
k1 0.398449 s 1
xverf 1 1 exp k 1 1( ) xverf1 0.114367
4.2 Calculul sectorului 2 din zona de radiatie
tr 1 736.8 C tr 2 786.8 C
T r1 273.15 tr 1 T r2 273.15 tr
7/25/2019 Gogu Catalina
18/24
T r1 1.01 103 K T r2 1.06 103 K
x2 0.45897 pr1 5.139995 bar pr2
x1 0.11436 x2 x1 x2 x2 0.57334
I2 1.909490T r2 T r1 ( ) 0.0174227110 3 T r22 T r12 0.103761910 6 T r23 T r13
I2 113.953844 k J
k g
Intrare : Hsp1 si QfurR2
Hsp1 1.460749 103
k J
KgMp alim
Iesire :
2 1031.9472 1033.6167x 2 0.42552368 0.14655111x2( ) T r2 1.9110402 0.19126293x 2( ) 10 3 T r22 0.3733196810 6 x 2 T r
H sp2 2.11632 103
k J
KgMp alim
Qabur 2 Gm p r s I2
Qabur2 1.002794 105
kJ
h
Bilant termic
Qfu r R2 Gmp Hsp2 Qabur2 Gmp Hsp1 QfurR2 1.149192 106 k J
h
Lungimea pe sectorul 2 din ZR
L R2QfurR2
de tr L R2 16.458912 m
ncR2 2
L ech R2 L R2 n cR2 75 d L echR2 27.507912
Masa molara cu abur la conversia x2
M g x 2 ( ) 28.3904979657725899
7/25/2019 Gogu Catalina
19/24
Caderea de presiune pe sectorul 2:
dpR2 1000
Give
2 R L echR2 Gm p2 1 rs( )2
36002 2 di5
p r1 M g x 1 ( ) 100
T r1
pr1 105 dpR2 Mg x 2 ( )
T r21000
T r1
pr1 100 M g x 1 ( )
T r21000
pr1 105 dpR2 M g x 2 ( )
p2 F in d dpR2 ( ) p2 1.239569 105
pR2 p2 10 5
pr2 pr1 pR2 pr2 3.900425 bar
pr2mpr1 pr2
2 pr2m 4.52021 bar
Densitatea sector 2 din ZRk g
m3
r1
pr1 M g x 1 ( ) 100
R Tr1 r1 2.197227
k g
m3
r2pr2 M g x 2 ( ) 100
R Tr2 r2 1.256442
r2mr1 r2
2 r2m 1.726834
k g
m3
wr1 Gm p 1 r s( )
r 1 di2
4
1
3600 wr1 73.573421
m
s
wr2 Gm p 1 rs( )
r2 di2
4
1
3600 wr2 128.662926
m
s
wr2mwr1 wr2
2 wr2m 101.118174 m
s
7/25/2019 Gogu Catalina
20/24
2L R2
wr2m 2 0.162769
Tm r2T r1 T r2
2 T m r 2 1.035 103
k2 1011.25
exp210000
R T m r 2
k2 4.486778 s 1
x2 1 x1( ) 1 exp k 2 2( )( ) x2 0.45897
x 2ver if x 1 x2 x2verif 0.573338 x2 0.57334
4.3 Calculul sectorului 3 din zona de radiatie
tr 2 786.8 C tr 3 849.4 CT r3 273.15 tr 3
T r2 273.15 tr
T r2 1.06 103 K T r3 1.123 103 K
pr2 3.900425 bar pr3
x2 0.5733 x3 0.8
I3 1.909490T r3 T r2 ( ) 0.0174227110 3 T r32 T r22 0.103761910 6 T r33 T r23
I3 145.174884 k J
k g
Qabur 3 Gm p r s I3
Qabur3 1.277539 105
kJ
h
H sp2 2.11632 103
k J
KgMp alim
3 1031.9472 1033.6167x 3 0.42552368 0.14655111x3( ) T r3 1.9110402 0.19126293x 3( ) 10 3 T r32 0.3733196810 6 x 3 T r
Hsp3 2.637656 103
k J
KgMp alim
7/25/2019 Gogu Catalina
21/24
Bilant termic
Qfu r R3 Gmp Hsp3 Qabur3 Gmp Hsp2
QfurR3 9.618918 105 k J
h
Lungimea pe sectorul 3 din ZR
L R3QfurR3
de tr
L R3 13.776363 m
ncR3 2
L ech R3 L R3 n cR3 75 d
L echR3 24.825363
pr2 3.900425
Caderea de presiune pe sectorul 3 :
dpR3pr1 10
5
2
Give
3 R L echR3 Gm p2 1 rs( )2
36002 2 di5
p r2 M g x 2 ( ) 100
T r2
pr2 105 dpR3 Mg x 3 ( )T r31000
T r2
pr2 100 M g x 2 ( )T r31000
pr2 105 dpR3 M g x 3 ( )
dpR3 0 d pR 3 p r2 105
p3 M i n err dpR3 ( )
pR3 160764.09 Pa
pr3 pr2 pR3 10 5 pr3 2.292784 bar
pr3mpr2 pr3
2 pr3m 3.096605 bar
Densitatea gaz cracat sector 3 din ZR
r2p r 2 M g x 2 ( ) 100
R Tr2 r2 1.256442
k g
m3
7/25/2019 Gogu Catalina
22/24
r3
p r 3 M g x 3 ( ) 100
R Tr3 r3 0.616809
k g
m3
r3mr2 r3
2 r3m 0.936625
k g
m
3
wr2 128.662926 m
s
wr3 Gm p 1 rs( )
r3 di2
4
1
3600 wr3 262.08692
m
s
wr3m
wr2 wr3
2 wr3m 195.374923
m
s
T m r 3 T r2 T r3
2 T m r 3 1.09131 103
3L R3
wr3m 3 0.070512
k3 1011.25
exp210000
R T m r 3
k3 15.804445 s 1 constanta de viteza
x2 0.57334
x 3 1 x2( ) 1 e k3 3 x3 0.286669
x 3ver if x 2 x3 x3verif 0.860009
Verificare MACH
1.0 Coef adiabatic
V s R T r3
M g x 3 ( ) 10 3
Vs 624.742416 m
s
Mawr3
Vs Ma 0.419512
7/25/2019 Gogu Catalina
23/24
xT
0
x 1
x 2
x 3
0
0.114368
0.57334
0.86
L R
0
L R1
L R2 L R1
L R3 L R2 L R1 ( )
0
19.047122
35.506033
49.282396
0 10 20 30 40 50550
610
670
730
790
850
0
0.2
0.4
0.6
0.8
1
tR xT
LR
tR
T r0
T r1
T r2
T r3
273.15
600
736.85
786.85
849.47
7/25/2019 Gogu Catalina
24/24
Compozitie efluent (gaz cracat "uscat" la iesire din ZR)
gf i
0.012
0.072
0.004
0.275
0.04
0.139
0.036
0.087
0.066
0.049
0.14
0.08
M i
2.016
16.043
26.038
28.054
30.070
42.081
44.097
53.084
56.108
58.123
M m p
250
n figf i
M i
gfi
1
xf 0.86
gigf i
x f
n fi
0
0
1
2
3
4
5
6
7
8
9
10
11
-35.95238110
-34.48793910
-41.53621610
-39.80252410
-31.33022910
-33.30315310
-48.16382110
-31.63891210
-31.17630310
-48.43039810
-31.3723710
-43.210
gi
0
0
1
2
3
4
5
6
7
8
9
10
11
0.013953
0.083721
-34.65116310
0.319767
0.046512
0.161628
0.04186
0.101163
0.076744
0.056977
0.162791
0.093023
n spigi
M i
0.55
18.1020.030383
gi
1.162791
n spi
0
0
1
2
3
4
5
6
7
8
9
10
-36.92137310
-35.21853310
-41.78629810
0.011398
-31.54677810
-33.84087610
-49.49281510
-31.90571210
-31.36779410
-49.80278810
-31.5957810
n spi 0.036275
0.036275 3.7209310 4 0.035903