Расчет конвейерной сверхвысокочастотной нагревательной установки для термообработки тонкопленочных и линейных материалов на основе прямоугольного волновода с Т-ребром

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

Рубрика Физика и энергетика
Вид курсовая работа
Язык русский
Дата добавления 29.12.2012
Размер файла 1,0 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

Размещено на http://www.allbest.ru/

Министерство образования и науки Российской Федерации

Саратовский государственный технический университет имени Гагарина Ю.А.

Кафедра радиотехники

Курсовая работа

По дисциплине Устройства СВЧ нагрева и обработки материалов

"Расчет конвейерной СВЧ нагревательной установки для термообработки тонкопленочных и линейных материалов на основе прямоугольного волновода с Т - ребром (ПВТР)"

Саратов 2012

Цель работы: Создание СВЧ нагревательных установок, конвейерных волноводных установок, на основе волноводов сложного сечения для равномерной обработки тонкослойного и линейного материала.

Исходные данные:

· мощность генератора 850 Вт

· камера на основе прямоугольного волновода с Т - ребром (ПВТР)

· рабочая частота 2450 МГц ± 2%

· скорость конвейерной установки v0 = 2см/с

· характеристики материала: L=0,3 м; W=5 мм.

Электро- и теплофизические параметры мышечной ткани говядины приведены в таблице 1.

Таблица 1 - Электро- и теплофизические параметры мышечной ткани говядины

t, °С

20

30

40

50

60

70

80

90

100

е'

53,5

51,8

50,0

48,0

46,5

45,0

43,0

42,0

40,0

е"

33,5

37,0

38,0

38,0

37,8

37,0

35,0

32,0

28,5

tgд

0,63

0,71

0,76

0,79

0,81

0,82

0,81

0,79

0,71

лт, Вт/м°С

0,44

0,43

0,42

0,45

0,47

0,48

0,51

0,52

0,54

ст, кг/м3

1003

1008

1012

1002

1002

1005

1004

1003

1002

Ст, Дж/кг·с

3607

3525

3442

3600

3481

3310

2865

2271

1667

ат, м3/с·10-4

1,22

1 22

1,22

1,25

1,35

1,45

1,78

2,29

3,24

tgд/е'

0,011

0,013

0,015

0,016

0,017

0,018

0,019

0,019

0,018

Стсm, Дж·кг/°С

3617

3546

3496

3600

3492

3329

2912

2292

1678

Таблица 2 - Электро- и теплофизические параметры полиметилметакрилата

t, °С

0

20

40

60

80

е'

3,8

3,9

3,98

4,2

4,4

е"

3,853

3,991

4,071

4,293

4,495

ст, кг/м3

1197

1190

1183

1177

1170

Ст, Дж/кг·с

1391

1427

1449

1466

1483

лт, Вт/м°С

0,194

0,195

0,197

0,199

0,2

ат, м3/с·10-4

0,1395

0,1367

0,136

0,1357

0,1349

tgд

1

1,03

1,05

1,11

1,16

Таблица 3 - Электро- и теплофизические параметры полиэтилена

t, °С

0

20

40

60

80

е'

3

3

3

3

3

е"

3,005

3,005

3,005

3,005

3,004

ст, кг/м3

955

950

944

938

932

Ст, Дж/кг·с

1707

2109

2376

2598

2819

лт, Вт/м°С

0,4

0,45

0,45

0,48

0,5

ат, м3/с·10-4

0,234

0,213

0,189

0,185

0,177

tgд

1

1

1

1

0,99

Таблица 4 - Электро- и теплофизические параметры картофеля

25

40

50

60

70

80

90

1

0,95

0,932

0,924

0879

0,841

0,818

0,492

0,571

0,618

0,623

0,621

0,64

0,602

1

1.161

1.256

1.266

1.262

1.301

1.224

1.505

1,169

1,345

1,327

1,532

1,682

1,677

1

0,995

0,988

0,981

0,953

0,953

0,951

1

0,777

0,894

0,882

1,018

1,118

1,114

Таблица 5 - Электро- и теплофизические параметры воды

20

30

40

50

60

70

80

90

100

1

0,987

0,942

0,897

0,839

0,795

0,737

0,692

0,603

0,177

0,143

0,113

0,098

0,084

0,071

0,064

0,048

0,017

1

0,807

0,638

0,554

0,474

0,401

0,361

0,271

0,096

1

0,998

0,998

0,998

0,999

1,001

1,003

1,006

1,009

1

0,997

0,994

0,99

0,985

0,98

0,974

0,967

0,96

998,2

995,7

992,2

988,1

983,1

977,8

971,8

965,3

958,4

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

Одним из существенных достоинств квазистационарных волноводов сложного поперечного сечения является наличие области в поперечном сечении, в которой электрическое поле доминантной волны однородно (область емкостного зазора). Размещение обрабатываемого материала в емкостном зазоре ВСС обеспечивает необходимое условие достижения однородной удельной плотности тепловых источников в сечении образца и соответственно создает условия равномерного его нагрева. Достижение однородной плотности тепловых источников по длине рабочей камеры в направлении распространения основной волны, путем соответствующего измерения продольной геометрии волновода, позволяет обеспечить однородную удельную плотность тепловых источников в объеме обрабатываемого материала. То есть на основе квазистационарных волноводов сложного поперечного сечения может быть создан класс СВЧ установок равномерного нагрева произвольных диэлектрических материалов, в том числе и материалов, электрофизические и тепловые параметры которых изменяются в процессе нагрева (термопараметрические материалы). При этом профильная геометрия рабочей камеры определяется требуемым законом изменения коэффициента затухания доминантной волны в направлении распространения, при котором обеспечивается qn=const в обрабатываемом материале, которое имеет следующий вид:

