2018-02-14
295
Дифференциальное уравнение нестационарной теплопроводности в случае одномерной задачи имеет вид
. (1)
В лабораторной работе №2 был изложен способ решения дифференциального уравнения теплопроводности (1) методом конечных разностей по явной схеме. Как было указано, решение задач по явной схеме будет устойчивым, если выполняется условие
. (2)
Последнее условие весьма обременительно. Как показывают практические расчеты, при достаточно малых ∆ х и конкретных значениях а величина ∆ τ оказывается очень малой и приходится медленно продвигаться по времени, т.е. делать большое число шагов по времени τ. Все это чрезвычайно повышает трудоемкость решения.
Рассмотрим теперь принципиально другое сеточное уравнение. Если взять приближенное значение производной по времени «назад» (см.рис.1.),
Рис.1. Пространственно-временная область одномерной задачи теплопроводности
то получим следующее конечно-разностное соотношение:
. (3)
Уравнение (3) решается труднее, поскольку в него входят три неизвестные температуры: 
Поэтому в данном случае нужно решать сразу всю систему разностных уравнений типа (3) для всех точек i, k сетки. Здесь так же, как и при решении уравнения Лапласа методом сеток, можно применить метод итераций.
Неявные разностные уравнения решаются сложнее, чем явные, но они абсолютно устойчивы при любом шаге по времени. Это позволяет выбирать шаг ∆ τ значительно большим, чем в явных схемах, и соответственно уменьшать общее время счета всей задачи.
Для реализации численного метода решения задач теплопроводности на ЭВМ разработан специальный метод прогонки, более эффективный, чем метод итераций.
Изложим метод прогонки на примере предыдущей задачи, предложенной в лабораторной работе №2. Уравнение (3) можно, очевидно, записать еще и в виде
(4)
при начальных и граничных условиях:
в момент времени k =0
, (5)
и задана в точках: i =1,2,…, n -1;
(6)
. (7)
Уравнение (4) для удобства дальнейших выводов запишем в виде
, (8)
где 
.
Идея метода прогонки заключается в следующем: связь между двумя соседними узлами 
и 
представляется в виде
, (9)
где 
и 
некоторые коэффициенты, подлежащие определению. Если они будут известны, то делая «прогонку» в направлении справа-налево, начиная с правого граничного условия 
, можно по (9) последовательно найти все температуры в (k +1)-м слое по времени.
Из (9), заменяя i на i -1, имеем
. (10)
Подставляя (10) в формулу (8), получим
. (11)
Отсюда
. (12)
Сравнивая (12) с (9), получим следующие формулы:
(13)
При i =1 из (8) найдем
. (14)
Или, используя граничное условие (6), из (14) получим
. (15)
С другой стороны, из (9) найдем
. (16)
Сравнивая (16) и (15), получим
(17)
Пользуясь системами (А) и (В), производя «прогонку» слева-направо в прямом направлении, последовательно найдем все коэффициенты 
.
Затем, используя «обратный ход», т.е. прогонку справа-налево, начиная с 
, как уже было указано, по (9) найдем все температуры в (k +1) слое, если известно распределение температур в k -м слое. Таким образом, указан переход от k -го слоя по времени к (k +1) слою.
Следовательно, отправляясь от известного начального (нулевого) слоя, можно построить решение 
во всех точках сетки i, k.
Постановка задачи
Найти температурное поле в теплоизолированном с боковой поверхности ограниченном стержне, разделенным на 
частей, при следующих краевых условиях:
1) в начальный момент времени τ =0, t = t 0;
2) при х =0, t = tw 1= const;
3) при х = l, где l – длина стержня, t = tw 2= const.
Шаг по времени рассчитать по формуле: 
, где 
– шаг по координате.
Числовые значения исходных данных взять из таблицы 11.
Часть таблицы 11
| Номер варианта | а∙ 106, м2/с | l, м | t 0, oC | tw1, oC | tw2, oC | 
|
| 1 | 7,5 | 0,01 | 600 | 50 | 50 | 9 |
Содержание отчета
1. Представить расчет температурного поля для одного шага по времени с необходимыми пояснениями, применяя метод прогонки.
2. Представить расчет температурного поля для одного шага по времени с необходимыми пояснениями, применяя на выбор один из различных программных продуктов (Mahtcad, он-лайн калькуляторы системы Интернет или др.).
3. Представить график распределения температур в стержне в различные моменты времени (для 5 шагов по времени).
Таблица 13
Матрица соответствия оценочных средств запланированным результатам обучения по дисциплине (модулю)
«Применение электронно-вычислительных машин в тепловых расчетах»
| Компетенции | Оценочные средства | ||||||||
| Конспекты вопросов и тем лекционного курса (Раздел №1) | Конспекты вопросов и тем лекционного курса (Раздел №2) | Конспекты вопросов и тем лекционного курса (Раздел №3) | Конспекты вопросов и тем лекционного курса (Раздел №4) | Отчет по лабораторной работе № 1 (Раздел №2) | Отчет по лабораторной работе № 2 (Раздел №3) | Отчет по контрольной работе (Раздел №3) | Вопросы к зачету с оценкой (Раздел №1-4) | ||
|
Самостоятельная работа студента |
Лабораторные работы |
Контрольная работа | Промежуточная аттестация – зачет с оценкой | ||||||
| ПК-1 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 | В (ПК-1) – I2 У (ПК-1) – I2 З (ПК-1) – I2 |
