PDE Toolboxで熱解析

www.youtube.com

clear; close all; clc;

pdem = createpde('thermal');
gm = multicylinder([20 25 35], 20, 'Void', [1 0 0]);
pdem.Geometry  = gm;
generateMesh(pdem);
thermalProperties(pdem, 'cell', 1, 'ThermalConductivity', 40);
thermalProperties(pdem, 'cell', 2, 'ThermalConductivity', 0.15);
thermalBC(pdem, 'Face', 3, 'Temperature', 85);
thermalBC(pdem, 'Face', 7, 'Temperature', 4);
result = solve(pdem);

subplot(131);
pdegplot(gm);
subplot(132);
pdeplot3D(pdem.Mesh);
subplot(133);
pdeplot3D(pdem, 'ColormapData', result.Temperature);
sgtitle('Heat transfer analysis using PDE toolbox');