clc, clear

% Fizicki parametri valjka
P = 500; % Snaga zagrevanja [W]
L = 1; % Duzina valjka [m]
D = 0.1; % Precnik valjka [m]
m = 2; % Broj delova svakog materijala
d = L /(2 * m); %Duzina jednog delica [m]
alfa0 = 5; % Koeficijent prelaska toplote strujanjem [W/m2K]
alfa1 = 0.5;
tAlfa0 = 20; % Temperatura vazduha pri kojoj je odredjeno alfa
t0=20; %Temperatura vazduha
exp=0.8; %Eksponent za proracun alfa
Pgama_d = 100 / 2 / m; %Snaga generisana u jednom delicu [W]

% Materijal (Bakar)
lambda = 401; % Specificna toplotna provodnost [W/mK]
cp = 385; % Specificna toplotna kapacitivnost [J/kgK]
ro = 8960; % Gustina [kg/m3]

Spr = pi*(D^2)/4; % Povrsina poprecnog preseka valjka
Som = D*pi*d; % Povrsina omotaca jednog m-tog delica valjk

m = ro * d * Spr; % Masa jednog delica

Temp1_poc = t0; % Pocetna vrednost temperature (integratora) delica 1
Temp2_poc = t0; % Pocetna vrednost temperature (integratora) delica 2
Temp3_poc = t0; % Pocetna vrednost temperature (integratora) delica 3
Temp4_poc = t0; % Pocetna vrednost temperature (integratora) delica 4
