clc, clear

% Fizicki parametri valjka

P = 1000/2; % Snaga zagrevanja [W]
L = 1; % Duzina valjka [m]
D = 0.1; % Precnik valjka [m]
m = 2; % Broj delova svakog materijala
d= L/(2*m); % Debljina jednog delica valjka

alfa = 20; % Koeficijent prelaska toplote strujanjem [W/m2K]
t0 = 22; % Temperatura vazduha 

% Materijal 1 (Bakar)
lambda1 = 401; % Specificna toplotna provodnost [W/mK]
cp1 = 385; % Specificna toplotna kapacitivnost [J/kgK]
ro1 = 8960; % Gustina [kg/m3]

% % Materijal 1 (Aluminijum)
% lambda1 = 237; % Specificna toplotna provodnost [W/mK]
% cp1 = 897; % Specificna toplotna kapacitivnost [J/kgK]
% ro1 = 2700; % Gustina [kg/m3]

% % Materijal 2 (Bakar)
% lambda2 = 401; % Specificna toplotna provodnost [W/mK]
% cp2 = 385; % Specificna toplotna kapacitivnost [J/kgK]
% ro2 = 8960; % Gustina [kg/m3]

% Materijal 2 (Aluminijum)
lambda2 = 237; % Specificna toplotna provodnost [W/mK]
cp2 = 897; % Specificna toplotna kapacitivnost [J/kgK]
ro2 = 2700; % Gustina [kg/m3]

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

RTpr1 = (1/lambda1)*d/Spr; % Toplotni otpor provodjenju za materijal 1
RTpr2 = (1/lambda2)*d/Spr; % Toplotni otpor provodjenju za materijal 2

RTstr_om = 1 / (alfa * Som); % Toplotni otpor strujanju sa omotaca valjka
RTstr_b = 1 / (alfa * Spr); % Toplotni otpor strujanju sa gornjeg bazisa valjka

m1 = ro1 * d * Spr; % Masa m-tog delica valjka materijala 1
m2 = ro1 * d * Spr; % Masa m-tog delica valjka materijala 2

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

set_param('lab6Termal2014b', 'SimulationCommand', 'start'); % Pokretanje simulacije
pause(2); 

%____________________________________________________________________________________________%

% Deo koda koji omogucava automatsko pokretanje drugog dela proracuna 
% za uvazavanje promene snage zagrevanja 

Temperature0=Temperature; % Cuvanje temperatura iz prvog dela ogleda (Matrica iz Scope-a Temperature)

Temp1_poc = t0 + Temperature0(length(Temperature0(:,3)),3); % Alternativno: Temperature1(length(Temperature1),3)
Temp2_poc = t0 + Temperature0(length(Temperature0(:,4)),4); % Alternativno: Temperature1(length(Temperature1),4)
Temp3_poc = t0 + Temperature0(length(Temperature0(:,5)),5); % Alternativno: Temperature1(length(Temperature1),5)
Temp4_poc = t0 + Temperature0(length(Temperature0(:,6)),6); % Alternativno: Temperature1(length(Temperature1),6)

P = 2*P; % Nova snaga zagrevanja [W]

set_param('lab6Termal2014b', 'SimulationCommand', 'start'); % Pokretanje simulacije
pause(2);

Temperature(:,1)=Temperature(:,1)+Temperature0(length(Temperature0),1); % Pomeranje vremenske ose u cilju nastavljanja na prvu simulaciju

TemperatureUk = [Temperature0; Temperature]; % Kreiranje matrice sa podacima iz obe simulacije

figure(2) % Iscrtavanje temperatura svakog dela valjka
plot(TemperatureUk(:,1),TemperatureUk(:,2:7))
title('Promena porasta temperature pojedinih delova valjka')
xlabel('Vreme (s)')
ylabel('Porast temperature (K)') 
legend('Osnova valjka', 'Prvi deo valjka', 'Drugi deo valjka', 'Treci deo valjka', 'Cetvrti deo valjka', 'Gornji bazis valjka')


