Решение задачи. 1. Тепловая схема однородного стержня .. ..3

Оглавление

1. Тепловая схема однородного стержня……………………………..……..3

2. Двумерная пластина с теплообменом с поверхности в среду…………9

3. Нестационарный процесс распределения температур в однородном стержневом элементе…………………………………………………………19

5. Список использованной литературы…………………………………...25

Тепловая схема однородного стержня

Рассмотрим одномерный стержень, поперечное сечение которого столь мало, чтобы можно было пренебречь изменением температуры по его сечению. В этом случае, температурное поле стержня является одномерным и изменяется только по оси х, направленной по длине стержня.

Постановка задачи. Зададим граничные условия 3 рода, т.е. на границах стержня происходит теплообмен со средой. Для решения задачи примем следующие начальные данные:

· температура среды и коэффициент теплоотдачи со стороны левой границы стержня равны: Та1=400 ⁰С и α 1=5 Вт/м^{2 0}С, со стороны правой – Та2=100⁰С и α 2=20 Вт/м^{2 0}С;

· длина стержня L равна 200 мм;

· теплопроводность стержня λ = 0, 1 Вт/м ⁰С;

· радиус стержневого элемента r= 5 мм.

Решение задачи.

Разобьем длину стержня сечениями перпендикулярными оси x на три конечных объема длиной Δ x=0.01мм и площадью сечения A, равной площади поперечного сечения стержня. Соответствующая тепловая схема приведена на рис.1.1. В центре каждого объема поместим по одному узлу, при этом номер узла совпадает с номером объема. Пронумеруем узлы тепловой схемы так, как показано на рис.1. от узла 1 (на левом торце стержня) до узла 3 (на правом торце стержня).

а)

б)

Рисунок 1. Одномерный стержень, разбитый на конечные объемы (а) и его тепловая схема (б)

Рассмотрим баланс потока теплоты в i-ом выделенном объеме, воспользовав­шись интегральным уравнением теплового баланса

(1)

где Vi = A ∆ xi – объем i-го элемента;

S_i – площадь всей поверхности выде­ленного i-го объема. Индекс i относит рассматриваемые переменные к i-му выделенному объему.

Поверхностный интеграл в левой части уравнения (1.1) выражает суммар­ный тепловой поток, пересекающий поверхность выделенного i-ro объема. Учи­тывая, что тепловой поток вдоль стержня одномерен (не изменяется в направле­нии перпендикулярном оси х), а тепловой поток с боковой поверхности стержня отсутствует (поскольку рассматривается теплоизолированный с боковой по­верхности стержень) можно записать, например, для узла 2, что

(2)

где J₂ и J₃ - тепловые потоки на левой и павой границах выделенного объема.

За положительное направление вектора теплового потока принято направление, соответствующее вытеканию из объема теплового по­тока.

Так как, рассматриваемая задача является стационарной, т.е. dT/dt =0 и внутренние источники теплоты отсутствуют, то уравнение (1) принимает следующий вид:

(3)

Уравнения теплового баланса, записанные для узлов тепловой схемы, имеют вид:

(4)

где R1, R2, R3, R4 – тепловые сопротивления выделенных объемов между узлами.

Баланс потоков теплоты, протекающих в ветвях соединенных с узлами 1, 2, 3 выражаются следующими уравнениями

(5)

Ориентированный граф тепловой схемы представлен на рис.2. Номера ветвей указаны в кружках.

Рисунок.2. Ориентированный граф тепловой схемы

В стационарном случае, вектор-столбец температур Т описывается матричным уравнением, при , т.е.

(6)

Введя вектор-столбец тепловых потоков ветвей J=||J1 J2 J3 J4||^Т, систему уравнений можно записать в матричном виде

AJ=0 (7)

Матрица A называется матрицей инциденции, для рассматриваемого случая имеет размерность 3*4 и равна:

(8)

Уравнение (7) является, по существу, дискретным уравнением баланса тепловых потоков в тепловой схеме. Вид матрицы A нетрудно понять. Строки матрицы соответствуют узлам графа и расположены согласно их номерам от 1 до 3 сверху вниз, а столбцы матрицы соответствуют ветвям графа, причем номер столбца равен номеру ветви в тепловой схеме.

