CLS
10 REM
COLOR 2: PRINT "LABORATORIJ MRA "
PRINT "PREHRAMBENO BIOTEHNOLOSKI FAKULTET U ZAGREBU "
PRINT "_______________________________________________"
PRINT "DINAMIKA SMRZAVANJA PREHRAMBENIH PROIZVODA"
PRINT "-----------------------------------------------"
INPUT "Upi{ite naziv prehrambenog proizvoda:"; naziv$
PRINT "Upi{ite dimenzije proizvoda pakiranog u obliku kvadra: "
INPUT " kolika je visina kvadra (maksimalna vrijednost=0.5 m): D(m) ="; visina
IF d > .5 GOTO 10
INPUT " kolika je {irina kvadra (neka je maksimalno 0.1*D )   : L(m) ="; sirina
IF L > .125 GOTO 10
INPUT "Upi{ite naziv ambala`nog materijala "; ambalaza$
INPUT "Upi{ite debljinu ambala`e d(m) ="; damb
PRINT "Po~etna temperatura materijala ( 10<=T<=30 C )"; : INPUT Tpoc
IF Tpoc > 30 OR Tpoc < 10 THEN GOTO 10
PRINT "Upi{ite temperaturu zraka za hla|enje  od 0 C do -50 C:"
INPUT " temperatura zraka To (C) ="; t0
IF t0 > 0 OR t0 < -50 GOTO 10
PRINT "Upi{ite temperaturu na kojoj proizvod smrzava Ts (C),"
IF Ts > 0 OR Ts < -50 GOTO 10
PRINT "Ako Vam nije poznata vrijednost iz literature upi{ite proizvoljnu "
PRINT "vrijednost, napr. -5 C"
INPUT " Temperatura smrzavanja Ts (C) ="; Ts
INPUT "Upi{ite brzinu rashladnog zraka v (m/s) ="; brzina
PRINT "Upi{ite vrijednost koeficijenta toplinske vodljivosti,"
PRINT "prehrambenog proizvoda u nesmrznutom stanju: "
INPUT " k-proizvoda (W/mK) ="; kproizvodat
PRINT "Upi{ite vrijednost specifi~ne topline nesmrznutog proizvoda:";
INPUT " Cp (J/kg) = "; cpp
PRINT "Upi{ite gusto}u nesmrznutog proizvoda"
PRINT "(ako ne raspola`ete s podatkom upi{ite vrijednost 1000)"
INPUT " gusto}a u nesmrznutom stanju ro (kg/m3) ="; gustocap
PRINT "Upi{ite maseni udio vode (x) u proizvodu  (0< x <1)"
INPUT " Maseni udio vode x ="; x
IF x < 0 OR x > 1 GOTO 10
PRINT "Upi{ite koeficijent toplinske vodljivosti ambala`e:"
INPUT " k-ambala`e (W/mK) ="; kambalaze
cpsuhet = cpp - x * 4180: IF cpsuhet < 0 THEN cpsuhet = 0
cpps = cpsuhet + x * 4180 * .48
rosmrz = gustocap * (1 - x * .1)
CLS
PRINT "SMRZAVANJE "; naziv$
PRINT "-------------------------------"
PRINT "DIMENZIJE PAKETA KOJI IMA OBLIK KVADRA"
PRINT " Visina (m) ="; visina
PRINT " [irina (m) ="; sirina
PRINT "PODACI ZA AMBALA@U "; ambalaza$
PRINT " Debljina ambala`e (m)="; damb
PRINT " Koeficijent toplinske vodljivosti ambala`e k (W/mK) ="; kambalaze
PRINT "Proizvod se smrzava zrakom temperature To (C) ="; t0
PRINT "Brzina zraka za hla|enje v (m/s)="; brzina
PRINT "Temperatura na kojoj proizvod smrzava Ts (C)="; Ts
PRINT "Koeficijent toplinske vodljivosti proizvoda"
PRINT "u nesmrznutom stanju k (W/mK) ="; kproizvodat
PRINT "Gusto}a proizvoda u nesmrznutom stanju ro (kg/m3) ="; gustocap
PRINT "Gustoca smrznutog proizvoda ro (kg/m3)="; rosmrz
PRINT "Specifi~na toplina nesmrznutog proizvoda Cp (J/kg) ="; cpp
PRINT "Specificna toplina smrznutog proizvoda Cp (J/kg) ="; cpps
PRINT "Maseni udio vode u proizvodu je x ="; x
PRINT "-----------------------------------------"
INPUT "da li su podaci ispravno upisani (da/ne) ="; r$
IF r$ = "ne" GOTO 10
CLS
PRINT "Odre|ivanje koeficijenta prijenosa topline u hidrodinami~kom"
PRINT "grani~nom sloju"
' proracun Reynoldsovog broja
mizraka = (17.19 + .0518 * t0) * .000001
kzraka = .02373 + .0000642# * T
cp = 1004.8
tlak = 101325!
r = 8.314
m = .029
rozraka = tlak * m / r / (273 + t0)
PRINT "Gusto}a rashladnog zraka na temperaturi To (C) "; t0
PRINT "            ro (kg/m3) ="; rozraka
PRINT "Viskozitet zraka  mi (Pas) ="; mizraka
PRINT "Koeficijent toplinske vodljivosti zraka k (W/mK) ="; kzraka
PRINT "Specifi~na toplina zraka Cp (J/kgK) ="; cp
Re = rozraka * brzina * visina / mizraka
PRINT "Reynoldsov broj Re ="; Re;
IF Re < 3500 THEN PRINT "  laminarno strujanje"
IF Re > 3500 THEN PRINT "  turbulentno strujanje"
Pr = mizraka * cp / kzraka
PRINT "Prandtlov broj Pr ="; Pr
IF Re < 3500 THEN Nu = .664 * Re ^ .5 * Pr ^ (1 / 3)
IF Re > 3500 THEN Nu = .027 * Re ^ .8 * Pr ^ .33
PRINT "Nusseltov broj Nu ="; Nu
hg = Nu * kzraka / visina
PRINT "Koeficijent prijenosa topline u grani~nom sloju  hg (W/m2K) ="; hg
PRINT "------------------------------"
PRINT "Ako ste prepisali rezultate pritisnite tipku ENTER "
INPUT r$
CLS
PRINT "USPOREDBA OTPORA PRIJENOSU TOPLINE U (m2K/W)"
hstar = 3000
r1 = 1 / hg
r2 = damb / kambalaze
r3 = 1 / hstar
r4 = .1 / 2.442
PRINT " Otpor u grani~om sloju zraka "; r1
PRINT " Otpor ambala`e ="; r2

