Управление температурными полями объектов с распределёнными параметрами Текст научной статьи по специальности «Математика»

УДК 681.5

УПРАВЛЕНИЕ ТЕМПЕРАТУРНЫМИ ПОЛЯМИ ОБЪЕКТОВ С РАСПРЕДЕЛЁННЫМИ ПАРАМЕТРАМИ

А.Б. Чернышев

ГОУ ВПО «Северо-Кавказский государственный технический университет», филиал в г. Кисловодске E-mail: chalbor@rambler.ru

Проведено исследование нелинейной системы управления температурным полем. Рассмотрен процесс формирования температурного поля под воздействием тепловых источников, включение которых осуществляется с помощью релейных элементов. Проводится анализ влияния источников в зависимости от их координат и времени включения. Приведена оценка отклонения температуры от заданного значения. Предложена методика оценки шага дискретизации распределённого управляющего воздействия. Исследование проводится на основе метода функций Грина.

Ключевые слова:

Температурное поле, управляющие воздействия, секционный нагреватель, шаг дискретизации, функция Грина.

Многие технологические процессы требуют обеспечения управления температурными полями, поэтому задача разработки методик синтеза регуляторов для систем управления объектами с распределёнными параметрами становится всё более актуальной. Одним из методов описания распределённых объектов является способ на основе импульсных переходных функций или функций Грина [1-3]. В качестве пространственного объекта рассмотрим тонкий однородный цилиндрический стержень длины l, теплоизолированный с боков. На концах стержня поддерживается нулевая температура, в начальный момент времени температура во всех точках равна нулю. Будем полагать, что управляющим воздействием является тепловой поток, создаваемый источниками, равномерно распределёнными по границе боковой поверхности цилиндра. Включение источников осуществляется с помощью релейных элементов. Значение температуры в точке х в момент времени t, вызванной действием мгновенного точечного источника, действующего в момент тв точке C, определяется с помощью функции Грина [1, 2]. Функция Грина для отрезка может быть представлена в виде ряда Фурье [4].

G (х, t ,£,т) =

^ да

= i ¿ exp

1 n=1

Поставим задачу стабилизации температуры на уровне некоторого значения Тзад. В начальный момент времени т0=0 включаются все источники. Формируется функция «начального нагрева» (рис. 1). Предполагается, что мощность каждого источника равна единице, в этой связи размерности правых и левых частей некоторых уравнений формально могут не совпадать. Каждый из единичных мгновенных точечных источников оказывает воздействие на все точки стержня. Значение температуры в произвольной фиксированной точке x¡ может быть представлено в виде соотношения:

Т (Xj, t,T0) =

nnaY . nn . nn „

----- I t sin---Xj sin Ci,

l J l 1 l

(t -т)

. nn . nn

sin-------X sin-------<

l l

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

При достижении в какой-то точке ху значения Тзад, в момент времени т включается источник, соответствующий датчику с номером у, и оказывает влияние на все точки (рис. 3). При этом продолжается действие всех источников, включённых при значении т0=0.

Таким образом, получим общий вид выходной функции:

Т (х,, г ,т) =

d k л

=хх 2 “p

/=1 n=1 1

l

. nn . nn „ sin----X,. sin-----------C +

l 1 l '

+¿ ¿2 exp

^ n=1 ^

nna Y ,

—J(t-т)

+exp

exp

na

t

П ^ . П

sin—X > sin —C, +

l tí r

2na

l

+exp

3na

l

. 3n ^ . 3n sin—x> sin—+...

l ^ '

l

При чётных значениях п положительные и отрицательные составляющие суммы

. пп

151П “Г4 /

/=1 Г

взаимно исключают друг друга. При нечётных значениях п

^ . пп . пп „ 1 П „

2-™ Т4^511' 7^'= п Т4'

/=1

Тогда:

Т (Xt) = у ¿ ^^exp

l n=1 2n 1

. (2n -1)n

x sin--------------x

l

(2 n -1)na l

(1)

При увеличении времени t каждое из слагаемых формулы (1) будет убывать. Влияние на итоговую

/ \ 2 _