б(о)=-б0(1-о)/о, (1.1)

где а0 - коэффициент затухания основной волны во входном сечении волновода, о =z/L - приведенная безразмерная координата в направлении распространения волны; L - длина рабочей камеры СВЧ нагревательной установки. Соотношение (1.1) носит общий для волноводов характер и является основным критерием определения продольной геометрии рабочей камеры микроволновой установки с однородной продольной плотностью тепловых источников на основе решения прямой внутренней краевой задачи электродинамики для квазистационарных ВСС, частично заполненных произвольным диэлектрическим материалом. Заметим, что обеспечить выполнение условия (1.1) и достичь однородной плотности тепловых источников - qn=const путем изменения геометрии камеры в направлении распространения волны можно только на основной (доминантной) волне волновода. Условие (1.1) может быть выполнено путем уменьшения геометрии поперечного сечения в направлении распространения доминантной волны, причем характер этого изменен напрямую определяется соотношением (1.1). Это достаточно сложная задача, поскольку коэффициент затухания зависит от многочисленных параметров: геометрических параметров, определяющих формулу поперечного сечения волновода - (q1(b),...,qn(b)); заполняющего волновод материала (q1(m),...,qn(m)) электрофизических параметров обрабатываемого материала - е', tgд и динамикой изменения данных параметров в рабочем диапазоне температур:

б = б(q1(b),...,qn(b)); (q1(m),...,qn(m)); е'(t),tgд л0, (1.2)

где п и к - число геометрических параметров, определяющих форму поперечного сечения волновода и обрабатываемого материала, л0 - рабочая длина волны. Подробно расчет продольной геометрии рабочей камеры на основе прямоугольного волновода с Т-ребром, П и Н-волноводов на основе соотношений (1.1), (1.2). Данная задача по сути является обратной внутренней краевой задачи электродинамики для квазистационарных волноводов сложного поперечного сечения, частично заполненных произвольным диэлектрическим, поглощающим СВЧ мощность, материалом и для каждого материала, базового волновода и электротехнического процесса термообработки расчет проводится на основе решения ВКЗ ЭиТ с учетом упрощающих моментов, соответствующих данному процессу термообработки. Это довольно сложная задача и для своего решения требует учета всех особенностей электротехнологического процесса термообработки данного материала.

СВЧ нагревательные установки на основе квазистационарных волноводов сложного поперечного сечения подразделяются на установки стационарного и конвейерного типов. Установки стационарного типа, предназначенные для термообработки неподвижных диэлектрических материалов, недостаточно широко используются в технике и энергетике СВЧ в силу невозможности равномерной термообработки крупногабаритных, нестандартных образцов. Наибольшие перспективы имеют конвейерные установки для термообработки движущихся диэлектрических материалов таких, как влажные ткани, листовые материалы, тонкопленочные материалы, различные жидкости, смолы и т.д. В данной работе рассматриваются СВЧ нагревательные установки конвейерного типа на основе квазистационарных волноводов сложного поперечного сечения, имеющих четко выраженный емкостной зазор, в котором и производится термообработка листовых и жидких материалов. Рассматриваются волноводные конвейерные установки поперечного типа , обрабатываемый листовый материал в которых движется в направлении, перпендикулярном направлению распространения доминантной волны, а также продольного типа, в которых обрабатываемый материал перемещается в направлении распространения электромагнитной волны. Характерной особенностью данных установок является то, что распределение температурного поля в установках продольного типа носит нестационарный характер () в силу затухающего характера волны в направлении распространения, ВТО время как в установках поперечного типа тепловое поле в образце в области емкостного зазора (пространство взаимодействия) при постоянстве скорости протяжки материала является стационарным полем, но данное тепловое поле не является установившимся, поскольку время нахождения образца в емкостном зазоре много меньше времени установления температуры в обрабатываемом материале. При этом под временем фи , нахождения образца в пространстве взаимодействия понимается время, за которое отрезок обрабатываемого материала длиной L=d проходит область емкостного зазора ВСС со скоростью - х0и =d/х0 ).

