Найдите тепловые граничные условия присвоенными геометрической области
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'