Разности температур в ветвях графа можно представить в виде вектора- столбца Δ T:

(9)

Введя вектор столбец температур узлов графа

(10)

простым перемножением легко убедится, что вектор-столбец (9) можно записать в следующем матричном виде:

(11)

(12)

где - вектор-столбец известных температур в ветвях 1 и 4 тепловой схемы.

Сравнение матрицы инциденций A (8) и матрицы в соотношении (11) показывает, что последняя матрица является транспонированной по отношению к матрице А, т.е. равна А^Т, поэтому вектор-столбец Δ Т (9) в соотношении (11) можно записать в матричном виде через транспонированную матрицу инциденций А^Т, т.е.

Δ Т=А^ТТ (13)

Полученные матрично-топологические соотношения (7) и (9) устанавливают связь между тепловыми потоками в ветвях тепловой схемы и преобразование узловых температур в разности температур в ветвях. Матрица инциденций А отображает структуру тепловой схемы.

Матрица А, естественным образом, была получена из системы уравнений баланса тепловых потоков в узлах графа.

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

Выше было показано, что тепловой поток J_i в i-ой ветви равен J_i=g_iΔ T_i.

Тогда связь векторов-столбцов J и Δ Т может быть записана в следующем матричном виде:

(14)

где G – квадратная матрица проводимостей ветвей размерностью М*М, М – количество ветвей графа.

Матрица проводимостей G формируется следующим образом: если ветвь i представляет собой тепловую проводимость g_i (кондуктивную или конвективную), то элемент ii матрицы G равен g_i.

Конвективные тепловые проводимости с торцов стержня, которые входят в ветви с номерами 1 и 4 равны соответственно и .

Тепловые кондуктивные проводимости, которые входят в ветви с номерами 2 и 3 равны и

Строим матрицу проводимостей G:

(15)

Введем матрицу B, которая находится по формуле:

(16)

где А^Т – транспонированная матрица А.

(17)

Подставляя выражения (8), (12), (15) и (16) в уравнение (18)

(18)

находим искомые температуры в узлах стержневого элемента.

(19)

Отсюда температуры равны: Т1=373.3 ⁰С; Т2=240 ⁰С; Т3=106.7 ⁰С.

Распределение температуры в стержневом термодинамическом элементе представлены на рис.3

Рисунок.3. График распределения температур по длине одномерного стержня

2. Двумерная пластина с теплообменом с поверхности в среду Рассмотрим двумерную пластину, пренебрегем изменением температуры по ее толщине. В этом случае, температурное поле пластины является двумерным и изменяется только по оси х и у.

Зададим граничные условия 3 рода, т.е. на границах пластины происходит теплообмен со средой. Для решения задачи примем следующие начальные данные:

· температуры сред и коэффициенты теплоотдачи со сторон границ пластины соответственно равны: Та1=400 ⁰С, Та2=50 ⁰С, Та3=500 ⁰С, Та4=60 ⁰С, α 1=3000 Вт/м^{2 0}С, α 2=60 Вт/м^{2 0}С, α 3=100 Вт/м^{2 0}С, α 4=2000 Вт/м^{2 0}С;

· длина пластины L равна 80 мм;

· теплопроводности материалов λ 1 = 50 Вт/м ⁰С, λ 2 =30 Вт/м ⁰С,

λ 3 = 390 Вт/м ⁰С, λ 4 = 400 Вт/м ⁰С;

· А0=0.0004 м²;

а) Та4

У₄ У₇ У₁₀

У₁ 1 У₂ 2 У₅ 3 У₈

R1 R2

R3 R4 R5

λ 1 λ 3

б)

g1 g2 g5 g8

g3 g6 g9

g12 g14

g11 g16

g13 g15 g17

g18 g19 g21 g23

g20 g22 g24

Рисунок 4. Двумерная пластина (а) и ее тепловая схема (б)

Рассмотрим баланс потока теплоты в i-ом выделенном объеме, воспользовав­шись интегральным уравнением теплового баланса

(20)