Необходимо заметить, что в установках поперечного типа, тепловое поле в обрабатываемом материале в области емкостного зазора является стационарным, однако, если рассматривать температуру нагрева в фиксированной точке образца, то при прохождении данной точки емкостного зазора температура в ней будет нарастать во времени. При этом, если электрофизические и тепловые свойства обрабатываемого материала неизменны в процессе нагрева, то нарастание температуры нагрева обрабатываемого материала будет линейно. Это связано с тем, что в данном случае при однородности электрического поля в емкостном зазоре волновода удельная плотность тепловых источников в области расположения образца постоянна в направлении перемещения материала и, следовательно, при движении материала с постоянной скоростью температура нагрева в заданной точке материала будет нарастать по линейному закону. При изменении же электрофизических и тепловых параметров обрабатываемого материала в процессе нагрева, то есть при обработке термопластических материалов, температура нагрева будет изменяться по нелинейному закону как во времени, при рассмотрении теплового процесса в заданной точке образца, так и в направлении перемещения обрабатываемого материала (стационарное тепловое поле). Аналогичное положение наблюдается и в конвейерных установках продольного типа, теплое поле в которых стационарно () в направлении распространения электромагнитной волны. При термообработке в данных установках материалов, физические свойства которых изменяются в рабочем диапазоне температур, температура нагрева материала в направлении его перемещения будет изменяться по нелинейному закону, при этом характер этого изменения определяется зависимостью электрофизических и тепловых параметров обрабатываемого материала от температуры нагрева. Если же изменение геометрии рабочей камеры в направлении распространения волны обеспечивает выполнение условия (1.1), а физические свойства обрабатываемого материала неизменны в процессе нагрева, то изменение температуры образца в направлении его перемещения будет происходить по линейному закону. Основной задачей при этом является определение характера нелинейного изменения температуры нагрева при термообработке термопараметрических материалов, что принципиально важно для обеспечения заданного электротехнологического процесса термообработки. Так же как и для установок поперечного типа, температура в фиксированной точке обрабатываемого материала является нестационарной и нарастает в направлении перемещения образца, при этом время термообработки при постоянстве скорости перемещения обрабатываемого материала в установках продольного типа много больше времени - фн в установках продольного типа. То есть время термообработки в данном случае будет сравнимо со временем установления теплового процесса в данном нагреваемом материале. Кроме того, в установках продольного типа обрабатываемый материал в области поперечного сечения волновода однороден, в то время как в установках поперечного типа неоднороден, что связано с изменением электрофизических и тепловых параметров обрабатываемого материала в поперечном сечении квазистационарного волновода в направлении перемещения образца, при этом по длине физические свойства материала не изменяются. Данное положение влияет на методику расчета теплового поля соответствующего заданному электротехнологическому процессу термообработки конкретного диэлектрического, поглощающего СВЧ мощность, материала. При этом наиболее трудоемкую часть расчета составляет определение удельной плотности тепловых источников в объеме обрабатываемого термопараметрического материала, а также определение теплового поля при изменении электрофилических и тепловых параметров материала в поперечном сечении образца для установок поперечного типа и в направлении распространения доминантной волны в установках продольного типа для заданных экспериментальных характеристик - е'(t), tgд(t), c(t), с(t), л(t), a(t), б(t). Необходимо отметить, что в справочной литературе по электродинамике и теплопередаче указанные температурные зависимости приведены для очень узкого класса термопараметрических материалов, что значительно затрудняет создание волноводных СВЧ нагревательных установок конвейерного типа для равномерной термообработки широкого класса термопараметрических материалов с целью улучшения их выходных параметров и характеристик.

Рассмотрим процесс нагрева диэлектрической пластины, расположенной вертикально в центре емкостного зазора квазистационарного волновода прямоугольной формы (ПВТР, П или Н-волноводы). В данных волноводах электрическое поле однородно в области емкостного зазора, причем в нашем случае линии вектора напряженности электрического поля параллельны поверхности обрабатываемого материала и постоянны в области его поперечного сечения. Следовательно, удельная плотность тепловых источников в объеме обрабатываемого материала является постоянной величиной:

(1.2)

