Построение математической модели процесса охлаждения потока движущей среды в пространстве состояний

An application of the spectral theory for construction of mathematical model for cooling process of moving medium flow in state space is considered. The process of flow cooling is described by a first-order hyperbolic equation. Comparison results of the realized model transient and the reference model transient received by an operator method are given.

Keywords: cooling of flow, spectral theory, state space, mathematical model, vector-matrix Cauchy form, vector of a spectral characteristic.

Основная часть

Теплообменные аппараты различных конструкций широко применяются в различных отраслях промышленности. В качестве базовой математической модели тепловых процессов, протекающих внутри теплообменного аппарата перекрестного тока, может использоваться уравнение теплового баланса, согласно которому распределение температуры охлаждаемого потока по длине трубки теплообменника в зависимости от координаты точки z и времени t описывается уравнением в частных производных [1]

; , (1)

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

,. (2)

Здесь V - скорость потока; - функция изменения коэффициента теплообмена; - температура среды, принятая одинаковой по всей длине теплообменника; L - общая длина трубки теплообменника; - начальное распределение температуры; - функция изменения температуры потока на входе теплообменника.

Управление температурой потока осуществляется за счет изменения скорости течения охлаждающего потока, что учтено в модели (1)-(2) коэффициентом . Для синтеза системы автоматического управления температурой потока на выходе теплообменного аппарата необходимо получить передаточную функцию по каналу «коэффициент теплообмена» - «температура на выходе», - . Аналитического решения уравнения (1) не существует, однако поведение температурного распределения потока можно рассмотреть в отклонениях от установившегося режима [2, 3].

Представим уравнение (1) в следующем виде:

. (3)

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

, , (4)

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

(5)

.

В установившемся режиме при и будем иметь

. (6)

После решения уравнения (6) получим следующее:

. (7)

Подставив (7) в (5) и упростив, получим

(8)

или, введя обозначение , получим

. (9)

Уравнение (9) представляет собой математическую модель теплообменного аппарата в отклонениях от установившегося режима, определяемого параметрами , , , .

Модель (9) можно преобразовать к бесконечномерной системе обыкновенных дифференциальных уравнений в форме Коши, используя свойства спектральных характеристик [4]. Для этого каждую из функций, зависящую от пространственной координаты , и следует представить в виде обобщенного ряда Фурье по пространственной переменной на основе выбранной ортонормированной системы функций , :

, , ,

где - бесконечномерная квадратная операционная матрица дифференцирования первого порядка [4],

, .

На основании свойств спектральных характеристик уравнение (9) в спектральной форме по пространственной переменной примет вид

, (10)

где - бесконечномерная единичная матрица.

Представление объекта (10) в векторно-матричной форме Коши имеет вид

(11)

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

, . (12)

- матрица, состоящая из нулей; - матрица, составленная в соответствии с [4] из элементов ортонормированной системы разложения, вычисленных для фиксированных значений пространственных координат на интервале .

В качестве системы разложения выберем систему функций

; . (13)

Для выполнения вычислительных процедур ограничимся первыми пятью членами разложения в системе (13) и соответственно числом уравнений в системе (11). Для параметров модели м/с, 1/с, , м матрицы ,, и будут иметь следующий вид:

; ; (14)

; .

Вектор управляющих воздействий принимает вид

.

Для оценки качества модели, полученной на основе спектрального метода, сравним ее переходный процесс (см. рисунок) с переходным процессом эталонной модели, полученной на основе операторного метода,

.

В момент времени t=0 происходит скачкообразное изменение управляющего сигнала на входе исследуемой и эталонной моделей с нуля на -0.01.

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

Анализ переходных процессов показал достаточную для задач синтеза систем управления точность исследуемой модели. Представление модели процесса охлаждения потока в пространстве состояний (10) в дальнейшем будет использоваться для синтеза систем автоматического управления режимами работ теплообменных аппаратов.

Библиографический список

1. Данилушкин И.А., Россеев Н.Н. Синтез системы автоматического управления температурным полем трубчатого теплообменника. Самара: СамГТУ.

2. Данилушкин А.И., Рапопорт Э.Я. Алгоритмы функционирования процесса непрерывно-последовательного индукционного нагрева // Алгоритмизация и автоматизация технологических процессов и промышленных установок: Межвуз. сб. научн. тр. Вып. VII. Куйбышев: КПтИ, 1976. С. 118-124.

3. Россеев Н.Н., Данилушкин И.А., Кузнецов П.К. Модель распределения температуры масла в аппарате воздушного охлаждения // Математическое моделирование и краевые задачи: Труды III Всеросс. научн. конференции. Ч. 2. Моделирование и оптимизация динамических систем и систем с распределенными параметрами. Самара: СамГТУ, 2006. С. 142-144.

4. Коваль В.А. Спектральный метод анализа и синтеза распределенных управляемых систем. Саратов: Изд-во Сарат. гос. техн. унт-та, 1997.

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