I think you should use the Simscape language and customize the equations of the Thermal Mass block to get what you want. The build process to get the custom block is at this link.
Just double-click the Thermal Mass block, click the "View Source" link and copy the code over to a new file. I got this working as shown in the attached image. The code is here:
component customMass
% Custom Thermal Mass
nodes
A = foundation.thermal.thermal; % :top
B = foundation.thermal.thermal; % :bottom
end
parameters
mass = { 1, 'kg' }; % Mass
sp_heat = { 447, 'J/(kg*K)' }; % Specific heat
end
variables
% Differential variables
T = {0, 'K'}; % Temperature difference between A and B
Q = { 0, 'J/s' }; % Heat flow
end
function setup
% Parameter range checking
if mass <= 0
pm_error('simscape:GreaterThanZero','Mass')
end
if sp_heat <= 0
pm_error('simscape:GreaterThanZero','Specific heat')
end
end
branches
Q : A.Q -> B.Q;
end
equations
T == A.T - B.T;
Q == mass * sp_heat * T.der;
end
end
- Sebastian
