Найдите тепловые начальные условия присвоенными геометрической области
tica = findThermalIC(thermalmodel.InitialConditions,RegionType,RegionID)
Создайте переходную тепловую модель, которая имеет три поверхности.
thermalmodel = createpde('thermal','transient'); geometryFromEdges(thermalmodel,@lshapeg); pdegplot(thermalmodel,'FaceLabels','on') ylim([-1.1 1.1]) axis equal
Установите начальные температуры для каждой поверхности.
thermalIC(thermalmodel,10,'Face',1); thermalIC(thermalmodel,20,'Face',2); thermalIC(thermalmodel,30,'Face',3);
Проверяйте первоначальную спецификацию условия на поверхность 1.
ticaFace1 = findThermalIC(thermalmodel.InitialConditions,'Face',1)ticaFace1 =
GeometricThermalICs with properties:
RegionType: 'face'
RegionID: 1
InitialTemperature: 10
Проверяйте первоначальные температурные технические требования на поверхности 2 и 3.
tica = findThermalIC(thermalmodel.InitialConditions,'Face',[2 3]);
ticaFace2 = tica(1)ticaFace2 =
GeometricThermalICs with properties:
RegionType: 'face'
RegionID: 2
InitialTemperature: 20
ticaFace3 = tica(2)
ticaFace3 =
GeometricThermalICs with properties:
RegionType: 'face'
RegionID: 3
InitialTemperature: 30
Создайте геометрию, которая состоит из трех вложенных цилиндров, и включайте геометрию в переходную тепловую модель.
gm = multicylinder([5 10 15],2)
gm =
DiscreteGeometry with properties:
NumCells: 3
NumFaces: 9
NumEdges: 6
NumVertices: 6
thermalmodel = createpde('thermal','transient'); thermalmodel.Geometry = gm; pdegplot(thermalmodel,'CellLabels','on','FaceAlpha',0.5)
Установите начальные температуры для каждой ячейки.
thermalIC(thermalmodel,0,'Cell',1); thermalIC(thermalmodel,100,'Cell',2); thermalIC(thermalmodel,0,'Cell',3);
Проверяйте первоначальную спецификацию условия на ячейку 1.
ticaCell1 = findThermalIC(thermalmodel.InitialConditions,'Cell',1)ticaCell1 =
GeometricThermalICs with properties:
RegionType: 'cell'
RegionID: 1
InitialTemperature: 0
Проверяйте первоначальную спецификацию условия на ячейки 2 и 3.
tica = findThermalIC(thermalmodel.InitialConditions,'Cell',[2:3]);
ticaCell2 = tica(1)ticaCell2 =
GeometricThermalICs with properties:
RegionType: 'cell'
RegionID: 2
InitialTemperature: 100
ticaCell3 = tica(2)
ticaCell3 =
GeometricThermalICs with properties:
RegionType: 'cell'
RegionID: 3
InitialTemperature: 0
Создайте тепловую модель и включайте квадратную геометрию.
thermalmodel = createpde('thermal','transient'); gm = @squareg; geometryFromEdges(thermalmodel,gm); pdegplot(thermalmodel,'FaceLabels','on') ylim([-1.1 1.1]) axis equal
Задайте свойства материала, источник тепла, установите начальные и граничные условия.
thermalProperties(thermalmodel,'ThermalConductivity',500,... 'MassDensity',200,... 'SpecificHeat',100); internalHeatSource(thermalmodel,2); thermalBC(thermalmodel,'Edge',[1 3],'Temperature',100); thermalIC(thermalmodel,0);
Сгенерируйте mesh и решите задачу.
generateMesh(thermalmodel); tlist = 0:0.5:10; result1 = solve(thermalmodel,tlist)
result1 =
TransientThermalResults with properties:
Temperature: [1541x21 double]
SolutionTimes: [1x21 double]
XGradients: [1541x21 double]
YGradients: [1541x21 double]
ZGradients: []
Mesh: [1x1 FEMesh]
Проверяйте в настоящее время активную первоначальную температурную спецификацию.
tica = findThermalIC(thermalmodel.InitialConditions,'Face',1)tica =
GeometricThermalICs with properties:
RegionType: 'face'
RegionID: 1
InitialTemperature: 0
Теперь возобновите анализ и решите задачу в течение многих времен с 10 до 15 секунд. Используйте ранее полученное решение в течение 10 секунд как начальное условие. С 10 секунд последний элемент в tlist, вы не должны задавать индекс времени решения. По умолчанию, thermalIC использует последний индекс решения.
ic = thermalIC(thermalmodel,result1);
Решите задачу
tlist = 10:0.5:15; result2 = solve(thermalmodel,tlist);
Проверяйте в настоящее время активную первоначальную температурную спецификацию.
tica = findThermalIC(thermalmodel.InitialConditions,'Face',1)tica =
NodalThermalICs with properties:
InitialTemperature: [1541x1 double]
pdeplot(thermalmodel,'XYData',tica.InitialTemperature)
Свойства GeometricThermalICs | Свойства NodalThermalICs | thermalIC