Найдите тепловые свойства материала присвоенными геометрической области
tmpa = findThermalProperties(thermalmodel.MaterialProperties,RegionType,RegionID)tmpa присвоенный заданной области.
Создайте переходную тепловую модель, которая имеет три поверхности.
thermalmodel = createpde('thermal','transient'); geometryFromEdges(thermalmodel,@lshapeg); pdegplot(thermalmodel,'FaceLabels','on') ylim([-1.1,1.1]) axis equal
Для поверхности 1, задайте следующие тепловые свойства:
Теплопроводность
Массовая плотность
Удельная теплоемкость
thermalProperties(thermalmodel,'ThermalConductivity',10,... 'MassDensity',1,... 'SpecificHeat',0.1,... 'Face',1);
Для поверхности 2, задайте следующие тепловые свойства:
Теплопроводность
Массовая плотность
Удельная теплоемкость
thermalProperties(thermalmodel,'ThermalConductivity',20,... 'MassDensity',2,... 'SpecificHeat',0.2,... 'Face',2);
Для поверхности 3, задайте следующие тепловые свойства:
Теплопроводность )
Массовая плотность
Удельная теплоемкость
thermalProperties(thermalmodel,'ThermalConductivity',30,... 'MassDensity',3,... 'SpecificHeat',0.3,... 'Face',3);
Проверяйте спецификацию свойств материала на поверхность 1.
mpaFace1 = findThermalProperties(thermalmodel.MaterialProperties,'Face',1)mpaFace1 =
ThermalMaterialAssignment with properties:
RegionType: 'face'
RegionID: 1
ThermalConductivity: 10
MassDensity: 1
SpecificHeat: 0.1000
Проверяйте спецификацию источника тепла на поверхности 2 и 3.
mpa = findThermalProperties(thermalmodel.MaterialProperties,'Face',[2,3]);
mpaFace2 = mpa(1)mpaFace2 =
ThermalMaterialAssignment with properties:
RegionType: 'face'
RegionID: 2
ThermalConductivity: 20
MassDensity: 2
SpecificHeat: 0.2000
mpaFace3 = mpa(2)
mpaFace3 =
ThermalMaterialAssignment with properties:
RegionType: 'face'
RegionID: 3
ThermalConductivity: 30
MassDensity: 3
SpecificHeat: 0.3000
Создайте геометрию, которая состоит из трех сложенных цилиндров, и включайте геометрию в тепловую модель.
gm = multicylinder(10,[1 2 3],'ZOffset',[0 1 3])gm =
DiscreteGeometry with properties:
NumCells: 3
NumFaces: 7
NumEdges: 4
NumVertices: 4
thermalmodel = createpde('thermal'); thermalmodel.Geometry = gm; pdegplot(thermalmodel,'CellLabels','on','FaceAlpha',0.5)
Теплопроводность цилиндра C1 .
thermalProperties(thermalmodel,'ThermalConductivity',10,'Cell',1);
Теплопроводность цилиндра C2 .
thermalProperties(thermalmodel,'ThermalConductivity',20,'Cell',2);
Теплопроводность цилиндра C3 .
thermalProperties(thermalmodel,'ThermalConductivity',30,'Cell',3);
Проверяйте спецификацию свойств материала на ячейку 1:
mpaCell1 = findThermalProperties(thermalmodel.MaterialProperties,'Cell',1)mpaCell1 =
ThermalMaterialAssignment with properties:
RegionType: 'cell'
RegionID: 1
ThermalConductivity: 10
MassDensity: []
SpecificHeat: []
Проверяйте спецификацию источника тепла на ячейки 2 и 3:
mpa = findThermalProperties(thermalmodel.MaterialProperties,'Cell',2:3);
mpaCell2 = mpa(1)mpaCell2 =
ThermalMaterialAssignment with properties:
RegionType: 'cell'
RegionID: 2
ThermalConductivity: 20
MassDensity: []
SpecificHeat: []
mpaCell3 = mpa(2)
mpaCell3 =
ThermalMaterialAssignment with properties:
RegionType: 'cell'
RegionID: 3
ThermalConductivity: 30
MassDensity: []
SpecificHeat: []