Вывод дифференциального уравнения теплопроводности при нагреве стержня

Лекция 4. Вывод уравнения теплопроводности

При построении математической модели распространения тепла в стержне сделаем следующие предположения:
1) стержень сделан из однородного проводящего материала с плотностью ρ;
2) боковая поверхность стержня теплоизолирована, то есть тепло может распространяться только вдоль оси ОХ;
3) стержень тонкий — это значит, что температура во всех точках любого поперечного сечения стержня одна и та же.

Рассмотрим часть стержня на отрезке [х, х + ∆х] (см. рис. 6) и воспользуемся законом сохранения количества тепла:

Общее количество тепла на отрезке [х, х + ∆х] = полному количеству тепла, прошедшему через границы + полное количество тепла, образованного внутренними источниками.

Общее количество тепла, которое необходимо сообщить участку стержня, чтобы повысить его температуру на ∆U, вычисляется по формуле: ∆Q= CρS∆x∆U, где С — удельная теплоемкость материала ( = количеству тепла, которое нужно сообщить 1 кг вещества, чтобы поднять его температуру на 1°), S — площадь поперечного сечения.

Количество тепла, прошедшее через левый конец участка стержня за время ∆t (тепловой поток) вычисляется по формуле: Q1 = -kSUx(x, t)∆t, где k — коэффициент теплопроводности материала ( = количеству тепла, протекающего в секунду через стержень единичной длины и единичной площади поперечного сечения при разности температур на противоположных концах, равной 1°). В этой формуле особого пояснения требует знак минус. Дело в том, что поток считается положительным, если он направлен в сторону увеличения х, а это, в свою очередь, означает, что слева от точки х температура больше, чем справа, то есть Ux CpS∆x∆U = kSUx(x + ∆х, t) ∆t — kSUx(x, t)∆t.

Если это равенство поделить на S∆x∆t и устремить ∆х и ∆t к нулю, то будем иметь:

Отсюда уравнение теплопроводности имеет вид

Ut = a 2 Uxx,
где — коэффициент температуропроводности.

В случае, когда внутри стержня имеются источники тепла, непрерывно распределенные с плотностью q(x,t), получится неоднородное уравнение теплопроводности

Начальные условия и граничные условия.

Для уравнения теплопроводности задается только одно начальное условие U|t=0 = φ(х) (или в другой записи U(x,0) = φ(х)) и физически оно означает, что начальное распределение температуры стержня имеет вид φ(х). Для уравнений теплопроводности на плоскости или в пространстве начальное условие имеет такой же вид, только функция φ будет зависеть, соответственно, от двух или трех переменных.

Граничные условия в случае уравнения теплопроводности имеют такой же вид, как и для волнового уравнения, но физический смысл их уже иной. Условия первого рода (5) означают, что на концах стержня задана температура. Если она не изменяется со временем, то g1(t) ≡ Т1 и g2(t) ≡ Т2, где Т1 и Т2 — постоянные. Если концы поддерживаются все время при нулевой температуре, то Т1= Т2 = 0 и условия будут однородными. Граничные условия второго рода (6) определяют тепловой поток на концах стержня. В частности, если g1(t) = g2(t) = 0, то условия становятся однородными. Физически они означают, что через концы не происходит теплообмен с внешней средой (эти условия еще называют условиями теплоизоляции концов). Наконец, граничные условия третьего рода (7) соответствуют случаю, когда через концы стержня происходит теплообмен с окружающей средой по закону Ньютона (напомним, что при выводе уравнения теплопроводности мы считали боковую поверхность теплоизолированной). Правда, в случае уравнения теплопроводности условия (7) записываются немного по-другому:

Физический закон теплообмена со средой (закон Ньютона) состоит в том, что поток тепла через единицу поверхности в единицу времени пропорционален разности температур тела и окружающей среды. Таким образом, для левого конца стержня он равен Здесь h1 > 0 — коэффициент теплообмена с окружающей средой, g1(t) — температура окружающей среды на левом конце. Знак минус поставлен в формуле по той же причине, что и при выводе уравнения теплопроводности. С другой стороны, в силу теплопроводности материала поток тепла через этот же конец равен Применив закон сохранения количества тепла, получим:

