Расчет конвейерной сверхвысокочастотной нагревательной установки для термообработки тонкопленочных и линейных материалов на основе прямоугольного волновода
Необходимо заметить, что в установках поперечного типа, тепловое поле в обрабатываемом материале в области емкостного зазора является стационарным, однако, если рассматривать температуру нагрева в фиксированной точке образца, то при прохождении данной точки емкостного зазора температура в ней будет нарастать во времени. При этом, если электрофизические и тепловые свойства обрабатываемого… Читать ещё >
Расчет конвейерной сверхвысокочастотной нагревательной установки для термообработки тонкопленочных и линейных материалов на основе прямоугольного волновода (реферат, курсовая, диплом, контрольная)
Министерство образования и науки Российской Федерации Саратовский государственный технический университет имени Гагарина Ю.А.
Кафедра радиотехники.
Курсовая работа По дисциплине Устройства СВЧ нагрева и обработки материалов.
" Расчет конвейерной СВЧ нагревательной установки для термообработки тонкопленочных и линейных материалов на основе прямоугольного волновода с Т — ребром (ПВТР)" .
Саратов 2012.
Цель работы: Создание СВЧ нагревательных установок, конвейерных волноводных установок, на основе волноводов сложного сечения для равномерной обработки тонкослойного и линейного материала.
Исходные данные:
· мощность генератора 850 Вт.
· камера на основе прямоугольного волновода с Т — ребром (ПВТР).
· рабочая частота 2450 МГц ± 2%.
· скорость конвейерной установки v0 = 2см/с.
· характеристики материала: L=0,3 м; W=5 мм.
Электрои теплофизические параметры мышечной ткани говядины приведены в таблице 1.
Таблица 1 — Электрои теплофизические параметры мышечной ткани говядины.
t, °С. | ||||||||||
е'. | 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. | ||||||||||
Ст, Дж/кг· с. | ||||||||||
ат, м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, Дж· кг/°С. | ||||||||||
Таблица 2 — Электрои теплофизические параметры полиметилметакрилата.
t, °С. | ||||||
е'. | 3,8. | 3,9. | 3,98. | 4,2. | 4,4. | |
е" . | 3,853. | 3,991. | 4,071. | 4,293. | 4,495. | |
ст, кг/м3. | ||||||
Ст, Дж/кг· с. | ||||||
лт, Вт/м°С. | 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,03. | 1,05. | 1,11. | 1,16. | ||
Таблица 3 — Электрои теплофизические параметры полиэтилена.
t, °С. | ||||||
е'. | ||||||
е" . | 3,005. | 3,005. | 3,005. | 3,005. | 3,004. | |
ст, кг/м3. | ||||||
Ст, Дж/кг· с. | ||||||
лт, Вт/м°С. | 0,4. | 0,45. | 0,45. | 0,48. | 0,5. | |
ат, м3/с· 10-4. | 0,234. | 0,213. | 0,189. | 0,185. | 0,177. | |
tgд. | 0,99. | |||||
Таблица 4 — Электрои теплофизические параметры картофеля.
0,95. | 0,932. | 0,924. | 0,841. | 0,818. | ||||
0,492. | 0,571. | 0,618. | 0,623. | 0,621. | 0,64. | 0,602. | ||
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. | ||
0,995. | 0,988. | 0,981. | 0,953. | 0,953. | 0,951. | |||
0,777. | 0,894. | 0,882. | 1,018. | 1,118. | 1,114. | |||
Таблица 5 — Электрои теплофизические параметры воды.
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. | ||
0,807. | 0,638. | 0,554. | 0,474. | 0,401. | 0,361. | 0,271. | 0,096. | |||
0,998. | 0,998. | 0,998. | 0,999. | 1,001. | 1,003. | 1,006. | 1,009. | |||
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) представление бинома Ньютона при 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) q12
. Легко видеть, что в случаях 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. | ||||||
ц1(t). | 1,02. | 1,04. | 1,105. | 1,15. | ||
ц2(t). | 1,03. | 1,05. | 1,11. | 1,16. | ||
ц4(t). | 0,979. | 0,973. | 0,97. | 0,965. | ||
ц3(t). | 1,025. | 1,04. | 1,05. | 1,06. | ||
Ш (t). | 1,006. | 0,99. | 0,9866. | 0,963. | ||
Таблица 7 — полиэтилен.
t,?C. | ||||||
ц1(t). | ||||||
ц2(t). | 0,99. | |||||
ц4(t). | 0,91. | 0,8. | 0,79. | 0,75. | ||
ц3(t). | 1,23. | 1,39. | 1,52. | 1,65. | ||
Ш (t). | 0,893. | 0,899. | 0,833. | 0,8. | ||
Таблица 8 — филе мяса говядины.
t,?C. | ||||||||||
ц1(t). | 0,97. | 0,93. | 0,9. | 0,87. | 0,84. | 0,8. | 0,78. | 0,75. | ||
ц2(t). | 1,13. | 1,21. | 1,25. | 1,29. | 1,3. | 1,3. | 1,29. | 1,13. | ||
ц4(t). | 10,05. | 1,01. | ||||||||
ц3(t). | 0,98. | 0,95. | 0,99. | 0,97. | 0,92. | 0,79. | 0,62. | 0,46. | ||
Ш (t). | 1,19. | 1,37. | 1,45. | 1,54. | 1,69. | 2,04. | 2,7. | 3,33. | ||
Таблица 9 — картофель.
t,?C. | ||||||||
Ш (t). | 1,581. | 1,526. | 1,584. | 1,48. | 1,452. | 1,412. | ||
Таблица 10 — вода.
t,?C. | ||||||||
Ш (t). | 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,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. | 28,8. | 34,6. | 47,5. | ||||||
tн. | 24,32. | 28,64. | 48,8. | 54,6. | 67,5. | ||||||
Рассчитаем чему равен коэффициент :
Найдем обратную величину:
Температурный напор на этом итерационном участке:
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,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. | ||
tн. | 75,6. | 77,45. | 79,4. | 81,6. | 86,65. | 89,5. | 92,7. | 96,2. | |||
Теперь найдем где у нас будет изгиб, т. е. найдем 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. | 48,8. | 54,6. | 67,5. | ||||||
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. | 86,65. | 89,5. | 92,7. | 96,2. | |||
Тем самым получаем график температуры:
На основе полученного графика делаем вывод, что температура изменяется не линейно..
Рассчитаем объем прохождения мышечной ткани говядины:
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бLо Вт.
б = 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. | |
Рпот(о), Вт. | |||||||||||
ДРпот(о), Вт. | |||||||||||
б = 2.
о. | 0,1. | 0,2. | 0,3. | 0,4. | 0,5. | 0,6. | 0,7. | 0,8. | 0,9. | 1,0. | |
е-2о. | 0,819. | 0,67. | 0,549. | 0,449. | 0,368. | 0,301. | 0,247. | 0,202. | 0,165. | 0,135. | |
Рпот(о), Вт. | |||||||||||
ДРпот(о), Вт. | |||||||||||
б = 3.
о. | 0,1. | 0,2. | 0,3. | 0,4. | 0,5. | 0,6. | 0,7. | 0,8. | 0,9. | 1,0. | |
е-3о. | 0,741. | 0,549. | 0,407. | 0,301. | 0,223. | 0,165. | 0,122. | 0,091. | 0,067. | 0,05. | |
Рпот(о), Вт. | |||||||||||
ДРпот(о), Вт. | |||||||||||
Проведем аналогичный расчет, выделив один итерационный этап для функции Ш (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. | ||
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. | ||
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.