1.
CLS
OPEN "c:\teplo1.txt" FOR OUTPUT AS #1
'
Pt = 30: Pm = 30
'
tt =.017
'
mg = 37.6
mh = 37.6
'
uh = 70.11
ug = 79.21
'
a1 =.028
a2 =.192
'
t11 = 35
'
t12 = 185,5
'
th = 65
'
th2 = 210
'
hn =.64
DEF fnc2 (t, c) = (3.385148 + 7.466386E-02 * c - 2.860758E-03 * c ^ 2 + 7.927979E-04 * t - 4.596554E-04 * c * t + 1.707478E-05 * t * c ^ 2 - 8.930671E-07 * t ^ 2 + 5.633727E-07 * t ^ 2 * c - 2.216021E-08 * t ^ 2 * c ^ 2) * 1000
DEF fnla2 (t, c) = (970.5798 + 1.731223 * c +.176437 * c ^ 2 + 2.125481 * t +.215597 * c * t -.188884 * c ^ 2 * t - 3.949857E-04 * t ^ 2 - 6.791743E-04 * c * t ^ 2 + 5.350655E-05 * c ^ 2 * t ^ 2) *.0001
DEF fnmj2 (t, c) = (1837.506 - 40.59476 * c + 1.507156 * c ^ 2 + 1.548409 * t +.19233 * c * t - 4.712997E-03 * c ^ 2 * t + 2.525072E-03 * 6 - 3.551515E-04 * c * t ^ 2 + 5.197373E-06 * c ^ 2 * t ^ 2) * 1E-08
DEF fnc1 (t, c) = (3.385148 + 7.466386E-02 * c - 2.860758E-03 * c ^ 2 + 7.927979E-04 * t - 4.596554E-04 * c * t + 1.707478E-05 * t * c ^ 2 - 8.930671E-07 * t ^ 2 + 5.633727E-07 * t ^ 2 * c - 2.216021E-08 * t ^ 2 * c ^ 2) * 1000
DEF fnla1 (t, c) = (970.5798 + 1.731223 * c +.176437 * c ^ 2 + 2.125481 * t +.215597 * c * t -.188884 * c ^ 2 * t - 3.949857E-04 * t ^ 2 - 6.791743E-04 * c * t ^ 2 + 5.350655E-05 * c ^ 2 * t ^ 2) *.0001
DEF fnmj1 (t, c) = (1837.506 - 40.59476 * c + 1.507156 * c ^ 2 + 1.548409 * t +.19233 * c * t - 4.712997E-03 * c ^ 2 * t + 2.525072E-03 * 6 - 3.551515E-04 * c * t ^ 2 + 5.197373E-06 * c ^ 2 * t ^ 2) * 1E-08
DEF fnAm (P, t) =.9928 +.000756 * 10.2 * P +.00367 * t + 1.37E-07 * 10.2 * P * t
tc2 = (th2 + th) / 2
tc1 = (t12 + t11) / 2
Q = mg * fnc2(tc2, a2) * (th2 - th)
t12 = t11 + Q / (fnc1(tc1, a1) * mh)
dtm = th2 - t12
dtb = th - t11
IF dtm > dtb THEN SWAP dtm, dtb
dtc = (dtb - dtm) / LOG(dtb / dtm)
tc1 = (t11 + t12) / 2
tc2 = tc1 + dtc
ug = ug *.1 * fnAm(Pt, tc2) / Pt
uh = uh *.1 * fnAm(Pm, tc1) / Pm
ro = mh / ug
la = fnla2(tc2, a2)
c = fnc2(tc2, a2)
mj = fnmj2(tc2, a2)
'
dz =.012
'
dv =.09
'
del =.0015
'
k = 1000
'
ltr = 10
a1: F = Q / (dtc * k)
ftr = 3.141 * dz * ltr
ntr = F / ftr
ntr = INT(ntr + 1)
PRINT " Ntr="; ntr
aa = (3 + SQR(9 + 12 * (ntr - 1))) / 6
aa = INT(aa)
Dtr = tt * (2 * aa - 1)
PRINT " D="; Dtr
|
|
|
INPUT " D"; D
w2 = ug / ((3.14 * dv ^ 2 / 4) * ntr)
Re2 = w2 * dv * ro / mj
Pr2 = mj * c / la
Nu2 =.021 * Re2 ^.8 * Pr2 ^.43
alf2 = Nu2 * la / dv
ro1 = mh / uh
la1 = fnla1(tc1, a1)
c1 = fnc1(tc1, a1)
mj1 = fnmj1(tc1, a1)
l = 2 * SQR(2 * hn * (D / 2) - hn ^ 2)
n = 2 * aa - INT((hn - D / 2) / (tt * SIN(3.141 / 3)))
fmtr = (l - dz * n) * h
w1 = uh / fmtr
Re1 = w1 * dz * ro1 / mj1
Pr1 = mj1 * c1 / la1
IF Re1 < 1000 THEN Nu1 =.336 * Re1 ^.5 * Pr1 ^.36 ELSE Nu1 =.132 * Re1 ^.65 * Pr1 ^.36
alf1 = Nu1 * la1 / dz
k1 = 1 / (1 / alf1 + del / 46.5 + 1 / alf2)
IF ABS(k1 - k) <.01 THEN GOTO a2
k = k1
GOTO a1
a2: F = Q / (dtc * k1)
PRINT #1, " :"
PRINT #1, " -"; dtb; ", "
PRINT #1, " -"; dtm; ", "
PRINT #1, " -"; dtc; ", "
PRINT #1, " tc2="; tc2; ", "
PRINT #1, " tc1="; tc1; ","
PRINT #1, " "
PRINT #1, " ", "", ""
PRINT #1, " ", c1, c, "/*"
PRINT #1, "", la1, la, "/*"
PRINT #1, " . ", mj1 * 1000000!, mj * 1000000!, "**10^6"
PRINT #1, " , Q="; Q; ""
PRINT #1, " "
PRINT #1, " alf2"
PRINT #1, " n="; ntr
PRINT #1, " ug="; ug; "^3/c"
PRINT #1, " ro2="; ro; "/^3"
PRINT #1, " ftr="; ftr; "^2"
PRINT #1, " w2="; w2; "/c"
PRINT #1, " Re2="; Re2
PRINT #1, " Pr2="; Pr2
PRINT #1, " Nu2="; Nu2
PRINT #1, " alf2="; alf2; "/2*"
PRINT #1, " alf1"
PRINT #1, " n="; ntr
PRINT #1, " uh="; uh; "^3/c"
PRINT #1, " ro1="; ro1; "/^3"
PRINT #1, " fmtr="; fmtr; "^2"
PRINT #1, " w1="; w1; "/c"
PRINT #1, " Re1="; Re1
PRINT #1, " Pr1="; Pr1
PRINT #1, " Nu1="; Nu1
PRINT #1, " alf1="; alf1; "/2*"
|
|
|
PRINT #1, " k="; k1; "/^2*"
PRINT #1, " F="; F; "^2"
PRINT #1, " L="; ltr; ""
END
:
- 35 ,
- 34.07678 ,
- 34.53633 ,
tc2= 137.4979 ,
tc1= 102.9616 ,
3458.24 3481.483 /*
.1185689.1258007 /*
. 19.96261 20.46453 **10^6
Q= 1.898105E+07
alf2
n= 709
ug=.4064007 ^3/c
ro2= 92.51952 /^3
ftr=.78525 ^2
w2= 2.852322 /c
Re2= 206324.1
Pr2=.5663478
Nu2= 293.5521
alf2= 2308.065 /2*
alf1
n= 709
uh=.3336158 ^3/c
ro1= 112.7045 /^3
fmtr=.134 ^2
w1= 2.48967 /c
Re1= 351403.3
Pr1=.5822394
Nu1= 437.3047
alf1= 2074.03 /2*
k= 987.9561 /^2*
F= 556.2966 ^2
L= 10
2.
CLS
'
sigt = 240: sigv = 385
'
nt = 1.5: nv = 2.4
'
P = 30
'
fi =.9
'
C1 =.002: C2 =.0012: C3 = 0
'
d = 1
sig1 = sigt / nt: sig2 = sigv / nv
IF sig1 > sig2 THEN sig = sig1 ELSE sig = sig2
btr = EXP(P / (sig * fi)): SS =.5 * d * (btr - 1) + C1
'PRINT "SS="; SS
'INPUT " SS,"; SS
A = (C2 + C3) / SS
IF A >.05 THEN C = C1 + C2 + C3 ELSE C = C1
S =.5 * d * (btr - 1) + C
PRINT " [sig]="; sig; ""
PRINT " ="; C; ""
PRINT " S="; S; ""
INPUT " S,"; S
bt = 1 + 2 * (S - C) / d
PRINT " bt="; bt
Pd = sig * fi * LOG(bt)
'PRINT " , S="; S
'INPUT " S"; S
PRINT " ,[P]="; Pd
END
[sig]= 160.4167
=.002
S=.1369645
S,?.14
bt= 1.276
,[P]= 35.18855
3.
CLS
'
P = 30
'
D = 1.6
'
ds =.1
'
z = 24
'
sig = 160
'
C =.002
'
fi = 1
'
|
|
|
S =.14
'
lb =.12
A = D + 2 * ds: B = 2.2 * ds / SIN(3.14 / z)
IF A > B THEN Db = A ELSE Db = B
PRINT "Db>="; Db; ""
INPUT " Db,"; Db
PRINT "D2>="; Db + 2 * ds; ""
INPUT " D2,"; D2
btr = EXP(P / (sig * fi))
PRINT " btr="; btr
Sfr =.5 * D * (btr - 1)
PRINT " Sfr="; Sfr; ""
PRINT "rn>="; Sfr; ""
INPUT " rn,"; rn
'INPUT "Sfr"; sfr
IF Sfr + C > S THEN Sf = S ELSE Sf = S
PRINT "Sf="; Sf; ""
INPUT " Sf,"; Sf
bt = (D + 2 * (Sf - C)) / D
PRINT " bt="; bt
Pd = sig * fi * LOG(bt)
PRINT " [p]="; Pd; ""
h3 = lb +.25 * ds
PRINT "h3>="; h3; ""
INPUT " h3,"; h3
h4 =.7 * Sfr
PRINT "h4>="; h4; ""
INPUT " , h4,"; h4
END
Db>= 1.686335
Db,? 1.7
D2>= 1.9
D2,? 1.9
btr= 1.240637
Sfr=.1203185
rn>=.1203185
rn,?.12
Sf=.14
Sf,?.14
bt= 1.276
[p]= 38.99683
h3>=.145
h3,?.145
h4>= 8.422297E-02
, h4,?.085
4.
CLS
'
DR = 1.1
'
p = 34.5
'
sig = 160
'
Db = 1.7
'
F = 36.39
'
FQ = 32.77
'
D4 = 1.7
'
d0 =.11
'
di =.22
'
C =.002
H3P =.41 * DR * SQR(p / sig * (DR + 3 * (Db - DR) * F / FQ) / (D4 - 2 * d0 - di))
PRINT " H3P"; H3P; ""
HS3 = H3P + C
PRINT " HS3="; HS3; ""
INPUT "HS3,"; HS3
INPUT " D5,"; D5
INPUT " D6,"; D6
INPUT " H1,"; H1
H2R =.41 * DR * SQR(p / sig * (3 * (Db - DR) * F / FQ) / DR + 1.5 * (p / sig) ^ 2 - 6 * H3 ^ 2 * (D4 - DR - 2 * d0) / DR)
PRINT " H2="; H2R + C; ""
INPUT "H2,"; H2
H3R =.41 * DR * SQR(p / sig * (DR + 3 * (Db - DR) * F / FQ) / ((D4 - 2 * d0 - D6) + ((H2 / HS3) ^ 2) * (D6 - D5) + ((H1 / HS3) ^ 2) * (D5 - di)))
PRINT " H3="; H3R + C; ""
|
|
|
INPUT "H3,"; H3
END
H3P.3284282
HS3=.3304282
HS3,?.33
D5,? 1.08
D6,? 1.15
H1,?.35
H2=.3084011
H2,?.31
H3=.3182192
H3,?.33
5.
CLS
'
sigts = 315: sigBS = 530
'
sigBK = 385
'
D = 1
'
P = 34.5
'
Dr = 1.1
'
nt = 1.5
'
H0 =.012
'
gm = 45 * 3.14 / 180
'
h1 =.07: h2 =.035
'
z = 24
IF P < 40 THEN PF = 40
PRINT " PF="; PF; ""
hr = (h1 + h2) / 2
PRINT " hr="; hr; ""
Fp =.5 * 3.14 * PF * Dr * hr * TAN(gm)
Fq = 3.14 * Dr ^ 2 / 4 * P
F = Fp + Fq
PRINT " F="; F; ""
PRINT " Fp="; Fp; ""
PRINT " Fq="; Fq; ""
INPUT " K1s="; K1s
INPUT " K2s="; K2s
sigs = sigts / nt
PRINT " sigs="; sigs; ""
d1sr = SQR(4 * K1s * K2s * F / (3.14 * z * sigs) + d0s ^ 2)
PRINT " d1sr="; d1sr; ""
INPUT " d2s="; d2s
INPUT " d3s="; d3s
PRINT ".5*(d2s+d3s)=";.5 * (d2s + d3s); ""
INPUT " d1s="; d1s
INPUT " ds="; ds
IF sigBS / sigBK > 1.5 THEN lr =.5 * ds * (1 + sigBS / sigBK) ELSE lr =.5 * ds
END
PF= 40
hr=.0525
F= 36.39364
Fp= 3.623814
Fq= 32.76983
K1s=? 1
K
2s=? 1.1
sigs= 210
d1sr=.100591
d2s=?.0961
d3s=?.935
.5*(d2s+d3s)=.51555
d1s=?.096
ds=?.1
6. 䳿
CLS
' ,
L = 10
' ,
a = 1
' ,
G =.8
'
n = 3
' ,
l1 = 4
'
fi = 1
' ,
D = 1
' ,
dlt1 = 120 * 3.1415 / 180
' ,
S =.14
' ,
C =.0025
'
nc = 2.4
' ,
dlt2 = 150 * 3.141 / 180
' ,
b =.22
' ,
sig = 160
' ,
E = 199000
' ,
b2 =.45
'
k2 = 1.2
'
nt = 1.5
' , /
q = G / (L + 4 / 3 * H)
'
psi1 =.9: psi2 = 1.2
PRINT " q ="; q; "/"
F1 = psi1 * G / n: F2 = psi2 * G / n
|
|
|
IF F1 > F2 THEN Fbi = F1 ELSE Fbi = F2
PRINT " Fbi ="; Fbi; ""
M = q * l1 ^ 2 / 8
PRINT " ="; M; "*"
Qi =.5 * fi
PRINT " Qi ="; Qi; ""
IF fm1 > SQR(D * (S - C)) THEN fi = 1 ELSE INPUT "fi"; fi
gm = 2.83 * a / D * SQR((S - C) / D)
PRINT " gm ="; gm
bt =.91 * b / SQR(D * (S - C))
PRINT " bt ="; bt
sigmx = 4 * M / (3.14 * D ^ 2 * (S - C))
'PRINT " sigmx = "; sigmx; ""
'Pdm = 2 * sig * (S - C) / (D + S - C)
'PRINT " Pdm = "; Pdm; ""
'PRINT " :"
'Mdm = 3.14 / 4 * D * (D + S - C) * (S - C) * sig
'PRINT " Mdm = "; Mdm; "*"
'Mde =.000089 * E / nc * D ^ 3 * (100 * (S - C) / D) ^ 2.5
'PRINT " Mde = "; Mde; "*"
'Md = Mdm / SQR(1 + (Mdm / Mde) ^ 2)
'PRINT " Md="; Md; "*"
'PRINT " :"
'Fdm = 3.14 * (D - S - C) * (S - C) * sig
'PRINT " Fdm = "; Fdm; ""
'Fde =.00031 * E / nc * D ^ 2 * (100 * (S - C) / D) ^ 2.5
'PRINT " Mde = "; Mde; ""
'Fd = Fdm / SQR(1 + (Fdm / Fde) ^ 2)
'PRINT " Fd="; Fd; ""
'PRINT " :"
'Qdm =.25 * 3.14 * D * (S - C) * sig
'PRINT " Qdm = "; Qdm; ""
'Qde = 2.4 * E * (S - C) ^ 2 / nc * (.18 + 3.3 * D * (S - C) / L ^ 2)
'PRINT " Qde = "; Qde; ""
'Qd = Qdm / (1 + (Qdm / Qde) ^ 2)
'PRINT " Qd="; Qd; ""
'y = D / (S - C): u = L / D
'PRINT " y,u"; y; u
k10 = EXP(-3.58) * SIN(3.58) / 3.58
IF k10 <.25 THEN k10 =.25
k11 = (1 - EXP(-bt) * COS(bt)) / 3.58
k12 = (1.15 -.1432 * dlt1) / SIN(.5 * dlt1)
IF 1.7 - 2.1 * dlt1 / 3.14 < 0 THEN k131 = 0 ELSE k131 = 1.7 - 2.1 * dlt1 / 3.14
k13 = k131 / SIN(.5 * dlt1)
k14 = (1.45 -.43 * dlt1) / SIN(.5 * dlt1)
k15 = (.8 * SQR(gm) + 6 * gm) / dlt1
IF k15 > 1 THEN k15 = 1
k16 = 1 -.65 / (1 + (6 * gm) ^ 2) * SQR(3.14 / (3 * dlt1))
k17 = 1 / (1 +.6 * (D / (S - C)) ^ (1 / 3) * b / D * dlt1)
k19 = 5 / (6 * (D / (S - C)) ^ (1 / 3) * dlt1)
Lon = (D + 2 * S) * dlt1 / 2
PRINT " dlt1,Lon="; Lon
Lnn = Lon * dlt2 / dlt1
PRINT " dlt2,Lnn="; Lnn
fm = (Lnn - Lon) / 2
PRINT " ";
IF (k19 * D + 1.5 * b) <= b2 THEN PRINT "" ELSE PRINT ""
Sef = (S - C) * SQR(1 + (S / (S - C)) ^ 2)
gmn = 2.83 * a / D * SQR((Sef - C) / D)
btn =.91 * b2 / ((Sef - C) * D)
sigmxn = 4 * M / (3.14 * D ^ 2 * (Sef - C))
lm = L
Pdmn = 2 * sig * (Sef - C) / (D + Sef - C)
B1n = 9.45 * D / lm * SQR(D / (100 * (Sef - C)))
IF B1n > 1 THEN B1n = 1
Pden =.0000208 * E * D / (nc * B1n * lm) * ((100 * (Sef - C)) / D) ^ 2.5
PRINT " Pden = "; Pden; ""
PRINT " Pdmn = "; Pdmn; ""
Pdn = Pdmn / (1 + (Pdmn / Pden) ^ 2)
PRINT " Pdn = "; Pdn; ""
PRINT " Mdm = "; Mdm; "*"
Mdmn = 3.14 / 4 * D * (D + Sef - C) * (Sef - C) * sig
PRINT " Mdmn = "; Mdmn; "*"
Mden =.000089 * E / nc * D ^ 3 * (100 * (Sef - C) / D) ^ 2.5
PRINT " Mden = "; Mden; "*"
Mdn = Mdmn / (1 + (Mdmn / Mden) ^ 2)
PRINT " Mdn="; Mdn; "*"
PRINT " :"
Fdmn = 3.14 * (D - Sef - C) * (Sef - C) * sig
PRINT " Fdmn = "; Fdmn; ""
Fden =.00031 * E / nc * D ^ 2 * (100 * (Sef - C) / D) ^ 2.5
PRINT " Fden = "; Fden; ""
Fdn = Fdmn / (1 + (Fdmn / Fden) ^ 2)
PRINT " Fdn="; Fdn; ""
PRINT " :"
Qdmn =.25 * 3.14 * D * (Sef - C) * sig
PRINT " Qdmn = "; Qdmn; ""
Qden = 2.4 * E * (Sef - C) ^ 2 / nc * (.18 + 3.3 * D * (Sef - C) / lm ^ 2)
PRINT " Qden = "; Qden; ""
Qdn = Qdmn / (1 + (Qdmn / Qden) ^ 2)
PRINT " Qdn="; Qdn; ""
yn = D / (Sef - C): u = lm / D
PRINT " y,u"; yn; u
k10n = EXP(-btn) * SIN(btn) / btn
IF k10n <.25 THEN k10n =.25
k11n = (1 - EXP(-btn) * COS(btn)) / btn
k12n = (1.15 -.1432 * dlt2) / SIN(.5 * dlt2)
IF 1.7 - 2.1 * dlt2 / 3.14 < 0 THEN k131 = 0 ELSE k131 = 1.7 - 2.1 * dlt2 / 3.14
k13n = k131 / SIN(.5 * dlt1)
k14n = (1.45 -.43 * dlt2) / SIN(.5 * dlt2)
k15n = (.8 * SQR(gmn) + 6 * gmn) / dlt2
IF k15n > 1 THEN k15n = 1
k16n = 1 -.65 / (1 + (6 * gmn) ^ 2) * SQR(3.14 / (3 * dlt2))
k17n = 1 / (1 +.6 * (D / (Sef - C)) ^ (1 / 3) * b2 / D * dlt2)
PRINT " ";
IF P * D / (4 * (Sef - C)) <= sig * fi THEN PRINT " " ELSE PRINT " "
v21n2 = -sigmxn * k2 / (nt * sig)
v22n2 = ((P * D / (4 * (Sef - C))) - sigmxn) * k2 / (nt * sig)
v21n4 = 0
v22n4 = P * D / (2 * (Sef - C)) * k2 / (nt * sig)
IF v21n2 >= 0 THEN v11n2 = -.23 * k13n * k15n / (k12n * k10n) ELSE v11n2 =.23 * k13n * k15n / (k12n * k10n)
IF v22n2 >= 0 THEN v12n2 = -.23 * k13n * k15n / (k12n * k10n) ELSE v12n2 =.23 * k13n * k15n / (k12n * k10n)
IF v21n4 >= 0 THEN v11n4 = -.53 * k11n / (k14n * k16n * k17n * SIN(.5 * dlt2)) ELSE v11n4 =.53 * k11n / (k14n * k16n * k17n * SIN(.5 * dlt2))
IF v22n4 >= 0 THEN v12n4 = -.53 * k11n / (k14n * k16n * k17n * SIN(.5 * dlt2)) ELSE v12n4 =.53 * k11n / (k14n * k16n * k17n * SIN(.5 * dlt2))
v21n2 = ABS(v21n2)
v22n2 = ABS(v22n2)
v22n4 = ABS(v22n4)
v2 = k2 / nt * sigmxn / (sig * fi)
PRINT ", "
PRINT "v11n2="; v11n2
PRINT "v21n2="; v21n2
PRINT ", "
PRINT "v11n4="; v11n4
PRINT "v21n4="; v21n4
INPUT " k1n2"; k1n2
INPUT " k1n4"; k1n4
PRINT " "
sigin2 = k1n2 * sig * nt / k2
sigin4 = k1n4 * sig * nt / k2
PRINT " 2 sigin2 ="; sigin2; ""
PRINT " 4 sigin4 ="; sigin4; ""
Fdn2 =.7 * sigin2 * SQR(D * (Sef - C)) * (Sef - C) / (k10n * k12n)
Fdn4 =.9 * sigin4 * SQR(D * (Sef - C)) * (Sef - C) / (k14n * k16n * k17n)
PRINT " Fdn2 ="; Fdn2; ""
PRINT " Fdn4 ="; Fdn4; ""
Fen = Fbi * 3.14 / 4 * SQR(D / (Sef - C)) * k13n * k14n
PRINT " Fen = "; Fen; ""
PRINT "Fbi="; Fbi
PRINT "Fdn2="; Fdn2
PRINT "Fdn4="; Fdn4
PRINT " ";
IF Fen <= Fdn2 AND Fen <= Fdn4 THEN PRINT "" ELSE PRINT ""
PRINT "M / Mdn + Fen / Fdn + (Qi / Qdn) ^ 2= "; M / Mdn + Fen / Fdn + (Qi / Qdn) ^ 2
PRINT " ";
IF M / Mdn + Fen / Fdn + (Qi / Qdn) ^ 2 <= 1 THEN PRINT "" ELSE PRINT ""
END
q =.08 /
Fbi =.32
=.16 *
Qi =.5
fi? 1
gm = 1.049392
bt =.5398993
dlt1,Lon= 1.340373
dlt2,Lnn= 1.6752
Pden = 1326.981
Pdmn = 51.9327
Pdn = 51.85329
Mdm = 0 *
Mdmn = 29.04644 *
Mden = 12190.59 *
Mdn= 29.04627 *
:
Fdmn = 77.9876
Fden = 42461.59
Fdn= 77.98734
:
Qdmn = 24.3325
Qden = 1392.121
Qdn= 24.32507
y,u 5.16182 10
,
v11n2= 0
v21n2= 5.260454E-03
,
v11n4=-1.836972
v21n4= 0
k1n2? 1.477
k1n4?.387
2 sigin2 = 295.4
4 sigin4 = 77.39999
Fdn2 = 87.87797
Fdn4 = 39.56281
Fen = 0
Fbi=.32
Fdn2= 87.87797
Fdn4= 39.56281
M / Mdn + Fen / Fdn + (Qi / Qdn) ^ 2= 5.930956E-03