в силу того, что обрабатываемый материал однороден и его физические свойства неизменны во времени (е' = const; tgд = const; щ = const). Будем считать, что процесс нагрева является интенсивным, при котором в силу большой инерционности тепловых процессов можно пренебречь теплоотдачей с поверхности нагреваемого тела в металлические пластины емкостного зазора и в воздушную среду внутри волновода. При выполнении данного предположения температура в образце равномерно распределена по поверхности поперечного сечения. (У^=0, где V] - оператор Лапласа в поперечном сечении образца). Уравнение теплопроводности для данного случая будет иметь следующий вид:

, (1.3)

где ст и ст- удельные теплоемкость и плотность обрабатываемого материала. Решение уравнения (1.3) определяет линейное изменение температуры нагрева в обрабатываемом материале:

(1.4)

Соотношение может быть представлено через температурный напор (tн=t-tср) в виде:

, (1.5)

где е" = е'·tgд - параметр, характеризующий влияние ЭМГ мощности в образце.

Соотношения (1.4) - (1.5) определяют температуру нагрева в любой момент времени и конечную температуру нагрева на заданном интервале времени. Данные соотношения подтверждают тот факт, что при перемещении данного образца в направлении, перпендикулярном направлению распространения доминантной волны, температура в нем будет нарастать по линейному закону (установки поперечного типа) при постоянстве скорости перемещения образца (ф = У/?0). Для движущегося с постоянной скоростью материала соотношение (1.5) преобразуется к виду:

(1.6)

Данное соотношение справедливо для материалов, электрофизические и тепловые параметры которых неизменны в процессе термообработки. Естественно при температурной зависимости электрофизических - е'(t), tgд(t), тепловых - C(t),с(t) параметров от температуры нагрева соотношение (1.6) не действительно и температура обрабатываемого материала, как для фиксированной точки образца, так и в направлении движения изменяются по нелинейному закону. Определение характера изменения температуры нагрева представляет в общем случае чрезвычайно сложную задачу в технике и энергетике СВЧ.

Рассмотрим процесс нагрева термопараметрического диэлектрического поглощающего СВЧ мощность материала в конвейерной волноводной установке поперечного типа на основе прямоугольного волновода с Т-ребром (ПВТР). Также как и в рассмотренном выше случае, обрабатываемый материал транспортируется через центр емкостного зазора с постоянной скоростью ?0. Геометрия продольного профиля рабочей камеры обеспечивает выполнение условия (1.7), то есть на единице длины образца в направлении распространения доминантной волны выделяется постоянная электромагнитная энергия. То есть температура нагрева диэлектрической пластины изменяется лишь в направлении ее перемещения, то есть в направлении оси оY. При решении задачи теплопроводности пренебрежем передачей тепла в объеме обрабатываемого материала вследствие теплопроводности (q = -лТ·grad(t) = 0), а также теплоотдачей с поверхности нагреваемого материала в окружающую среду в пространстве взаимодействия, что обусловлено малым временем нахождения обрабатываемого материала в области емкостного зазора, а также использованием диэлектрических теплоизоляционных вставок, обладающих высоким тепловым сопротивлением (RТ = д/лТ "1, где д - толщина теплоизоляционной пластины, RТ - тепловое сопротивление). Относительная диэлектрическая проницаемость вставок - е'всm должна быть меньшее е' обрабатываемого материала, что в режиме распространения медленных волн в рабочей камере позволит повысить величину напряженности электрического поля в области расположения обрабатываемого материала, то есть удельную плотность тепловых источников в образце. Необходимо заметить, что термообработка термопараметрического материала проводится в электромагнитном поле доминантной волны, электрическое поле которой в области емкостного зазора однородно, причем электрическое поле в основном сосредоточено в емкостном зазоре.

Определим распределение электрического поля в нагреваемом термопараметрическом материале, что является необходимым условием нахождения плотности тепловых источников в объеме обрабатываемого материала. В отличие от термообработки однородных диэлектрических материалов, физические свойства которых не зависят от температуры нагрева, в конвейерных установках поперечного типа электрофизические и тепловые параметры обрабатываемого материала изменяются в направлении перемещения образца, что оказывает существенное влияние на распределение электрического поля и источники тепла - qv. Представим термопараметрический материал в виде слоистого материала -е' и tgд в каждом слое которого и постоянны, при этом в разных слоях электрофизические параметры различны. Данная модель позволяет смоделировать нагрев термопараметрического материала, то есть определить е' и tgд в каждом слое, соответствующие заданному изменению температуры обрабатываемого материала в направлении его перемещения и, соответственно, экспериментальному изменению - е'(t) и tgд(t) в процессе термообработки. Поскольку нормальная составляющая вектора напряжен-ности электрического поля на границе раздела сред терпит разрыв, то для определения характера изменения электрического поля в емкостном зазоре воспользуемся граничным условием для вектора электрической индукции (нормальная составляющая) на границе раздела каждого слоя:

, (1.7)

где Dn - нормальная составляющая вектора электрической индукции, k - число слоев в модели.

Из соотношения (1.7) следует следующее выражение для компонент вектора напряженности электрического поля в образце -Еу:

, (1.8)

где ; ; - абсолютная диэлектрическая проницаемость каждого диэлектрического слоя; ,, - составляющая вектора напряженности электрического поля в каждом слое в направлении перемещении обрабатываемого материала. Соотношение (1.8). позволяет определить составляющую напряженности электрического поля через составляющую поля произвольного слоя. В качестве эталонного слоя выберем слой диэлектрика той же толщины, находящийся на входе емкостного зазора (, - абсолютное значение диэлектрической проницаемости при температуре среды; соответственно и - значения нормальных составляющих векторов электрической индукции и напряженности электрического поля при температуре ). Нормировку соотношения (1.8) проведем относительно слоя диэлектрика на входе рабочей камеры:

(1.9)

Соотношение (1.9) определяет распределение Еу составляющей вектора напряженности электрического поля в термопараметрическом материале в области емкостного зазора (при d00 , где d0 - толщина каждого из диэлектрических слоев

(1.10)

,

где -напряженность электрического поля в емкостном зазоре, частично заполненного однородным диэлектрическим материалом с расположенным между обкладками конденсаторного зазора. Соотношение (1.10) позволяет определить удельную плотность тепловых источников в объеме термопараметрического материала:

, (1.11)

где ; и относительная диэлектрическая проницаемость однородного диэлектрического материала соответствующая и термопараметрического материала на входе в рабочую камеру; - функция, определяющая изменение от температуры нагрева:

, (1.12)

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

экспериментально; q0- удельная плотность тепловых источников в объеме обрабатываемого материала, на входе в область взаимодействия электромагнитного поля с термопараметричеcим материалом, которая определяется соотношением:

(1.13)

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

(1.14)

и, соответственно, соотношение (1.11) с помощью выражения (1.14) можно представить в следующем виде:

(1.15)

Тогда уравнение теплопроводности для движущегося с постоянной скоростью слоя диэлектрика (то есть температура в слое изменяется во времени) запишется в виде:

(1.16)

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

(1.17)

(1.18)

(1.19)

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

(1.20)

Интегрируя данное уравнение в области емкостного зазора, получим следующее выражение:

,

Уравнение (1.20) для заданных электрофизических и тепловых свойств конкретного термопараметрического материала может быть решено различными итерационными методами [1-3]. Необходимо учитывать, что температурные параметры-cT(t), pT(t) менее зависимы от температуры нагрева в рабочем диапазоне температур, нежели электрофизические параметры. Это касается твердых материалов, жидких сред и газов. Наиболее сильное изменение электрофизических и тепловых параметров наблюдается в таких электротехнологических процессах, как сушка влажных материалов, вулканизация резино-технических изделий, высокотемпературная обработка и спекание специальных брикетированных материалов. Наиболее частый случай может быть описан следующим соотношением:

(1.21)

В данной диссертационной работе будет проведено решение уравнения теплопроводности и интегральных уравнений (1.20) и (1.21) для конкретных термопараметрических материалов одним из приближенных методов. Необходимо отметить, что при расчете теплового поля в емкостном зазоре конвейерной СВЧ нагревательной установки поперечного типа на основе квазистационарных волноводов, в которых обрабатываемый материал перемещается перпендикулярно направлению распространения доминантной волны на основе температурных зависимостей электрофизических и тепловых параметров, метод решения и последовательность итерации привязана к функции , которая для термопараметрических материалов изменяется нелинейно.

Воспользуемся для решения уравнения (1.20) линейной аппроксимацией функции . Будем считать, что на малом температурном интервале функция изменяется по линейному закону:

, (1.22)

где -угловой коэффициент (тангенс угла касательной в данной точке кривой к оси абсцисс) на первом этапе итерационного процесса решение уравнения теплопроводности, в точке . Для i-го итерационного этапа решение уравнения (1.20), линейная аппроксимация функции может быть представлена в виде:

, (1.23)

где =const и определяется значением tк на i-ом интервале разбиения искомой функции. Для упрощения процесса решения обозначим входящую в уравнение (1.23) функцию при tkl, через . Тогда соотношение (1.23) можно записать в следующем виде:

(1.24)

Подставляя соотношение (1.23) в интегральное выражение (1.20), получим:

(1.25)

Проводя интегрирование в левой части соотношения (1.25) и учитывая, что = const и = const, получим:

(1.26)

удовлетворяя уравнение пределам интегрирования, получим квадратное алгебраическое уравнение относительно температуры нагрева на i-ом интервале итерации функции :

, (1.27)

(1.28)

Решение уравнения (1.27) будет иметь следующий вид:

(1.29)

Соотношение (1.29) справедливо для i-го интервала аппроксимации функции и соответствующего i-го диапазона изменения температуры нагрева . Решение уравнения (1.30) согласно соотношению (1.29) для первого итерационного этапа будет иметь вид:

(1.30)

Учитывая, что на первом итерационном шаге решение уравнения теплопроводности , то при малых значениях временного интервала , величина . Используя для подкоренного выражения решение (1.30) представление бинома Ньютона [4] при u<<1, получим окончательное выражение для температуры нагрева:

(1.31)

Это вполне естественный результат, поскольку на первом итерационном этапе решения уравнения теплопроводности для термопараметрического материала, нагреваемого в квазистационарном волноводе величина определяется в точке y=0, в окрестности которой температура изменяется линейно. При дальнейшем движении обрабатываемого материала будут проявляться нелинейные эффекты в распределении теплового поля в образце.

Предложенный в работе метод решения внутренней краевой задачи теплопроводности для термопараметрических материалов, термообработка которых осуществляется в конвейерных установках поперечного типа на основе квазистационарных волноводов, в которых нагреваемый материал пропускается в центре емкостного зазора, то есть в области существования однородного электрического поля, является наиболее универсальным и позволяющим получить распределение теплового поля в области пространства взаимодействия с высокой степенью точности. При этом уменьшение температурного интервала при линейной итерации функции , определяющей изменение электрофизических и тепловых параметров обрабатываемого материала в рабочем диапазоне температур, позволяет повысить точность вычислений. Данный подход определения распределения температуры нагрева в термопараметрическом материале требует определения температурных зависимостей электрофизических и тепловых параметров в рабочем диапазоне температур. Это позволит наиболее полно реализовать заданный электротехнологический процесс термообработки конкретного диэлектрического, поглощающего СВЧ мощность, материала. Заметим, что конвейерные СВЧ нагревательные установки волноводного типа, предназначенные для нагрева термопараметрических материалов, не являются установками универсального действия, поскольку изменение электрофизических и тепловых параметров в процессе термообработки накладывают жесткие условия на геометрию рабочей камеры установки, скорость протяжки диэлектрического материала, уровень подводимой в камеру СВЧ мощности. Кроме того, существенным ограничительным фактором является то, что термообработка термопараметрических материалов проводится в электромагнитном поле доминантной (основной) волны квазистационарных волноводов сложного поперечного сечения.

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

В данной работе предложен и апробирован эффективный способ определения теплового поля в термопараметрическом материале, транспортировка которого осуществляется в центре емкостного зазора квазистационарного волновода прямоугольного внешнего профиля с постоянной скоростью . Как было показано выше, в конвейерных установках поперечного типа изменение температуры нагрева образца происходит в направлении перемещения обрабатываемого материала. Продольная однородность тепловых источников и равномерность нагрева листового материала обеспечивается соответствующим изменением продольной геометрии рабочей камеры, при которой обеспечивается выполнение условия (1.21). В связи с этим обрабатываемый материал можно представить как многослойный образец, температура нагрева в каждом из слоев постоянна и ее изменение происходит только в направлении оси oY, то есть в направлении перемещения. Толщину образца (каждого слоя) будем считать постоянной и равной -. На входе в рабочую камеру (пространство взаимодействия) температура поверхности образца - . Будем считать, что в первом слое удельная плотность тепловых источников определяется соотношением:

, (1.32)

где L - ширина, W - толщина и d - высота обрабатываемого материала. - поглощенная образцом СВЧ - мощность при за единицу времени. Заметим, что величина определяется экспериментально для неподвижного материала. Изменение температуры нагрева по толщине слоя определяется (1.31) (линейный закон изменения температуры нагрева в i - м слое), которое может быть представлено в виде:

, (1.33)

где . Для первого этапа итерационного процесса решения краевой задачи теплопроводности получим:

(1.34)

при этом ( определяется соотношением (1.32)).

Подставляя в соотношение (2.54) , получим температуру , являющуюся конечной температурой первого слоя, которая определяет значение на втором этапе решения задачи на основе кривых , для температуры tlk . Решение на втором этапе, таким образом, определяется следующим соотношением:

(1.35)

Подставляя в (1.35) значение ,получим значение , определяющее значение на последующем этапе решения. Проводя последовательно решение для каждого слоя, получим уравнение, определяющее изменение температуры нагрева на последнем, n-ом, этапе решения задачи теплопроводности:

, (1.36)

что, соответственно, позволяет определить конечную температуру на выходе образца . Таким образом, по значениям - tcp,tl,t2..... tn мы построим искомую кривую, определяющую изменение температуры нагрева в термопараметрическом материале в области емкостного зазора в СВЧ нагревательной установке поперечного типа. Увеличение числа слоев в исследуемом образце позволяет повысить точность определения - образца, при этом время термообработки в каждом случае определяется величиной - , а всего образца - , то есть:

(1.37)

Распределение температуры образца в области емкостного зазора, у которого , а величины электрофизических параметров изменяются линейно в рабочем диапазоне температур:

(1.38)

При этом значения данных параметров при температуре среды (начальное условие) равны: . Рассмотрены три случая: 1) q1>q2; 2) ; 3) q1<q2. Легко видеть, что в случаях 1 и 3 температура t(y) в области емкостного зазора в термопараметрическом материале изменяется по нелинейному закону, а в случае 2 - по линейному закону, также как и для материала, электрофизические и тепловые параметры которого неизменны в процессе термообработки. Данные результаты получены описанным выше способом.