I па

сумму второго слагаемого будет в 3 exp

na

~Т

t

где р=1,2,3,... - порядковый номер включения источника, 1(р) - один из источников, тр - момент времени включения источника под номером 1(р). Или

Т(х,,г, т) = Т(х,,(,т0) + ^Т(х,,г-тр ).

Р

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

раз меньше первого. С заданной точностью е можно определить значение I *, начиная с которого слагаемым с номером п можно пренебречь.

* 1 ( I } \ 1

г >—2— I — I 1п—

п -1 \па) пе

Таким образом, в момент времени

г > -( —I 1п— функция, определяющая распре-8 \па ) 3е

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

ч2 " г

Т (х, t) = — exp

na

l

JL ' . JL ^

sin—x> sin—c.

l t1 r

(2)

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

т1 =| — I ln

. na

( _ d п

2sin — х ¿ sin—C l Kpj=i l '

lT„

при выполнении условия:

~ • n 2sin—xm

l 1

(3)

(4)

Включение источников осуществляется с помощью релейных элементов, поэтому т1 - время первого переключения реле, Тза() - пороговое значение. В момент времени т1 источники, соответствующие крайним датчикам, вырабатывают сигнал Т (х, г-т1) =

2 V

= т ¿ exp

n=1

i na , , ч

-\ — I (t-T1)

(5)

Преимущество полученной формулы в том, что действие источников выражается с помощью только первой гармоники ряда Фурье. Величина

0 d _

2 X-1 • П е

sin—ti является постоянной относительно x и

1 j=1 I

t. Её значение зависит только от количества источников. При значениях t=0, к=3 и

0 d -г

2 • П ^

—^ sin —t, t = 1 формула (1) примет вид:

1 i=l 1

T (х,0) = sin—X + isln3^ x + isln5^ x.

l 3 l 5 l

Этот сигнал добавляется к имеющемуся значению Т(х,/). В результате наложения сигналов, при значении От1 для середины отрезка, получим:

т (х, г,т1) = т (х, г)+2Т (х, г-т1).

Импульс, приложенный в момент т1 в точке 4, с течением времени От1 оказывает воздействие на все точки отрезка. Время «прихода» максимума этого сигнала для всех точек различно. Это время зависит от расстояния между точкой включения источника и точкой наблюдения, а также от близости источника к концу отрезка. От длины отрезка зависит значение максимальной температуры, «дошедшей» от источника до точки наблюдения. В связи с тем, что на концах стержня во время всего

процесса температура поддерживается равной нулю, максимальная температура с течением времени будет стремиться к середине стержня. В качестве точки наблюдения будем рассматривать точку x=l/2. Выражение функции источника для неограниченной прямой имеет вид [4]:

G(x,4,t-т) =—, 1 -exp (x 4)

I—2--------

2^па (t-т)

4а (t -т)

Используя это выражение, получены формулы, выражающие время достижения максимального значения температуры в середине отрезка при воздействии источника 4 е ( 0;Г).

2а

- + Т

при

10 4) -+

2а

10

-4

4а2

■ + т1; при

— <4 < -10 1 2

— <4 < —.

10 10

+Т; при 0<4 < 10 (6)

В момент «прихода» максимума сигнала значение температуры в середине отрезка будет определяться формулой

Т (2-т )=Т (2,)+2Т (2-).

Отклонение температуры от заданного значения можно определить из равенства

¿1 =

