Найдите тепловые граничные условия, присвоенные геометрической области
tbca = findThermalBC(thermalmodel.BoundaryConditions,RegionType,RegionID)
Создайте тепловую модель и включите квадратную геометрию.
thermalmodel = createpde('thermal'); geometryFromEdges(thermalmodel,@squareg); pdegplot(thermalmodel,'EdgeLabels','on') ylim([-1.1 1.1]) axis equal
Примените граничные условия температуры к ребрам 1 и 3 квадрата.
thermalBC(thermalmodel,'Edge',[1 3],'Temperature',100);
Примените граничное условие теплового потока к ребру 4 квадрата.
thermalBC(thermalmodel,'Edge',4,'HeatFlux',20);
Проверьте спецификацию граничных условий на ребре 1.
tbcaEdge1 = findThermalBC(thermalmodel.BoundaryConditions,'Edge',1)tbcaEdge1 =
ThermalBC with properties:
RegionType: 'Edge'
RegionID: [1 3]
Temperature: 100
HeatFlux: []
ConvectionCoefficient: []
Emissivity: []
AmbientTemperature: []
Vectorized: 'off'
Проверьте спецификации граничных условий на ребрах 3 и 4.
tbca = findThermalBC(thermalmodel.BoundaryConditions,'Edge',[3:4]);
tbcaEdge3 = tbca(1)tbcaEdge3 =
ThermalBC with properties:
RegionType: 'Edge'
RegionID: [1 3]
Temperature: 100
HeatFlux: []
ConvectionCoefficient: []
Emissivity: []
AmbientTemperature: []
Vectorized: 'off'
tbcaEdge4 = tbca(2)
tbcaEdge4 =
ThermalBC with properties:
RegionType: 'Edge'
RegionID: 4
Temperature: []
HeatFlux: 20
ConvectionCoefficient: []
Emissivity: []
AmbientTemperature: []
Vectorized: 'off'
Создайте тепловую модель и включите блочную геометрию.
thermalmodel = createpde('thermal','transient'); gm = importGeometry(thermalmodel,'Block.stl'); pdegplot(thermalmodel,'FaceLabels','on','FaceAlpha',0.5)
Примените граничное условие температуры к граням 1 и 3 блока.
thermalBC(thermalmodel,'Face',1,'Temperature',100); thermalBC(thermalmodel,'Face',3,'Temperature',300);
Примените граничное условие конвекции к граням 5 и 6 блока.
thermalBC(thermalmodel,'Face',[5,6],... 'ConvectionCoefficient',5,... 'AmbientTemperature',27);
Проверьте спецификацию граничных условий на гранях 1 и 3.
tbca = findThermalBC(thermalmodel.BoundaryConditions,'Face',[1,3]);
tbcaFace1 = tbca(1)tbcaFace1 =
ThermalBC with properties:
RegionType: 'Face'
RegionID: 1
Temperature: 100
HeatFlux: []
ConvectionCoefficient: []
Emissivity: []
AmbientTemperature: []
Vectorized: 'off'
tbcaFace3 = tbca(2)
tbcaFace3 =
ThermalBC with properties:
RegionType: 'Face'
RegionID: 3
Temperature: 300
HeatFlux: []
ConvectionCoefficient: []
Emissivity: []
AmbientTemperature: []
Vectorized: 'off'
Проверьте спецификации граничных условий на гранях 5 и 6.
tbcaFace5 = findThermalBC(thermalmodel.BoundaryConditions,'Face',5)tbcaFace5 =
ThermalBC with properties:
RegionType: 'Face'
RegionID: [5 6]
Temperature: []
HeatFlux: []
ConvectionCoefficient: 5
Emissivity: []
AmbientTemperature: 27
Vectorized: 'off'
tbcaFace6 = findThermalBC(thermalmodel.BoundaryConditions,'Face',6)tbcaFace6 =
ThermalBC with properties:
RegionType: 'Face'
RegionID: [5 6]
Temperature: []
HeatFlux: []
ConvectionCoefficient: 5
Emissivity: []
AmbientTemperature: 27
Vectorized: 'off'