Необходимо отметить, что точность определения теплового поля в области расположения обрабатываемого термопараметрического материала во многом зависит от точности определения удельной плотности тепловых источников в первом слое - q0. Уточнение q0 можно провести следующим образом - для полученного значения определяется средняя температура нагрева первого слоя:

(1.39)

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

Расчет.

Данные приведенные в таблице 1 позволяют однозначно определить функцию Ш(t) и соответственно функцию Ш1(t)=l/Ш(t), которые определяют температурное изменение электрофизических ц1(t), ц2(t) и тепловых ц3(t), ц4(t) параметров от температуры нагрева образца. Данные зависимости однозначно определяют температуру обрабатываемого материала в емкостном зазоре в направлении его перемещения, а так же начальные условия при расчете продольного профиля рабочей камеры, обеспечивающей однородную удельную плотность тепловых источников в направлении распространения доминантной волны квазистационарного волновода, то есть равномерную термообработку конкретного термопараметрического материала. Для нахождения температурных зависимостей Ш(t) и Ш1(t) определим температурные функции электрофизических и тепловых параметров - ц1(t), ц2(t), ц3(t), ц4(t), а также основополагающие функции Ш(t) и Ш1(t), которые позволяют определить распределение теплового поля в термопараметрическом материале в области емкостного зазора рабочей камеры установки.