T\{. tm Т |-Тшд

Полученная погрешность характеризует качество релейного управления. В выражении значение Т(1/2,т) вычисляется по формуле (2), значение т1 -по формуле (3), при выполнении условия (4), Т(//2,/м-г1) - по формуле (5), , - из соотношений (6). Таким образом, получено:

\2

¿1 =

4 V

+-X ехР

l п=1

где

2

- exp

ппа

I

(tm -Т1)

. пп . пп

~~2'$1 - Тзад

ln

(2 . п ^ . п4^

2sin ~x1 X sin T4

iTa

При условии

251п пгх^Р Ё *1пП4> >1Т^.

Г ¡=1 Г

Время второго включения источника (переключения реле) определяется по формуле:

2т,

-+т.

I--

2T I

jL1 зад1

п sin

I 1

Тогда т2-т1 - временной интервал переключения реле.

Результаты компьютерного моделирования описанного процесса согласуются с выводами и зависимостями, полученными с использованием математического аппарата.

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

• Задать значения: Д1 - допустимая погрешность, Тзад - требуемая температура, / - длина стержня, а — коэффициент температуропроводности материала.

• Изменять значение параметра п - количество секций нагревателя.

• Для каждого значения п вычислять: 4, Т, т2, ,

$1 = 2-; т1 = 2п

па

ln

( п d п ^ 2sin — x1 X sin—4 i -

ITad

2т1

I —

2T l +Т1’

• п

п sin—4

tm =

2 4) 3l < l

+ Т1; при 10 <4 < 2

10 4) f l < 31

_____+_+т1; при _<4<_

10-4) l' 0< 4 < l

_____+_+т1; при 0<4<_

Проверять выполнение условий:

tm =

T1 =

A1 + Tsaö---------------exp

• n e

siny £ =----------

l 4

l

na

exp

na

~T

(tm -T1 )

4 . п

-вШ-X! > Тзад. п I

При выполнении указанных условий зафиксировать значение п и определить размер секции: I

Г =- •

п

Проверять выполнение условий:

\2

. п „ l

sin—g1 =— exp l n

i na .

-l — I (2T T2)

(

1 - exp

na

(T2 -T1)

\

4 . n

-sm- X1 > Тзад. n l

• При выполнении указанных условий зафиксировать значение п и определить размер секции:

I

Г2 =~ • п

• Из полученных значений г1 и г2 выбрать меньшее:

г = ш1п(г1; г2}.

Если г!>г2, т. е. г=г2, то желаемый результат достигнут. Колебания, в пределах допустимой амплитуды, будут совершаться относительно заданного значения выходной функции.

Если г:<г2, т. е. г=г1, то амплитуда колебаний будет в пределах допустимой, но функция, совершая колебательные движения, будет медленно убывать и через определённое время выйдет за пределы допустимого интервала относительно заданного значения.

СПИСОК ЛИТЕРАТУРЫ

1. Рапопорт Э.Я. Структурное моделирование объектов и систем управления с распределёнными параметрами. - М.: Высшая школа, 2003. - 299 с.

2. Карташов Э.М. Аналитические методы в теории теплопроводности твёрдых тел. - М.: Высшая школа, 2001. - 550 с.

3. Першин И.М. Анализ и синтез систем с распределенными параметрами. - Пятигорск: РИА КМВ, 2007. - 244 с.

4. Тихонов А.Н., Самарский А.А. Уравнения математической физики. - М.: Наука, 1972. - 736 с.

Поступила 06.02.2009 г.

УДК 621.311

АНАЛИЗ ПЕРЕМЕННЫХ РЕЖИМОВ ВЕРТИКАЛЬНОГО СМЕШИВАЮЩЕГО ПОДОГРЕВАТЕЛЯ

Н.Н. Галашов

Томский политехнический университет Е^У: gal@tpu.ru

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

Ключевые слова:

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

Подогреватели смешивающего типа нашли широкое применение на тепловых электростанциях, но до сих пор их процессы еще слабо изучены. Наибольшее внимание было уделено вопросам надежности работы подогревателей в переходных режимах при сбросах нагрузки [1, 2]. Вопросы теплообмена в смешивающих подогревателях исследованы в узком диапазоне давлений, всего до 0,13 МПа [3], в то время как существующие подогреватели работают при давлениях до 0,2 МПа, и предполага-

ется использовать их при более высоких давлениях. Имеются единичные исследования по изменению параметров в этих подогревателях при изменении режимов работы турбоустановки [4], которые не позволяют оценить влияние отдельных факторов на характеристики отсеков подогревателя, что требуется при проектировании подогревателей.

В данной работе рассматриваются вопросы изменения по греющим отсекам подогревателя нагрева воды Д/ и ее недогрева до температуры насы-