1 CLS
PRINT "LABORATORIJ MRA"
PRINT "PREHRAMBENO BIOTEHNOLO[KI FAKULTET"
PRINT "------------------------------------"
PRINT
PRINT "MATEMATI^KO MODELIRANJE I VO\ENJE INDUSTRIJSKIH PROCESA"
PRINT "  KONCENTRIRANJE SOKA NARAN^E S TRI STUPNJA UPARAVANJA"
PRINT : PRINT "Upi{ite slijede}e podatke"
PRINT
10 INPUT "PROTOK SOKA OD NARAN^E NA ULAZU U PROCES (kg/h) ="; Qso
IF Qso < 15000 THEN PRINT "Q>=15000": GOTO 10
IF Qso > 20000 THEN PRINT "Q<=20000": GOTO 10
20 INPUT "PROTOK PARE U PRVI IZMJENJIVA^ TOPLINE   (kg/h) ="; Qp
IF Qp < 4000 THEN PRINT "Q>=4000": GOTO 20
IF Qp > 6200 THEN PRINT "Q<=6200": GOTO 20
30 INPUT "MASENI UDIO SUHE TVARI U SOKU NA ULAZU   (kg/kg) ="; xo
IF xo < .05 OR xo > .12 THEN PRINT "0.05 <= x <= 0.12": GOTO 30
40 INPUT "TEMPERATURA SOKA NA ULAZU                    (C) ="; Tso
IF Tso < 5 OR Tso > 25 THEN PRINT "5 <= T(C) <= 25 ": GOTO 40
50 INPUT "TEMPERATURA PARE NA ULAZU U PRVI IZMJENJIVA^ (C) ="; Tpo
IF Tpo < 115 OR Tpo > 125 THEN PRINT "115 >= T(C) <= 125": GOTO 50
CLS
PRINT "ZADANI SU SLIJEDE]I PODACI"
PRINT "-----------------------------"
PRINT "PROTOK SOKA OD NARAN^E NA ULAZU U PROCES  (kg/h)  ="; Qso
PRINT "PROTOK PARE U PRVI IZMJENJIVA^ TOPLINE    (kg/h)  ="; Qp
PRINT "MASENI UDIO SUHE TVARI U SOKU NA ULAZU    (kg/kg) ="; xo
PRINT "TEMPERATURA SOKA NA ULAZU                    (C)  ="; Tso
PRINT "TEMPERATURA PARE NA ULAZU U PRVI IZMJENJIVA^ (C)  ="; Tpo
PRINT
PRINT
INPUT "Da li su podaci ispravno zadani (*/n)"; r$
IF r$ = "n" OR r$ = "N" THEN GOTO 1
Qpul = Qp
Qsul = Qso
S = 266
ho = 1576 * 3.6
h1 = 1091 * 3.6
DEF fnH (x) = ho - h1 * x
co = 4.1865
c1 = 4.1865 - 1.491
DEF fncp (x) = co - c1 * x
io = 2519.74
i1 = 1.584
DEF fnip (T) = io + i1 * T
i2 = 4.1886
DEF fniv (T) = i2 * T

CLS
' prvi stupanj
alfao = Qso * fncp(xo) * Tso + Qp * fnip(Tpo)
bo = (ho * S * Tpo - Qp * fnip(Tpo)) / h1 / S
b1 = (Qp * i2 - ho * S) / h1 / S
naz = i1 - xo * c1 + Qp / Qso * i2
d0 = -io / (co - i1)
d1 = xo * (co - i1) / naz
d2 = (io - alfao / Qso) / naz
p1 = .5 * (d1 * Tpo - d0 * d1 + bo + d2 * b1) / (b1 - d1)
p2 = (d1 * d0 * Tpo + d2 * bo) / (b1 - d1)
dis = p1 * p1 - p2
IF dis < 0 THEN PRINT "dis<0": END
T = -p1 + SQR(dis)
x = (bo + b1 * T) / (Tpo - T)
QL = xo / x * Qso
QV = Qso - QL
CLS
PRINT "            BILANCE MASE I ENERGIJE ZA PRVI UPARIVA^"
PRINT "-------------------------------------------------------------"
PRINT "1. UKUPNA BILANCA MASE"
PRINT "Ulazni  protok soka naran~e            (kg/h) ="; Qso
PRINT "Izlazni protok koncentrata             (kg/h) ="; QL
PRINT "Izlazni protok para soka naran~e       (kg/h) ="; QV
PRINT "Ulazni  protok zasi~ene vodene pare    (kg/h) ="; Qp
PRINT
PRINT "2. BILANCA MASE SUHE TVARI"
PRINT "Maseni udio suhe tvari u ulaznom  toku (kg/kg) ="; xo
PRINT "Maseni udio suhe tvari u izlaznom toku (kg/kg) ="; x
PRINT
PRINT "3. BILANCA ENERGIJE"
PRINT "Ulazni  tok entalpije soka naran~e          (MJ/h) ="; Qso * fncp(xo) * Tso / 1000
PRINT "Ulazni  tok entalpije zasi~ene vodene pare  (MJ/h) ="; Qp * fnip(Tpo) / 1000
PRINT "Izlazni tok entalpije para naran~inog soka  (MJ/h) ="; QV * fnip(T) / 1000
PRINT "Izlazni tok entalpije koncentrata soka      (MJ/h) ="; QL * fncp(x) * T / 1000
PRINT "Izlazni tok entalpije ukapljene vodene pare (MJ/h) ="; Qp * fniv(T) / 1000

PRINT
PRINT "TEMPERATURA U UPARIVA^U T(C) ="; T
PRINT "EFIKASNOST UPARIVANJA (kg otparenog soka/ kg vodene pare )="; QV / Qp
PRINT : PRINT
INPUT "drugi stupanj"; r$

