CC BY
26
УДК 519.711.3
МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ ПОТОКА ЖИДКОСТИ ПРИ ПАРАМЕТРИЧЕСКИХ ВОЗМУЩЕНИЯХ
©2010 И. А. Данилушкин Самарский государственный технический университет
Разработана методика моделирования параметрических возмущений при исследовании температурного поля потока жидкости как объекта управления с распределёнными параметрами. В качестве модели объекта используется уравнение первого порядка в частных производных. Полученные результаты позволяют реализовать ступенчатое изменение параметров объекта в компьютерных пакетах численного моделирования динамических систем.
Объект с распределёнными параметрами, параметрическое возмущение, моделирование
Учёт пространственной
распределённости объектов управления при синтезе систем автоматического управления позволяет повысить точность моделей, а, в ряде случаев, является единственно возможным способом получения адекватной модели. Процесс нагрева потока жидкости за счёт теплообмена со стенкой трубки может быть описан дифференциальным уравнением в частных производных первого порядка:
ЩМ+„(0.ЩМ=р (,Нв1 (0 - 0( х,,))
аt ах
в(0, t) = Е0), в(х,0) = во(х) , 0 < х < Ь, (1) где в(х, 0 - температурное поле
жидкости, находящейся в трубке длиной Ь; вНО - температура стенки трубки, которая принимается одинаковой по всей длине; у(^) - скорость потока,
изменяющаяся во времени. Приведённый коэффициент теплообмена Р ^),
зависящий от радиуса трубки, физических свойств жидкости и коэффициента конвективного теплообмена между стенкой трубки и жидкостью также может изменяться во времени.
Решение уравнения (1) для у(^ и Р ^) не может быть получено аналитически. Возможный подход к получению передаточных функций по каналам скорость потока - температура потока и приведённый коэффициент
теплообмена - температура потока
предполагает рассмотрение процесса в отклонениях от установившегося режима
[1] с последующей линеаризацией
дифференциального уравнения и исключения из рассмотрения
произведения приращений, как
бесконечно малых высшего порядка. Подобный подход не обеспечивает параметрических изменений модели, а лишь привносит аддитивное возмущение на выход объекта [1 ].
Для оценки влияния
параметрического возмущения в виде
ступенчатого воздействия на
температурное распределение потока
предлагается в качестве начального распределения в0(х) в (1) принять температурное распределение,
соответствующее установившемуся
режиму для фиксированных значений параметров у(0=у0, Р ^) = Р0, вт = втс ■ Тогда, решение уравнения (1) при других значениях параметров у(^=у1 и/или Р ^) = Р1 будет описывать реакцию объекта на соответствующее
параметрическое возмущение.
Рассмотрим математическую
модель объекта, описываемого
следующим уравнением:
+V. ММ+Рв(хА = Рвт(t). (2)
дt ах
Решение неоднородного уравнения
(2) с учётом начальных и граничных
Управление, вычислительная техника и информатика
условий может быть найдено методами структурной теории распределённых систем [2, 3], как результат пространственного интегрирования
передаточной функции распределённого объекта:
Ж(х, £, р) = 1(х - X) • - • ех/- — + Ь
и стандартизирующей функции
х,р) = /З О,т(р) + °>0(х) + V• 8(х)• g(р), (4) где р - оператор преобразования Лапласа.
В установившемся режиме влияние начального температурного
распределения отсутствует. Поэтому,
приняв О0(х) = 0, получим
О( х р) = | Ж(x, £, р) • а(Х, р) • =
p
p + p
l - exp
p+p
+
(5)
//
+ g( p) exp
p+p
x
Из (5) можно получить температурное распределение потока, соответствующее установившемуся
режиму, при QT (p) = 1 • QTC , QTC ° const,
gC ° const. Согласно
g(p) = — • gc, p
предельной теореме
Q¥ (x) = Q( x ¥) = Hm[p •Q( x, p)]:
p®G
От
Ґ / l - exp
V
-1 x
v
\\ / + gC exp
V
//
-1 x
v
V
(б)
/
Приняв (6) начальным
температурным распределением потока при скорости V = ^, можно записать выражение для температурного распределения потока при скорости V = V!. Для этого необходимо выполнить процедуру пространственного
интегрирования произведения
передаточной функции (3) с выражением (6). Полученный результат позволит записать выражение для температурного распределения потока при ступенчатом изменении скорости:
p
p+p
1- exp
/ О W
p+p x
-— x
V
v
+
//
(3) + Q
+ g( p) exp l
p + p
x
+
тс , о p+p + (gc - °тс ^
l - exp
p+p
x
+
//
(В)
vG p + p (vG - vl )
x
x
exp
V
vG
- exp
/
/ n \
p + p x
-——x vl
Аналогично может быть получено выражение для температурного распределения потока при ступенчатом изменении приведённого коэффициента теплообмена Р с Р0 на Р1:
Q(x p) = Q (p>
pl
p + pl
l - exp
p + pl
W
x
+
//
+ g( p) exp
+ °тс “p
p+pl
+ (gc - °тс ^
x
+
l - exp l
p + pl
x
+
//
(9)
p + pl - pG
x
Г pg 1 г
exp x - exp
_ v V у V
p + pl
Выражения (8) и (9) могут быть реализованы в компьютерных пакетах численного моделирования
сосредоточенных динамических систем при моделировании поведения
температуры в фиксированной точке х.
Библиографический список
1. Данилушкин А.И., Рапопорт Э.Я. Алгоритмы функционирования процесса непрерывно-последовательного индукционного нагрева// Алгоритмизация и автоматизация технологических процессов и промышленных установок. Межвузовский сборник научных трудов. Вып. VII. - Куйбышев: КПтИ, 1976. -С.118-124.
v
v
v
x
v
v
v
v
v
x
v
2. Бутковский А.Г. Структурная теория распределённых систем. М.: Наука, l9ll.
3. Рапопорт Э.Я. Анализ и синтез
систем автоматического управления с распределёнными параметрами: Учеб.
пособие. М.: Высш. шк., 2GG5.
References
1. Danilushkin A.I, Rapoport E.Ya. Algorithms for the operation of the process of continuous-sequential induction heating// Algorithmization and automation of
technological processes and industrial plants. Interuniversity collection of scientific papers. - Vol. VII. - Kuybyshev: KPtI, 1976. - pp. 118-124.
2. Butkovskiy A.G. Structural theory of distributed systems. - Moscow: Nauka, 1977.
3. Rapoport E.Ya. Analysis and synthesis of automatic control systems with distributed parameters: tutorial.- Moscow: Vysshaya Shkola, 2005.
SIMULATION OF THERMAL FIELD OF FLUID FLOW IN THE PRESENCE OF
PARAMETRIC DISTURBANCES
©2010 I.A. Danilushkin
Samara State Technical University
The method of simulation of parametric disturbances during the research of thermal field of fluid flow as a plant with distributed parameters is developed. A first-order partial differential equation is used as a plant model. The obtained results allow simulating a stepped variation of plant parameters in program packages of computational modeling of dynamic systems.
Plant with distributed parameters, parametric disturbance, simulation
Информация об авторе
Данилушкин Иван Александрович, к.т.н., доцент, кафедра автоматики и управления в технических системах Самарского государственного технического университета, [email protected]. Область научных интересов: объекты и системы с распределёнными параметрами.
Danilushkin Ivan, PhD in Technical Sciences, associate professor, automation and control in technical systems department of Samara State Technical University, [email protected]. Research interests: plants and systems with distributed parameters.
