问圆柱形翅片(销鳍)的Simscape模型

EN

% Variables
r       = 0.01;             % Radius of the pin fin in m
l       = 0.2;              % Length of the pin fin in m
h_conv  = 500;              % Convective heat transfer in W/m²K
h_cond  = 500;              % Conductive heat transfer in W/mK
T_base  = 140;              % Temperature on the left end of the fin tip in °C
T_amb   = 40;               % Ambient temperature in °C    
n_elem  = 20;               % Number of division
x       = [0:l/n_elem:l]';  % Vector of division, necessary for evaluation

A       = r^2 * pi;         % Cross sectional area of the pin fin in m²
S       = 2 * pi * r * l;   % Surface area of the pin in m²

% Analytical solution
beta    = sqrt((h_conv*S)/(h_cond*A));
f       = cosh(beta*(l-x))/cosh(beta*l);

temp    = f*(T_base-T_amb)+T_amb;

% Plot of the temperature distribution
plot (x,temp);
axis([0 0.2 40 140]);

x_ss     = [0:0.05:0.2]';                    % Vector of division, necessary for evaluation of the Simscape results
temp_ss  = [T_base,temp_simscape(end,:)]';   % Steady state results of Simscape model at 1/4, 2/4, 3/4 and 4/4 of the length

% Plot analytical vs. Simscape solution
plot (x,temp);
hold on;
plot (x_ss,temp_ss,'-o');
axis([0 0.2 40 140]);