Аналогично получается условие (14) на правом конце стержня, только постоянная λ2 может быть другой, так как, вообще говоря, среды, окружающие левый и правый конец, бывают разные.

Граничные условия (14) являются более общими по сравнению с условиями первого и второго рода. Если предположить, что через какой-либо конец не происходит теплообмена со средой (то есть коэффициент теплообмена равен нулю), то получится условие второго рода. В другом случае предположим, что коэффициент теплообмена, например h1, очень большой.

Перепишем условие (14) при х = 0 в виде и устремим . В результате будем иметь условие первого рода:

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

Решение первой начально-краевой задачи для уравнения теплопроводности.

Рассмотрим однородную первую начально-краевую задачу для уравнения теплопроводности:

Найти решение уравнения

удолетворяющее граничным условиям

и начальному условию

Решим эту задачу методом Фурье.

Шаг 1. Будем искать решения уравнения (15) в виде U(x,t) = X(x)T(t).

Найдем частные производные:

Подставим эти производные в уравнение и разделим переменные:

По основной лемме получим

Теперь можно решить каждое из этих обыкновенных дифференциальных уравнений. Обратим внимание на то, что используя граничные условия (16), можно искать не общее решение уравнения б), а частные решения, удолетворяющие соответствующим граничным условиям:

Шаг 2. Решим задачу Штурма-Лиувилля

Эта задача совпадает с задачей Штурма-Лиувилля, рассмотренной в лекции 3. Напомним, что собственные значения и собственные функции этой задачи существуют только при λ>0.

Собственные значения равны

Собственные функции равны (См. решение задачи)

Шаг 3. Подставим собственные значения в уравнение а) и решим его:

Шаг 4. Выпишем частные решения уравнения (15):

В силу линейности и однородности уравнения (15) их линейная комбинация

Шаг 5. Определим коэффициенты An в (19), используя начальное условие (17):

Приходим к тому, что начальная функция φ(x) разлагается в ряд Фурье по собственным функциям задачи Штурма-Лиувилля. По теореме Стеклова такое разложение возможно для функций, удовлетворяющих граничным условиям и имеющих непрерывные производные второго порядка. Коэффициенты Фурье находятся по формулам

Вычислив эти коэффициенты для конкретной начальной функции φ(x) и подставив их значения в формулу (19), мы тем самым получим решение задачи (15), (16), (17).

Замечание. Используя формулу (19), можно также, как в лекции 3, получить решение первой начально-краевой задачи для уравнения Ut = a 2 Uxx. Оно будет иметь вид

где

Электронная библиотека

Рассмотрим однородный стержень (под стержнем в механике понимается тело с одним превалирующим линейным размером, например, столб можно рассматривать как стержень с переменным сечением) постоянного поперечного сечения S и длины l, теплоизолированный с боков, ось которого примем за ось Ох (рис. 4.3). Обозначим через U(x,t) температуру стержня в сечении с абсциссой х в момент времени t (предполагается, что во всех точках любого поперечного сечения стержня температура одна и та же.

Пусть – плотность стержня, – его удельная теплоемкость (это количество калорий, которое необходимо, чтобы единицу массы стержня нагреть на 1 ºС), – коэффициент теплопроводности (он представляет собой количество тепла в калориях, которое будет протекать за единицу времени через сечение стержня, если температура стержня падает на 1 ºС при перемещении вдоль стержня на единицу длины), – интенсивность теплового источника, находящегося в сечении х для момента t, отнесенная к единице массы и единице времени (т.е. количество тепла, создаваемого этим источником тепла за единицу времени и приходящегося на единицу массы; например, лампочку включенную в помещении, аппаратуру, работающую на космическом корабле и др. можно рассматривать как источники тепла;). Температуру же U(x, 0) стержня в начальный момент времени мы считаем известной и обозначим её через :

Тепловой режим на концах стержня может быть весьма разнообразен, т.е. температура может быть и постоянной, и изменяться, например, по законам: и др.). Мы рассмотрим два случая:

1) концы стержня поддерживаются при постоянной нулевой температуре:

