Статья:

Моделирование процесса температуропроводности в тонкой прямоугольной пластине

Конференция: L Студенческая международная научно-практическая конференция «Молодежный научный форум: технические и математические науки»

Секция: Физико-математические науки

Выходные данные
Сахибназарова В.Б. Моделирование процесса температуропроводности в тонкой прямоугольной пластине // Молодежный научный форум: Технические и математические науки: электр. сб. ст. по мат. L междунар. студ. науч.-практ. конф. № 10(50). URL: https://nauchforum.ru/archive/MNF_tech/10(50).pdf (дата обращения: 10.12.2024)
Лауреаты определены. Конференция завершена
Эта статья набрала 0 голосов
Мне нравится
Дипломы
лауреатов
Сертификаты
участников
Дипломы
лауреатов
Сертификаты
участников
на печатьскачать .pdfподелиться

Моделирование процесса температуропроводности в тонкой прямоугольной пластине

Сахибназарова Виктория Бахтиёровна
студент, СНИУ им. академика С.П. Королева, РФ, г. Самара

 

В данной работе в качестве объекта исследования рассматривается квадратная пластина со стороной равной 1. Полагается также, что процесс температуропроводности описывается однородным уравнением с коэффициентом a=1:

u't - ( u"xx + u"yy ) = 0.

Начальная температура пластины: u(x,y,0)=0, т.е. температура внутри тела (за исключением границы) везде равна 0. Длина каждого участка изоляции, заданной на границе, равна 0,5.

Граничные условия:

·     Участок A: .

·     Участок B: теплоизоляция.

·     Участок C: .

·     Участок D: .

·     Участок E: теплоизоляция.

·     Участок F: .

Положим, что процесс температуропроводности внутри пластины описывается однородным дифференциальным уравнением, с коэффициентом a=1[1]:

u't - ( u"xx + u"yy ) = 0.

Положим также, что температура вдоль границ пластины B, E со временем не меняется (статична), а значения температуры вдоль остальных границ образуются граничными условиями 1 рода:

·     Сторона A: .

·     Сторона F:

·     Сторона D:

·     Сторона C:  

Пусть на границе B присутствует термоизоляция. Угол α отклонения нормали  к границе B от оси x равен π. Тогда угол β отклонения нормали от оси y будет равен  - π = . Тогда граничное условие термоизоляция (условие 2-го рода) будет иметь вид:

·     Сторона B: .

·     Сторона E:

В качестве начальных условий положим, что при t=0 температура внутри пластины во всех точках равна 0:

·     ,

т.е. температура внутри тела (за исключением границы) везде равна 0.

Приведённое описание предмета моделирования не даёт оснований для введения условий 3-го рода.

Для исследования динамики температуропроводности внутри пластины найдём решение дифференциального уравнения методом конечных разностей. Для этого заменим все дифференциальные уравнения конечно-разностными аналогами.

Поскольку пластина квадратная, то шаги дискретизации выберем одинаковые по x и y, т.е. hx=hy=h , следовательно l=hx/hy=1. Так как сторона квадрата равна 1, то h=1/n, n – количество делений дискретизации. Тогда шаг дискретизации по времени будет равен t = mh2.

Конечно-разностная форма однородного дифференциального уравнения в данном случае примет вид:

.

Конечно-разностное уравнение граничного условия термоизоляции на границах B и E получим из полученного выше соотношения:

,

,

откуда следует, что на границах B и E соответственно:

, "k=0,1,2,…

, "k=0,1,2,…

Построим теперь начальную матрицу температур тела  для n=10, используя ранее полученные

уравнения граничных и начальных условий. Вид матрицы U0 приведен на рисунке 1.

 

Рисунок 1. Результаты вычислений в графической форме

 

«Внутренние» элементы матрицы отражают начальное условие и будут равны 0. Первая строка матрицы будет соответствовать температуре на границе A и условию термоизоляции на границе B. Температура на границе B соответствует температуре в ближайших внутренних точках, т.е. равна 0. Первый столбец матрицы будет соответствовать температуре на границе C. Последняя строка матрицы будет соответствовать температуре на границе D и условию термоизоляции на границе E, а последний столбец – на границе F.

Полученная матрица определяет температуру в точках пластины в нулевой момент времени (начальный временной слой). Следующие временные слои вычисляются уже с использованием полученных выше зависимостей.

 

Рисунок 2. Результаты вычислений в матричной форме

 

Рисунок 3. Результаты вычислений в графической форме

 

Исходя из проведенных исследований, можно сделать вывод, что на точность прогноза влияет не только предельное эвклидово расстояние между входным вектором и центром кластера и размер обучающей выборки, но также и размер плавающего окна. Наилучшие результаты по прогнозу дают обучающая выборка размером 240, предельное евклидово расстояние равное 0,005 и размер плавающего окна, равный 3.

 

Список литературы:

1.  Баландин А. В. Моделирование процесса температуропроводности в тонкой прямоугольной пластине. Методические указания. [Текст]. – Самара: СГАУ, 2014. – 27 с.