Таблица 6 - полиметилметакрилат

t,?C

0

20

40

60

80

ц1(t)

1

1,02

1,04

1,105

1,15

ц2(t)

1

1,03

1,05

1,11

1,16

ц4(t)

1

0,979

0,973

0,97

0,965

ц3(t)

1

1,025

1,04

1,05

1,06

Ш(t)

1

1,006

0,99

0,9866

0,963

Таблица 7 - полиэтилен

t,?C

0

20

40

60

80

ц1(t)

1

1

1

1

1

ц2(t)

1

1

1

1

0,99

ц4(t)

1

0,91

0,8

0,79

0,75

ц3(t)

1

1,23

1,39

1,52

1,65

Ш(t)

1

0,893

0,899

0,833

0,8

Таблица 8 - филе мяса говядины

t,?C

20

30

40

50

60

70

80

90

100

ц1(t)

1

0,97

0,93

0,9

0,87

0,84

0,8

0,78

0,75

ц2(t)

1

1,13

1,21

1,25

1,29

1,3

1,3

1,29

1,13

ц4(t)

1

10,05

1,01

1

1

1

1

1

1

ц3(t)

1

0,98

0,95

0,99

0,97

0,92

0,79

0,62

0,46

Ш(t)

1

1,19

1,37

1,45

1,54

1,69

2,04

2,7

3,33

Таблица 9 - картофель

t,?C

25

40

50

60

70

80

90

Ш(t)

1

1,581

1,526

1,584

1,48

1,452

1,412

Таблица 10 - вода

t,?C

25

40

50

60

70

80

90

Ш(t)

1

1,581

1,526

1,584

1,48

1,452

1,412

Формулы для расчёта функции Ш(t) для картофеля и воды:

На основе таблиц 3и 6 построим графики функции Ш(t) и Ш1(t) для полиметилметакрилата и полиэтилена соответственно.

Рисунок 1 - График функции Ш(t) для полиметилметакрилата

Рисунок 2 - График функции Ш1(t) для полиметилметакрилата

Рисунок 3 - График функции Ш(t) для полиэтилена

Рисунок 4 - График функции Ш1(t) для полиэтилена

Рисунок 5 - График функции Ш(t) для филе мяса говядины

Рисунок 6 - График функции Ш1(t) для филе мяса говядины

Рисунок 7 - График функции Ш1(t) для картофеля

сверхвысокочастотный нагревательный электродинамика теплопроводность

Рисунок 8 - График функции Ш1(t) для картофеля

Решать будем, используя функцию Ш(t). Представим функцию Ш(t) с угловым коэффициентом -q1*(t):

Ш(t)=1+ q1*t(ф)

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

Здесь под t будем понимать температурный напор.

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

Так как мы используем функцию Ш(t), то здесь можно выделить два итерационных этапа:

первый температурный напор равен 54;

а второй соответственно 26.