где Vi = A ∆ xi – объем i-го элемента;

S_i- площадь всей поверхности выде­ленного i -го объема. Индекс i относит рассматриваемые переменные к i-му выделенному объему.

Поверхностный интеграл в левой части уравнения (3.3.1) выражает суммар­ный тепловой поток, пересекающий поверхность выделенного i-ro объема. Так как, рассматриваемая задача является стационарной, т.е. dT/dt =0

Уравнения теплового баланса, записанные для узлов тепловой схемы, имеют вид:

где R – тепловые сопротивления выделенных объемов между узлами и рассчитываются и соответствуют кондуктивной проводимости

А на границе материалов тепловые сопротивления рассчитываются как , где λ i и λ g-коэффициенты теплопроводности граничащих материалов.

Конвективные тепловые проводимости, которые входят в ветви .

Соответственно тепловые проводимости равны

g1=1.2; g2=0.5; g3=0.5; g4=0.024; g5=0.3; g6=0.4; g7=0.024; g8=0.04; g9=0.3; g10=0.024; g11=1.2; g12=2.2; g13=3.9; g14=1.15; g15=2.95; g16=0.04; g17=2; g18=1.2; g19=3.9; g20=0.8; g21=2; g22=0.8; g23=0.04; g24=0.8

В стационарном случае, вектор-столбец температур Т описывается матричным уравнением, при , т.е.

(21)

Введя вектор-столбец тепловых потоков ветвей J=||J1 J2 J3 J4 J5 …… J24||^Т, систему уравнений можно записать в матричном виде

AJ=0 (22)

Матрица A называется матрицей инциденций, для рассматриваемого случая имеет размерность 9*24 и равна:

Она получается путем сложения по горизонтали отдельных матриц в программе MathCAD с помощью функции augment

А^Т – транспонированная матрица А

Уравнение (22) является, по существу, дискретным уравнением баланса тепловых потоков в тепловой схеме. Вид матрицы A нетрудно понять. Строки матрицы соответствуют узлам графа, а столбцы матрицы соответствуют ветвям графа, причем номер столбца равен номеру ветви в тепловой схеме.

Разности температур в ветвях графа можно представить в виде вектора- столбца Δ T:

.

Введя вектор столбец температур узлов графа

(23)

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

(24)

(25)

где Та - вектор-столбец известных температур в ветвях. Сравнение матрицы инциденций A и матрицы в соотношении (23) показывает, что последняя матрица является транспонированной по отношению к матрице А, т.е. равна

Δ Т=А^ТТ (26)

Полученные матрично-топологические соотношения (3.3.4) и (3.3.8) устанавливают связь между тепловыми потоками в ветвях тепловой схемы и преобразование узловых температур в разности температур в ветвях. Матрица инциденций А отображает структуру тепловой схемы.

Матрица А, естественным образом, была получена из системы уравнений баланса тепловых потоков в узлах графа.

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

Выше было показано, что тепловой поток J_i в i-ой ветви равен J_i=g_iΔ T_i.

Тогда связь векторов-столбцов J и Δ Т может быть записана в следующем матричном виде:

(27)

где G – квадратная матрица проводимостей ветвей размерностью М*М, М – количество ветвей графа.

Матрица проводимостей G формируется следующим образом: если ветвь i представляет собой тепловую проводимость g_i (кондуктивную или конвективную), то элемент ii матрицы G равен g_i.

Строим матрицу проводимостей G:

Она получается путем сложения по горизонтали отдельных матриц в программе MathCAD с помощью функции augment и сложения по вертикали с помощью функции stack.

Строим новую матричную функцию:

(28)

-вектор –столбец начальных температур в момент времени равный нулю в узлах тепловой схемы.

(29)

Подставляя в уравнение (29) находим искомые температуры в узлах стержневого элемента

Т1=402.253⁰С

Т2=399.665⁰С

Т3=379.652⁰С

Т4=426.676⁰С

Т5=431.819⁰С

Т6=429.166⁰С

Т7=435.115⁰С

Т8=441.048⁰С

Т9=442.65⁰С

Рисунок 5. График распределения температур полученный в программе MathCAD