PRINT " Otpor na granici smrznuto/nesmrznuto ="; r3
PRINT " Otpor leda debljine 0.1 (m) na -20 C ="; r4
PRINT "Otpor smrznutog dijela proizvoda raste proporcionalno debljini"
PRINT "zamrznutog proizvoda"
PRINT
PRINT "Ako ste prepisali podatke pritisnite tipku ENTER"
INPUT r$
CLS
PRINT "Koeficijent toplinske vodljivosti i gusto}a proizvoda"
PRINT "na temperaturi smrzavanja Ts (C)="; Ts
kleda = 2.21 - .0114 * Ts
kvode = .588
kpsmrz = kproizvodt + (kleda - kvode) * x
ropsmrz = gustocap - 100 * x
Hl = 335000!
PRINT " ks (W/mK)="; kpsmrz
PRINT " ro (kg/m3)="; ropsmrz
PRINT
PRINT "Ako ste prepisali podatake pritisnite ENTER"
INPUT r$
CLS

' dinamika hladjenja
DIM tos(65), yos(65)
PRINT " DINAMIKA HLADJENJA DO TEMPERATURE SMRZAVANJA "
tautekuce = gustocap * cpp * sirina / 2 / hg
tausmrz = rosmrz * cpps * sirina / 2 / hg
thladt = tautekuce * LOG((Tpoc - t0) / (Ts - t0))
thlads = tausmrz * LOG((Ts - t0) / 5)
time = 0
delt = thladt / 20
FOR i = 1 TO 21
tmp = t0 + (Tpoc - t0) * EXP(-time / tautekuce)
tos(i) = time: yos(i) = tmp
PRINT "vrijeme(min)="; time / 60, " T( C )="; tmp
time = time + delt
NEXT i
thladf = time - delt
ifin = 21
PRINT "--------------------------------------------"
30 INPUT "Da li ste prepisali podatke ?  (da/*) "; r$
IF r$ <> "da" GOTO 30
CLS

PRINT "DINAMIKA SMRZAVANJA VODE U MATERIJALU"
alfa = (Ts - t0) / ropsmrz / x / Hl
beta = r1 + r2 + r3
FOR i = 1 TO 21
del = (i - 1) * sirina / 20 / 2
time = thladf + (beta * del + .5 / kpsmrz * del * del) / alfa
tos(ifin + i) = time: yos(ifin + i) = Ts
PRINT "Vrijeme (min)="; time / 60, "Debljina smrznutog sloja (mm)="; del * 1000
NEXT i
20 INPUT "Da li ste prepisali podatke ?  (da/*) "; r$
IF r$ <> "da" GOTO 20

CLS
PRINT "DINAMIKA HLADJENJA SMRZNUTOG PROIZVODA "
if2 = ifin + 21
del = thlads / 20
FOR i = 1 TO 21
j = if2 + i
time = time + del
tmp = t0 + (Ts - t0) * EXP(-i * del / tausmrz)
tos(j) = time: yos(j) = tmp
PRINT "Vrijeme (min) ="; time / 60, "T (C) ="; tmp
NEXT i
PRINT "-------------------------------------------"
PRINT "Ako ste prepisali podatke pritisnite ENTER za grafi~ki prikaz";
INPUT r$
tmax = time / 60: PRINT "tmax="; tmax

FOR i = 1 TO 63
tos(i) = tos(i) / 60 / tmax * 650
yos(i) = (yos(i) - 30) / (-80) * 250
NEXT i
REM  graf
SCREEN 9
CLS
PSET (50, 0)
DRAW "D250": DRAW "R650": DRAW "U250": DRAW "L650"
LOCATE 2, 1: PRINT "T"
LOCATE 3, 1: PRINT "E"
LOCATE 4, 1: PRINT "M"
LOCATE 5, 1: PRINT "P"
LOCATE 6, 1: PRINT "E"
LOCATE 7, 1: PRINT "R"
LOCATE 8, 1: PRINT "A"
LOCATE 9, 1: PRINT "T"
LOCATE 10, 1: PRINT "U"
LOCATE 11, 1: PRINT "R"
LOCATE 12, 1: PRINT "A"
LOCATE 14, 1: PRINT "(C)"
LOCATE 19, 6: PRINT "0": LOCATE 19, 73: PRINT INT(tmax)
LOCATE 20, 30: PRINT "VRIJEME (min) "
LOCATE 21, 10: PRINT "DINAMIKA SMRZAVANJA "; naziv$
PSET (45, 94): DRAW "R655"
LOCATE 1, 2: PRINT "+30"
LOCATE 7, 4: PRINT "0"
LOCATE 18, 2: PRINT "-50"

REM " podaci
PSET (50, yos(1))
FOR i = 2 TO 63
LINE (50 + tos(i - 1), yos(i - 1))-(50 + tos(i), yos(i))
NEXT