Рассчитаем чему равен коэффициент -q1*(t):

Найдем обратную величину:

Примем, что , а , тогда перепишем формулу для температурного напора:

Приравняем о = 1, тогда:

Для первого итерационного участка имеем:

tн =72(ех -1) = 54

х -1) = 54/72= 0,75

ех= 1,75

x= 0,56

В итоге имеем:

tн =72(е0,56·о -1)

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

ех

1,058

1,12

1,18

1,25

1,32

1,4

1,48

1,57

1,66

1,75

ех -1

0,058

0,12

0,18

0,25

0,32

0,4

0,48

0,57

0,66

0,75

72(е0,56·о -1)

4,32

8,64

13

18

23

28,8

34,6

41

47,5

54

tн

24,32

28,64

33

38

43

48,8

54,6

61

67,5

74

Рассчитаем чему равен коэффициент :

Найдем обратную величину:

Температурный напор на этом итерационном участке:

tн =74+16,7(ехо -1) = 26

Зная, что о = 1, получим:

tн =16,7(ех -1) = 26

х -1) = 26/16,7= 1,557

ех= 2,557

x= 0,94

В итоге имеем:

tн =16,7(е0,94·о -1)

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

ех

1,099

1,2

1,33

1,46

1,6

1,758

1,9

2,12

2,33

2,56

ех -1

0,099

0,2

0,33

0,46

0,6

0,758

0,9

1,12

1,33

1,56

16,7(е0,94·о -1)

1,65

3,45

5,44

7,62

10,02

12,65

15,546

18,725

22,2

26

tн

75,6

77,45

79,4

81,6

84

86,65

89,5

92,7

96,2

100

Теперь найдем где у нас будет изгиб, т.е. найдем d1 и d2

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

Составим систему уравнений:

d1=2,6· d2

d1+ d2=6,6

Решая ее, получим, что d1 =4,8мм, d2 =1,8мм.

Составим таблицу, на основе которой построим график:

d

0,48

0,96

1,44

1,92

2,4

2,88

3,36

3,84

4,32

4,8

t

24,3

28,6

33

38

43

48,8

54,6

61

67,5

74

d

4,98

5,16

5,34

5,52

5,7

5,88

6,06

6,24

6,42

6,6

t

75,6

77,45

79,4

81,6

84

86,65

89,5

92,7

96,2

100

Тем самым получаем график температуры:

На основе полученного графика делаем вывод, что температура изменяется не линейно.

Рассчитаем объем прохождения мышечной ткани говядины:

V = L Wd = 0,3 ·5 · 10- 3 · 6,6 · 10- 3 = 9,9 · 10-6 м3 = 9,9см3

здесь L - ширина(м), W - толщина(м), d - высота обрабатываемого материала(м).

В нашей установке стоит магнетрон мощностью 850 Вт, примем, что наш образец поглощает 600 Вт.

Найдем количества тепла, необходимого для нагрева Q, используя

Ст,ср =3096 Дж/кг·с; ст,ср=1004 кг/м3

Q= Ст·с·V·Дt = 2462 Дж

Найдем время, в течение которого происходит нагрев всего образца:

Таким образом, получается, что наш образец будет проходить область емкостного зазора за 0,33 секунды.

Построим распределение уровня поглощенной мощности в рабочей камере, обеспечивающей q1 = const в объеме обрабатываемого материала.

Распределение Рпот(о),Вт по длине рабочей камеры при различных б

Р0 = 600 Вт

Р = Р0· е-2б Вт

б = 1

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1,0

е-о

0,905

0,819

0,741

0,67

0,607

0,549

0,497

0,449

0,407

0,368

Рпот(о),Вт

543

491

444

402

364

325

298

270

244

221

ДРпот(о),Вт

57

52

47

42

38

35

31

28

26

23

б = 2

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1,0

е-

0,819

0,67

0,549

0,449

0,368

0,301

0,247

0,202

0,165

0,135

Рпот(о),Вт

491

402

329

270

221

181

148

121

99

81

ДРпот(о),Вт

109

89

73

59

49

40

33

27

22

18

б = 3

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1,0

е-

0,741

0,549

0,407

0,301

0,223

0,165

0,122

0,091

0,067

0,05

Рпот(о),Вт

444

329

244

181

134

99

73

54

40

30

ДРпот(о),Вт

156

115

85

63

47

35

26

19

14

10

Проведем аналогичный расчет, выделив один итерационный этап для функции Ш(t) картофеля.

Рисунок 11 - Температурная зависимость функции Ш(t) для картофеля

Проведем расчет для температурной точки 25? С

Рассчитаем чему равен коэффициент -q1*(t):

Найдем обратную величину:

Примем, что , а , тогда перепишем формулу для температурного напора:

Приравняем о = 1, тогда:

Для одного итерационного участка имеем:

tн =25(ех -1) = 25

х -1) = 1

ех= 2

x= 0,693

В итоге имеем: tн =25(е0693,·о -1)

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

tн

1,794

3,717

5,777

7,986

10,353

12,89

15,608

18,522

21,645

24,993

Теперь найдем d1 :

Получим, что d1=6,6мм.

Составим таблицу, на основе которой построим график:

d

0,66

1,32

1,98

2,64

3,3

3,96

4,62

5,28

5,94

6,6

t

1,794

3,717

5,777

7,986

10,353

12,89

15,608

18,522

21,645

24,993

Произведем аналогичный расчет для точки 90? С

Рассчитаем чему равен коэффициент -q1*(t):

Найдем обратную величину:

