Поиск свойств теплового материала, назначенных геометрической области
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 задайте следующие тепловые свойства:
Теплопроводность 10 )
Массовая плотность 1
Удельная теплота составляет )
thermalProperties(thermalmodel,'ThermalConductivity',10,... 'MassDensity',1,... 'SpecificHeat',0.1,... 'Face',1);
Для грани 2 задайте следующие тепловые свойства:
Теплопроводность 20 )
Массовая плотность 2
Удельная теплота 0,2 )
thermalProperties(thermalmodel,'ThermalConductivity',20,... 'MassDensity',2,... 'SpecificHeat',0.2,... 'Face',2);
Для грани 3 задайте следующие тепловые свойства:
Теплопроводность 30 ))
Массовая плотность 3
Удельная теплота 0,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
Vertices: [4x3 double]
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: []