2) концы стержня находятся в той же теплоизолирующей оболочке, что и весь стержень. Это означает, что через концы не происходит протекание тепла.

Рассмотрим сечение х нашего стержня и найдем, какое количество Q тепла протечет (слева направо) через это сечение за элементарный промежуток времени . В момент t температура стержня в точке х будет равна . Возьмем отличную от точки х точку x + dx стержня (пусть для определенности dx > 0) согласно закону Фурье количество тепла, протекающего в направлении оси Ох за бесконечно малый промежуток времени dt через сечение S с абсциссой х, будет:

где k – коэффициент теплопроводности ( представляет здесь величину градиента температуры U). Во формуле (4.37) стоит знак минус, так как при , т.е. при росте U вместе с х, поток тепла направлен в противоположную сторону, и наоборот).

Составим тепловой баланс для элемента , заключенную между бесконечно близкими сечениями х и x + dx. Предположим для определенности, что U возрастает в направлении Ох.

Тогда через сечение с абсциссой х тепло выходит, а через сечение с абсциссой x + dx входит. Пусть dQ есть количество тепла, накопленное нашим элементом за время dt. Тогда количество тепла, созданное за dt источниками тепла в элементе равно:

Используя формулу (4.37) будем иметь:

Применяя формулу Лагранжа из дифференциального исчисления с точностью до бесконечно малых высшего порядка малости, получим:

Следовательно, формула (4.39) примет вид:

С другой стороны, есть скорость изменения температуры элемента , и поэтому представляет собой изменение его температуры. Так как масса элемента равна , то накопленное при этом количество тепла будет:

Итак, сравнивая (4.41) и (4.40), после сокращения получим:

или, обозначая (величина а называется коэффициентом температуропроводности), окончательно получим дифференциальное уравнение, описывающее распределение температуры U(x,t) в стержне и носит название уравнения теплопроводности:

Если источники тепла отсутствуют, то уравнение (4.43) принимает вид уравнения свободного теплообмена в стержне:

Аналогичным образом можно вывести уравнения для распространения температуры в пластине:

где – оператор Лапласа (см. разд. 1.4.6).

Вывод: уравнение теплопроводности есть уравнение с частными производными второго порядка, параболического типа, так как ; в процессе вывода уравнения теплопроводности мы пришли к двум математическим задачам:

1) Найти то решение уравнения (4.44), которое удовлетворяет начальному условию (4.35) и граничным условиям (4.36).

2) Найти то решение уравнения (4.44), которое удовлетворяет начальному условию (4.35) и граничным условиям:

т.е. когда концы стержня теплоизолированы (количество тепла, протекающее через них, равно нулю).

Срочно?
Закажи у профессионала, через форму заявки
8 (800) 100-77-13 с 7.00 до 22.00

Дифференциальное уравнение теплопроводности для стержня

Температурное поле стержня или пластины может быть получено на основании решения соответствующего дифференциального уравнения для стержня или пластины вместе с граничными условиями. Такое уравнение обычно выводят специально, анализируя процесс переноса тепла в стержне или пластине. Однако оно может быть получено из общего уравнения Фурье, так как стержень или пластина представляют собою частные случаи тел с двух — и трёхмерным полем температур. Рассмотрим наиболее общий случай — анизотропное тело, все три размера которого имеют один порядок; пусть в теле равномерно распределены источники тепла удельной

мощностью W Вт/м. Стационарное поле температур в таком теле описывается дифференциальным уравнением (2.16) для параллелепипеда

Для ограниченного цилиндра в стационарном режиме можно использовать

уравнение (2.8), положив в нём — = 0:

Рассмотрим сначала температурное поле параллелепипеда, на границах у = ±1 и z = ±l2 которого справедливы условия (2.24) третьего рода

Обозначим t-tc = 9 и, в силу симметрии температурного поля, перепишем граничные условия в виде

Условия на границах х = ±1х пока не оговариваем. Предположим, что градиент температур в данном теле в направлениях у и z мал и покажем, как можно учесть эту особенность температурного поля в дифференциальном уравнении (2.16), которое представим в форме

, d29 , d2& , d29 ттг п