CLS
' drugi stupanj
Qso = QL
xo = x
Tso = T
Qp = QV
Tpo = T

alfao = Qso * fncp(xo) * Tso + Qp * fnip(Tpo)
bo = (ho * S * Tpo - Qp * fnip(Tpo)) / h1 / S
b1 = (Qp * i2 - ho * S) / h1 / S
naz = i1 - xo * c1 + Qp / Qso * i2
d0 = -io / (co - i1)
d1 = xo * (co - i1) / naz
d2 = (io - alfao / Qso) / naz
p1 = .5 * (d1 * Tpo - d0 * d1 + bo + d2 * b1) / (b1 - d1)
p2 = (d1 * d0 * Tpo + d2 * bo) / (b1 - d1)
dis = p1 * p1 - p2
IF dis < 0 THEN PRINT "dis<0": END
T = -p1 + SQR(dis)
x = (bo + b1 * T) / (Tpo - T)
QL = xo / x * Qso
QV = Qso - QL
CLS
PRINT "            BILANCE MASE I ENERGIJE ZA DRUGI UPARIVA^"
PRINT "-------------------------------------------------------------"
PRINT "1. UKUPNA BILANCA MASE"
PRINT "Ulazni  protok soka naran~e            (kg/h) ="; Qso
PRINT "Izlazni protok koncentrata             (kg/h) ="; QL
PRINT "Izlazni protok para soka naran~e       (kg/h) ="; QV
PRINT "Ulazni  protok zasi~ene vodene pare    (kg/h) ="; Qp
PRINT
PRINT "2. BILANCA MASE SUHE TVARI"
PRINT "Maseni udio suhe tvari u ulaznom  toku (kg/kg) ="; xo
PRINT "Maseni udio suhe tvari u izlaznom toku (kg/kg) ="; x
PRINT
PRINT "3. BILANCA ENERGIJE"
PRINT "Ulazni  tok entalpije soka naran~e          (MJ/h) ="; Qso * fncp(xo) * Tso / 1000
PRINT "Ulazni  tok entalpije pare soka naran~e     (MJ/h) ="; Qp * fnip(Tpo) / 1000
PRINT "Izlazni tok entalpije para naran~inog soka  (MJ/h) ="; QV * fnip(T) / 1000
PRINT "Izlazni tok entalpije koncentrata soka      (MJ/h) ="; QL * fncp(x) * T / 1000
PRINT "Izlazni tok entalpije ukapljene pare soka   (MJ/h) ="; Qp * fniv(T) / 1000

PRINT
PRINT "TEMPERATURA U UPARIVA^U T(C) ="; T
PRINT "EFIKASNOST UPARIVANJA (kg otparenog soka/ kg vodene pare )="; QV / Qp
PRINT : PRINT
INPUT "tre}i stupanj"; r$

CLS
Qso = QL
xo = x
Tso = T
Qp = QV
Tpo = T

alfao = Qso * fncp(xo) * Tso + Qp * fnip(Tpo)
bo = (ho * S * Tpo - Qp * fnip(Tpo)) / h1 / S
b1 = (Qp * i2 - ho * S) / h1 / S
naz = i1 - xo * c1 + Qp / Qso * i2
d0 = -io / (co - i1)
d1 = xo * (co - i1) / naz
d2 = (io - alfao / Qso) / naz
p1 = .5 * (d1 * Tpo - d0 * d1 + bo + d2 * b1) / (b1 - d1)
p2 = (d1 * d0 * Tpo + d2 * bo) / (b1 - d1)
dis = p1 * p1 - p2
IF dis < 0 THEN PRINT "dis<0": END
T = -p1 + SQR(dis)
x = (bo + b1 * T) / (Tpo - T)
QL = xo / x * Qso
QV = Qso - QL
CLS
PRINT "            BILANCE MASE I ENERGIJE ZA TRE]I UPARIVA^"
PRINT "-------------------------------------------------------------"
PRINT "1. UKUPNA BILANCA MASE"
PRINT "Ulazni  protok soka naran~e            (kg/h) ="; Qso
PRINT "Izlazni protok koncentrata             (kg/h) ="; QL
PRINT "Izlazni protok para soka naran~e       (kg/h) ="; QV
PRINT "Ulazni  protok zasi~ene vodene pare    (kg/h) ="; Qp
PRINT
PRINT "2. BILANCA MASE SUHE TVARI"
PRINT "Maseni udio suhe tvari u ulaznom  toku (kg/kg) ="; xo
PRINT "Maseni udio suhe tvari u izlaznom toku (kg/kg) ="; x
PRINT
PRINT "3. BILANCA ENERGIJE"
PRINT "Ulazni  tok entalpije soka naran~e          (MJ/h) ="; Qso * fncp(xo) * Tso / 1000
PRINT "Ulazni  tok entalpije pare soka naran~e     (MJ/h) ="; Qp * fnip(Tpo) / 1000
PRINT "Izlazni tok entalpije para naran~inog soka  (MJ/h) ="; QV * fnip(T) / 1000
PRINT "Izlazni tok entalpije koncentrata soka      (MJ/h) ="; QL * fncp(x) * T / 1000
PRINT "Izlazni tok entalpije ukapljene pare soka   (MJ/h) ="; Qp * fniv(T) / 1000
PRINT
PRINT "TEMPERATURA U UPARIVA^U T(C) ="; T
PRINT "EFIKASNOST UPARIVANJA (kg otparenog soka/ kg vodene pare )="; QV / Qp
PRINT
PRINT "UKUPNA EFIKASNOST UPARIVANJA (kg otparenog soka/kg vodene pare)="; (Qsul - QL) / Qpul

END