Примем, что , а , тогда перепишем формулу для температурного напора:

Приравняем о = 1, тогда:

Для одного итерационного участка имеем:

tн =60.7(ех -1) = 90

х -1) = 1.48

ех= 2.48

x= 0,908

В итоге имеем:

tн =60.7(е0,908·о -1)

о

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

tн

5.77

12.087

19.006

26.582

34.878

43.963

53.911

64.805

76.734

89.797

Теперь найдем d1 :

Получим, что d1=6,6мм.

Составим таблицу, на основе которой построим график:

d

0,66

1,32

1,98

2,64

3,3

3,96

4,62

5,28

5,94

6,6

t

5.77

12.087

19.006

26.582

34.878

43.963

53.911

64.805

76.734

89.797

Тем самым получаем график температур:

Продольный профиль рабочей камеры СВЧ нагревательной установки поперечного типа на основе ВСС.

а) продольный профиль установки на ПВТР для а -1;

б) продольный профиль установки на ПВТР для а = 2 и а = 3;

в) поперечное сечение рабочей камеры на ПВТР.

Рис. а) и б) - вид сверху (двухстороннее возбуждение).

Заключение

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

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

3. Показано, что в конвейерных установках на базе квазистационарных волноводов сложного поперечного сечения можно осуществить равномерную термообработку материалов, электрофизические и тепловые свойства которых изменяются в процессе нагрева, при расположении термопараметрического материала в емкостном зазоре волновода и соответствующем изменении продольной геометрии рабочей камеры установки, что позволит значительно расширить область применения установок с бегущей волной в СВЧ-энергетике.

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

1. Коломейцев В.А. Взаимодействие электромагнитных волн с поглощающими средами и специальные СВЧ-системы равномерного нагрева.// Дисс. на соискание учёной степени д.т.н., Саратов, 1999г., 432с.

2. Железняк А.Р. СВЧ-устройства на основе волноводов сложного поперечного сечения для равномерного нагрева диэлектрических материалов.// Дисс. на соискание учёной степени к.т.н., Саратов, 2001г.

3. Цыганков А.В. Электротехнологические СВЧ установки равномерного нагрева диэлектрических материалов на основе волноводов сложных сечений.// Дисс. на соискание учёной степени к.т.н., Саратов, 2003г.

4. Коломейцев В.А., Яковлев В.В. Расчёт электромагнитных полей рабочей камеры СВЧ нагревательной установки на П- волноводе .// Изв . Вузов Сер. Радиоэлектроника , 1987г.-N.9. с.65-66.

5. Коломейцев В.А., Яковлев В.В. Синтез прямоугольного волновода с Т- ребром - Саратов: СГТУ, 1997г.

6. Коломейцев В.А., Комаров В.В. Расчёт параметров базовых элементов рабочих камер СВЧ нагревательных установок на волноводах сложных сечений. // Современные проблемы применения СВЧ- энергии: Сб. тр. науч.- техн. конф,- Саратов, 1993г.- С.61-62.

Размещено на Allbest.ru


Подобные документы

  • Уравнение теплопроводности: физический смысл и выводы на примере линейного случая. Постановка краевой задачи остывания нагретых тел, коэффициент теплопроводности. Схема метода разделения переменных Фурье применительно к уравнению теплопроводности.

    курсовая работа [245,8 K], добавлен 25.11.2011

  • Основные физические принципы волноводной фотоники. Классификация оптических волноводов. Геометрическая оптика планарных волноводов. Классификация мод планарного волновода. Волноводные моды тонкопленочного волновода. Эффективная толщина волновода.

    реферат [2,0 M], добавлен 16.06.2019

  • Особенность волновода как направляющей системы. Решение задачи распространения волн в волноводе круглого сечения с физической точки зрения. Структура поля в плоскости продольного сечения. Применение волны H01 круглого волновода для дальней связи.

    курсовая работа [279,6 K], добавлен 25.06.2013

  • Дифференциальное уравнение теплопроводности. Поток тепла через элементарный объем. Условия постановка краевой задачи. Методы решения задач теплопроводности. Численные методы решения уравнения теплопроводности. Расчет температурного поля пластины.

    дипломная работа [353,5 K], добавлен 22.04.2011

  • Содержание закона Фурье. Расчет коэффициентов теплопроводности для металлов, неметаллов, жидкостей. Причины зависимости теплопроводности от влажности материала и направления теплового потока. Определение коэффициента теплопередачи ограждающей конструкции.

    контрольная работа [161,2 K], добавлен 22.01.2012

  • Устройство прямоугольного объемного резонатора. Структура электромагнитного поля. Общая задача о собственных колебаниях в прямоугольном объемном резонаторе. Понятие основного типа колебаний. Структура электромагнитного поля в прямоугольном резонаторе.

    курсовая работа [356,3 K], добавлен 13.05.2011

  • Выбор размеров поперечного сечения волновода. Определение максимальной и пробивной мощности, затухания и длины волн, фазовой и групповой скорости волновода, характеристического сопротивления. Установление частотного диапазона, в котором можно работать.

    курсовая работа [6,0 M], добавлен 10.12.2012

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

    курсовая работа [1,1 M], добавлен 20.01.2014

  • Величина коэффициента и единица измерения теплопроводности. Расчет теплоотдачи у наружной поверхности ограждения. Сущность теплового излучения. Удельная теплоёмкость материала, её зависимость от влажности. Связь теплопроводности и плотности материала.

    контрольная работа [35,3 K], добавлен 22.01.2012

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

    курсовая работа [261,3 K], добавлен 20.02.2013

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.