X —— + X —т — + X —— + W — о а*2 ‘ ay2 az2

Проведем над уравнением (2.93) следующую операцию:

[Л = ту f J fix, у, z)dydz,

т. е. умножим почленно уравнение на dydz и проинтегрируем его в пределах от 0 до 1у и /г, а затем разделим результат на произведение lylz.

9(х) = 1У1 [5(х, у, z) = — J j Э(х, у, z)dydz.

t rз a’i9, і rf3 a2f, , Vf/r^.

^ a7] = Tllr дґ^ = m[lwdy]dz

Рис. 2.8. К выводу дифференциальных уравнений параллелепипеда,

пластины и стержня

Поскольку в последнем выражении интегрирование и дифференцирование происходит по одной переменной у, то эти операции нельзя менять

местами. Рассмотрим интеграл

Последнее преобразование сделано на основании условия (2.92 а). Следовательно,

f Wx, К >z)dz = —»■Я (х), >Л 2 ду2 ІілЛ і,

&z(x) = — j&(x, ly, z)dz — значение осреднённой по z температуры на

Аналогично проведём операцию lyz над третьим членом уравнения (2.93). Используем условия (2.92):

9y(x) = yWx, y,h) 1иТгч1. Тогда выражение для Ь2

Уравнение (2.99), в котором Ь определяется по формуле (2.100), является дифференциальным уравнением для стержня с источником тепла.

Если стержень имеет периметр U и площадь поперечного сечения S, а теплообмен на границах происходит по закону Еіьютона с коэффициентом теплообмена а, то можно обосновать следующее выражение для b:

Тепломассообмен

Водяной теплый пол: преимущества системы

Систем обогрева жилища много, на любой вкус. Одним из наиболее востребованных в последнее время вариантов является водяной теплый пол http://ukrakvabud.com.ua/vodyanye-teplye-poly.html. Важнейшая особенность этой конструкции та, что нагревательные элементы устанавливаются не …

Cуперпозиции температурных полей

Для линейных задач теплопроводности справедлив принцип суперпозиции (сложения) температурных полей: температурное поле тела, которое формируется в результате нескольких тепловых воздействий, может быть представлено в виде алгебраической суммы температурных полей, вызванных …

Анализ ошибки измерения температуры

Ошибки за счёт оттока тепла по проводам термопары могут быть полностью устранены только в случае равенства температур пластины и окружающей среды. Такой предельный случай встречается редко, поэтому рассмотрим реальный процесс …

Продажа шагающий экскаватор 20/90

Цена договорная
Используются в горнодобывающей промышленности при добыче полезных ископаемых (уголь, сланцы, руды черных и
цветных металлов, золото, сырье для химической промышленности, огнеупоров и др.) открытым способом. Их назначение – вскрышные работы с укладкой породы в выработанное пространство или на борт карьера. Экскаваторы способны
перемещать горную массу на большие расстояния. При разработке пород повышенной прочности требуется частичное или
сплошное рыхление взрыванием.
Вместимость ковша, м3 20
Длина стрелы, м 90
Угол наклона стрелы, град 32
Концевая нагрузка (max.) тс 63
Продолжительность рабочего цикла (грунт первой категории), с 60
Высота выгрузки, м 38,5
Глубина копания, м 42,5
Радиус выгрузки, м 83
Просвет под задней частью платформы, м 1,61
Диаметр опорной базы, м 14,5
Удельное давление на грунт при работе и передвижении, МПа 0,105/0,24
Размеры башмака (длина и ширина), м 13 х 2,5
Рабочая масса, т 1690
Мощность механизма подъема, кВт 2х1120
Мощность механизма поворота, кВт 4х250
Мощность механизма тяги, кВт 2х1120
Мощность механизма хода, кВт 2х400
Мощность сетевого двигателя, кВ 2х1600
Напряжение питающей сети, кВ 6
Более детальную информацию можете получить по телефону (063)0416788


источники:

http://libraryno.ru/4-5-vyvod-uravneniya-teploprovodnosti-spec_gl_vm/

http://msd.com.ua/teplomassoobmen/differencialnoe-uravnenie-teploprovodnosti-dlya-sterzhnya/