% Soojenemine soojuskandja temperatuuri mõjul:   c\cdot m\frac{\textrm{dT}\left(t\right)}{\textrm{dt}}=\textrm{UA}\left(T_i \left(t\right)-T\left(t\right)\right)
% T(t) on sisetemperatuur, Ti(t) on soojuskandja temperatuur, soojusvool on \frac{\mathrm{d}Q\left(t\right)}{\textrm{dt}} [W]

U = 1200 % soojusjuhtivus [W/(m2*K)]
A = 1 % soojusvahetuse pindaja [m2]
c = 4180; % erisoojus [J/(kg*K)]
m = 1000 % soojenev mass [kg]

%
% {\mathrm{U}}_k {\mathrm{A}}_k \left(T\left(t\right)-T_k \left(t\right)\right) on soojusvool keskkonda [W]
%

Ak = 6 % soojenevat massi ümbritseva kesta pindala [m2] h=2,D=0.8m
Uk = 12 % soojusjuhtivus keskkonda [W/(m2*K)] soojustusega 1.2
UAk = Uk*Ak; % soojuserikaod läbi kogu pinna A [W/K]

%
% Soojenemisprotsess kütmise vahendusel C\frac{\textrm{dT}\left(t\right)}{\textrm{dt}}=P\left(t\right)-{\mathrm{U}}_k A_k \left(T\left(t\right)-T_k \left(t\right)\right),
%  ja jahtumine väliskeskkonda; Tk(t) on väliskeskonna temperatuur, P(t) on küttevõimsus [W]

T1 = 0.97*3600 % soojuslik ajakonstant [s]
C = T1*U*A % soojuserimahtuvus [J/K]

% \frac{C}{{\mathrm{U}}_k A_k }\frac{\textrm{dT}\left(t\right)}{\textrm{dt}}=\frac{1}{{\mathrm{U}}_k A_k }P\left(t\right)-\left(T\left(t\right)-T_k \left(t\right)\right) ja
% I järku üldkuju: T_1 \frac{\textrm{dT}\left(t\right)}{\textrm{dt}}={\mathrm{K}}_o u\left(t\right)-\left(T\left(t\right)-T_k \left(t\right)\right)
% ehk ka \frac{\textrm{dT}\left(t\right)}{\textrm{dt}}=\frac{1}{C}P\left(t\right)-\frac{{\mathrm{U}}_k A_k }{C}\left(T\left(t\right)-T_k \left(t\right)\right) vt. simulatsiooniskeemi soojen211.slx
%
% Soojusmassilevi (aine hulk temperatuuriga Tin toob soojushulga Q [J])
% Q=c\cdot m\cdot T_{\textrm{in}}   ja  \frac{\mathrm{d}Q\left(t\right)}{\textrm{dt}}=c\frac{\textrm{dm}\;}{\textrm{dt}}\left(T_{\textrm{in}} -T\right)

