Создайте тепловую модель, импортируйте геометрию и установите начальную температуру равной 0 на всей геометрии.

thermalModel = createpde('thermal','transient');
geometryFromEdges(thermalModel,@lshapeg);
thermalIC(thermalModel,0)

ans = 
  GeometricThermalICs with properties:

            RegionType: 'face'
              RegionID: [1 2 3]
    InitialTemperature: 0

Установите различные начальные условия на каждом фрагменте L-образной мембранной геометрии.

Создайте модель и включите 2-D геометрию.

thermalModel = createpde('thermal','transient');
geometryFromEdges(thermalModel,@lshapeg);
pdegplot(thermalModel,'FaceLabels','on')
axis equal
ylim([-1.1 1.1])

Установите начальные условия.

thermalIC(thermalModel,0,'Face',1)

ans = 
  GeometricThermalICs with properties:

            RegionType: 'face'
              RegionID: 1
    InitialTemperature: 0

thermalIC(thermalModel,10,'Face',2)

ans = 
  GeometricThermalICs with properties:

            RegionType: 'face'
              RegionID: 2
    InitialTemperature: 10

thermalIC(thermalModel,75,'Face',3)

ans = 
  GeometricThermalICs with properties:

            RegionType: 'face'
              RegionID: 3
    InitialTemperature: 75

Используйте указатель на функцию, чтобы задать начальную температуру, которая зависит от координат.

Создайте тепловую модель для переходного анализа и включите геометрию. Геометрия представляет собой стержень с круглым сечением. Модель 2-D является прямоугольной полосой, размер y которой простирается от оси симметрии до внешней поверхности и размер X которой простирается по фактической длине стержня.

thermalmodel = createpde('thermal','transient');
g = decsg([3 4 -1.5 1.5 1.5 -1.5 0 0 .2 .2]');
geometryFromEdges(thermalmodel,g);

Установите начальную температуру в штоке в зависимость от координат y, например, ${10}^3\left(0.2-{\mathit{y}}^2\right)$.

T0 = @(location)10^3*(0.2 - location.y.^2);
thermalIC(thermalmodel,T0)

ans = 
  GeometricThermalICs with properties:

            RegionType: 'face'
              RegionID: 1
    InitialTemperature: @(location)10^3*(0.2-location.y.^2)

Создайте тепловую модель и включите квадратную геометрию.

thermalmodel = createpde('thermal','transient');
geometryFromEdges(thermalmodel,@squareg);
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]

pdeplot(thermalmodel,'XYData',result1.Temperature(:,end))

Теперь возобновите анализ и решите задачу за время от 10 до 15 секунд. Используйте полученный ранее решение в течение 10 секунд в качестве начального условия. Поскольку 10 секунд является последним элементом в tlist, вам не нужно указывать временной индекс решения. По умолчанию thermalIC использует последний индекс решения.

thermalIC(thermalmodel,result1)

ans = 
  NodalThermalICs with properties:

    InitialTemperature: [1541x1 double]

Решите задачу и постройте график решения.

result2 = solve(thermalmodel,10:0.5:15)

result2 = 
  TransientThermalResults with properties:

      Temperature: [1541x11 double]
    SolutionTimes: [10 10.5000 11 11.5000 12 12.5000 13 13.5000 14 14.5000 15]
       XGradients: [1541x11 double]
       YGradients: [1541x11 double]
       ZGradients: []
             Mesh: [1x1 FEMesh]

pdeplot(thermalmodel,'XYData',result2.Temperature(:,end))

Чтобы использовать ранее полученное решение для определенного времени решения вместо последнего, задайте индекс времени решения как третий параметр thermalIC. Например, используйте решение в момент времени 5 секунд, которое является 11-м элементом в tlist.

tlist(11)

ans = 5

thermalIC(thermalmodel,result1,11);

result2 = solve(thermalmodel,5:0.5:15)

result2 = 
  TransientThermalResults with properties:

      Temperature: [1541x21 double]
    SolutionTimes: [1x21 double]
       XGradients: [1541x21 double]
       YGradients: [1541x21 double]
       ZGradients: []
             Mesh: [1x1 FEMesh]

pdeplot(thermalmodel,'XYData',result2.Temperature(:,end))