Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

Зиганшин Г.З.

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

Формализация физических и различных технологических процессов очень часто сводится к построению дифференциальных уравнений, описывающих изменение какого-либо, например, выходного параметра (параметров) процесса при изменении других (воздействующих на него) параметров. При этом первичные уравнения подвергаются таким значительным преобразованиям с целью приведения их к линейным и обыкновенным, что полученные дифференциальные уравнения не удовлетворяют требованиям, предъявляемым к описывающим функциям управляемых объектов. Они описывают изменение выходного параметра в интервале времени только от момента окончания изменения и установления постоянных значений воздействующих (возмущающих и управляющих) параметров до момента достижения постоянного значения управляемого (зависимого) параметра. В них переменными являются отношения малых изменений параметров объекта к их принимаемым постоянным и начальным значениям, отсутствуют функции чувствительности управляемого параметра по воздействующим и скорости изменения воздействующих параметров. Такой подход возможен и рационален в тех случаях, когда τ2 >> τ1 где τ2 = t2 - t01 = t1 - t0, t0, t1 - моменты начала и окончания изменений воздействующих параметров; t0 , t2 - моменты начала и окончания изменения управляемого параметра соответственно. Однако реальные технологические процессы в промышленном производстве характеризуются непрерывным изменением воздействующих параметров и состояния. Как динамические они представляют собой бесконечную последовательность переходных процессов, при которых нет постоянных начальных значений. Каждое предыдущее значение параметров и состояния является начальным для последующего их изменения. При этом продолжительности исходных процессов соизмеримы (τ2 → τ1), управляемый параметр следует рассматривать как сложную функцию от функций, зависящих от одного независимого аргумента - времени.

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

з

Рис. 1. Графики изменения выходного параметра x(t) (1, 3) при изменении воздействующего u(t) (2, 4).

Для непрерывного процесса смешения двух потоков u и f [6- 1×10- 4 м3/сек] с плотностями ρ и ν [106 г/м3] одного и того же вещества уравнение баланса записывается: для t = t0 в виде

                                   uj0 - ρξ0 + fi0 × νη0 = (uj0 + fi0) xk0 ,                                   (1)

для t = t2 , где t2 - момент установления равновесия в системе после изменения воздействующих параметров, -

(uj0+u0β)(ρξ0+ρ0b)+ (fi0 + f0α)(γη0 + γ0s) =

                                                               = [(uj0 + u0β) + (fi0+ f0α)](xk0 + x0γ),                                      (2)

где uj0, ρξ0, fi0, γη0, xk0, j-ое, i-ое, η-ое, k-ое и ξ-ое начальные значения переменных; u0β, ρ0p, ƒ0α, ν0s, x0γ - β-ое, p-ое, α-ое, s-ое, γ-ое изменения переменных.  Подставляя (1) в (2), получим для  tt2

f               (3)            

Левая часть выражения (3) отражает потенциальное изменение управляемого параметра при мгновенном или опережающем изменении воздействующих параметров. Его можно назвать уравнением потенциального изменения, то есть уравнение (1) написано для tt0, а уравнение (2) - для tt2. На основании уравнения (3) и учитывая, что x0g(t) = x0g[ u(t), f0a(t), ρ0p(t), ν0s(t), t] для t0 < t< t2 можно записать                                                                 

               f

f                    (4)

Уравнение (4) отражает динамический баланс обобщенных сил воздействия и сопротивления, изменяющихся от начальных до новых установившихся значений. В тот момент, когда при  t = t1. изменение воздействующих параметров заканчивается, они принимают постоянные значения,

f

и из (4) для tt1, получаем x

f.  (5)

Для случаев постоянства концентраций в потоках ρξp = ρξ0, νηs = νη0 (4) и (5) можно записать в виде

f               (6)

f         (7)

 

В общем случае в (7)

x0g = x0g [ u(t), f0a(t), ρ0p(t), ν0s(t), t].            (8)  

Дифференцируя (8), получим

                  f

=f    (9)

Подставляя (9) в левую часть (7), получаем уравнение (4). Это означает, что обыкновенные дифференциальные уравнения (5) и (7) являются частными случаями дифферен-циального уравнения (4), представляющего обобщенное уравнение для описания динамики многих физических и технологических процессов в промежутке времени от момента начала изменения воздействующих параметров до момента установления постоянного значения управляемого параметра. Уравнения (4) и (6) удовлетворяют и условиям непрерывности изменения воздействующих и управляемых параметров.


Как видно, уравнение (4) получается путем замены левой части уравнения (3)

f     выражением

f

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

f

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

f (10)

f     (11)

Существование и единственность решения уравнения (4)определяется его решением при известных значениях функций чувствительности. Для упрощения  процедуры рассмотрим (6), отличающееся от (4) условиями pξpp= pξp, νηs=νη0..  Для этого обозначим правую часть (6)

f ff.

Тогда из (6) получим

f,

а при  f   -

f . (12)       

Обозначая f, из уравнения (12)

f.

Согласно [2] решение этого уравнения записывается в виде

f

где  f   и окончательно при xk0 = 0

f,

Принимая для примера  u = (1 - e- 0,8t) , f  ,

f,

f

Подставляя данные таблицы 5, где для одного случая γ = 1  F0γ = 0,0477, b0β = 0,4, a0a = 0,27, получаем x0γ = 0,0477 (1 - e- t) - 0,52 et (e0,2t - 1).

Определяя значения x01 при γ = 2 в области t ∈ [0;15], получим график 1 изменения x01 во времени (рис.1). График 2 характеризует изменение u(t). Отрицательная часть графика 1 обусловлена не процессом в системе, а только свойством аппроксимирующего решения, как например, полином, являясь универсальным средством приближения графиков, не может охватить абсолютно всю область значений функций.

Для случая  u = - (1 - e- 0,8t)  ,  f  получим

j(t) = F + (b0β - a0a) 0,8e- 0,8t

и решение для γ = 2

x02 = - 0,0477(1 - e t) + 0,52e- t (e 0,2t - 1).

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

Уравнение (1) применимо для описания баланса в одной любой зоне многозонного аппарата удаления технологической жидкости путем ступенчатого смешения (вымывания) ее при орошении твердых веществ менее концентрированной жидкостью. Статический баланс в четырехзонном аппарате, в котором потоки твердых веществ и орошающей жидкости двигаются навстречу (рис.2), для t = t0, описываются системой уравнений:

f

f

f

   f                     (13)  

а для  t0 < t< t2 - системой уравнений

f

f

f

f     (14)

Вычитая решение системы (13) относительно начального значения плотности на выходе любой зоны xδk0 из решения (14) относительно значений плотности xδkg на выходе той же зоны после установления равновесия в аппарате получим xδkg - xδk0 = xδ0γ, где δ - номер зоны аппарата. Тогда для плотности 4-ой зоны при возможных изменениях подачи воды в 4-ую зону, поступления сорбированной смеси в 1-ую зону и  плотностей поступающих в 4-ую зону орошающей и сорбированной смесей имеем

f       (15)  


где ρ40p - изменение независимого параметра орошающего потока; u40β -  изменение орошающего потока с целью управления процессом; f10a, ν10s - изменения параметров независимого потока, являющегося нагрузкой на аппарат. Здесь приняты изменения f10a и ν10s с индексами первой зоны, относящимися к первой зоне потому, что именно они определяют плотности на выходе всех зон. Точно также в уравнение первой зоны должны войти изменения u40b, ρ40p параметров потока четвертой зоны. Для любой зоны уравнение (15) записывается, опуская номера зон при u40b, f10a, ν10s, ρ40p, в виде         

    p


Рис. 2. Схема технологических потоков.

f     (16)          

Рассмотрим (16) для сокращения процедур вычисления при ρξp = ρξ0 = 1,  νηs = νη0 =1,8 и получим

f  ,     (17)

где

f              (18)          

Форма записи (17) обусловлена необходимостью указать на определяемость значения функции Fδ из системы (18) и невозможностью явного выражения ее из системы (18). Значение функции Fδ  из (18) определяется только путем вычисления.

Для проверки существования решения уравнения (17) (множества решений для каждого γ = 1, ..., t) необходимо решить две задачи: определение значений  Fδ  и определение значений функций чувствительности

f, f.

Определим значения  xδ0γ f0α. Поскольку значения плотностей смесей на выходе всех зон при различных значениях подачи воды и расхода смеси в 1-ую зону определяются из системы четырех уравнений (14), то решение задачи сводится к дифференцированию этой системы:

f             (19)       

f      (20)         

Решение задач (19) и (20) встречает значительные трудности, нужны другие методы.

Для решения задачи проведена следующая работа. По уравнению (14) вычислены равновесные значения плотностей по зонам для трех значений подачи воды u  и трех значений расхода смеси fia.

Результаты вычислений приведены в таблице 1. Для наглядности по табл. 1 построены 6 графиков (рис. 3). Данные таблицы аппроксимированы уравнениями, приведенными в табл. 2.. Дифференцируя по воздействующим параметрам уравнения табл. 2, получим дифференциальные уравнения в частных производных (табл. 3), представляющих собой функции чувствительности. При всех вычислениях значения расходов уменьшены в 100 раз для повышения точности вычислений. В таблицах приведены уменьшенные их значения.

Далее, продолжая подготовку данных к проверке уравнения (17), необходимо получить значения ∂xδ / ∂u4, ∂xδ / ∂f1. Для этого решаются уравнения, приведенные в табл. 3, подставляя в них значения u4 и f1. Во избежание громоздких вычислений выполним их только для плотности x4, т.е. для первых шести уравнений табл. 3. Результаты вычислений приведены в табл. 4.

Таблица 1.Значения плотности x4 при различных значениях  u4, f1

№№ кривых

1

2

3

№№ кривых

4

5

6

u4

f1 =1,08

f1 = 0,90

f1 = 0,72

f1

u4 = 1,08

u4 = 0,90

u4 = 0,72

1,08

1,16

1,1075

1,0606

1,08

1,16

1,2229

1,3071

1,02

1,1786

1,1226

1,0708

1,02

1,1422

1,2023

1,2853

0,96

1,1977

1,14

1,0870

0,96

1,1247

1,1813

1,2545

0,90

1,2229

1,16

1,0974

0,90

1,1075

1,16

1,238

0,84

1,2485

1,2828

1,1148

0,84

1,0909

1,1387

1,2127

0,78

1,2765

1,2087

1,1354

0,78

1,0752

1,1177

1,1866

0,72

1,3071

1,2380

1,16

0,72

1,0606

1,0975

1,16

Таблица 2 Уравнения плотностей при равновесных состояниях в круговом аппарате

№№ уравнений

Уравнения

1

x4 = 1,8727 - 1,036 u4 + 0,3482 f,

2

x4 = 1,804 - 1,07 u4 + 0,3938 f,

3

x4 = 1,6648 - 0,988 u4+0,3968 f,

4

x4 = 0,928 + 0,1214 f1+0,0866 f

5

x4 = 0,8503+0,338 f1+0,0066 f

6

x4 = 0,7813+0,605 f1-0,1095 f

Таблица 3. Дифференциальные уравнения по воздействующим параметрам

№№ уравнений

Уравнения

1

x4 / ∂u4 =0,69634 u4 - 1,036

2

x4 / ∂u4 =0,7876 u4 - 1,07

3

x4 / ∂u4 =0,7936 u4 - 0,988

4

x4 / ∂f1 =0,17326 f1 + 0,1214

5

x4 / ∂f1 =0,0134 f1 + 0,338

6

x4 / ∂f1 =0,219 f1 + 0,605

p

Рис. 3. Графики зависимости x4 от подачи воды (1 - 3) и расхода смеси (4 - 6)

Для определения значений F необходимо решить систему (14) при одновременном изменении расхода смеси f1 и подачи воды u4. Эти данные можно получить из табл. 2 следующим образом. Первые строки данных для кривых 1, 4 в табл. 2 приняты за начальные значения. Полагая изменения u и f0a на 0,06, соответствующее изменение x40γ получим вычитанием первой строки из второй [x40γ (1,08) - x40γ (1,02)].

Таблица 4. Значения функций чувствительности x40γ по  u40β и f10a.

№№ уравнение

1

2

3

№№ уравнений

4

5

6

u4

f1=1,08

f1=0,90

f1=0,72

f1

u4 =1,08

u4 =0,90

u4 =0,72

1,08

-0,28

-0,22

-0,1309

1,08

0,3084

0,3522

0,3685

1,02

-0,326

-0,2666

-0,1785

1,02

0,2980

0,3514

0,3816

0,96

-0,3678

-0,3139

-0,2261

0,96

0,2876

0,3505

0,3947

0,90

-0,41

-0,3611

-0,2737

0,90

0,2773

0,3499

0,4079

0,84

-0,451

-0,4085

-0,3214

0,84

0,2669

0,3491

0,4200

0,78

-0,4931

-0,4556

-0,369

0,78

0,2565

0,3483

0,4342

0,72

-0,535

-0,5029

-0,4166

0,72

0,2461

0,3475

0,4473

Далее, вычитанием первой строки из третьей при существующем изменении u40β и f10a на 0,06 и т.д. Результаты внесены в графы 3, 4, 6, 7 табл. 5. Значение F получается алгебраическим суммированием данных граф 4 и 7. Для кривых 2 и 5 (рис.3) за начальные значения fi0, ui0 приняты строки при x4 =1,16 (как в первом случае). Шаги изменения сделаны вверх и вниз от этой строки. Результаты вычислений тоже внесены в графы 4, 7. Значения F получаются также суммированием.

Для кривых 3, 6 за начальные значения приняты (при uj0 = 0,72) данные последних строк тоже при x4 = 1,16. Шаги изменения сделаны снизу вверх. Остальные приемы аналогичны. Остается указать, как пользоваться данными таблицы 5 для проверки уравнения (17). Для выяснения роли функций чувствительности в уравнении (17) необходимо взять данные кривых, исходящих из одной точки 1,16, сходящихся к ней или пересекающихся в той точке. Это кривые 1 и 4, 2 и 5, 3 и 6. Итак, берем кривые I и 4. Уравнение  (17) принимает вид       

f     (21)             

Таблица 5. Изменение F при одновременном изменении u  и f0a, x4k0 = 1,16

кри-

вых

fi0

ui0

 

u

x40γ

кри-

вых

f0a

x40γ

F

1

2

3

4

5

6

7

8

1

 

0,06

0,0188

4

-0,06

-0,0178

0,0010

1,08

0,12

0,0377

0,12

-0,0353

0,0024

1,08

0,18

0,0629

0,18

-0,0525

0,0104

 

0,24

0,0885

0,24

-0,0691

0,0194

 

0,30

0,1165

0,30

-0,0848

0,0317

 

0,36

0,1471

0,36

-0,0994

0,0477

2

 

0,18

-0,0525

5

0,18

0,629

0,0104

0,90

0,12

-0,0374

0,12

0,0423

0,0049

0,90

0,06

-0,02

0,06

0,0213

0,0013

 

-0,06

0,0228

-0,06

-0,0213

0,0015

 

-0,12

0,0487

-0,12

-0,0423

0,0064

 

-0,18

0,0780

-0,18

-0,0625

-0,0155

3

 

0,26

-0,0994

6

0,36

0,1471

0,0477

0,72

0,30

-0,0892

0,30

0,1253

0,0361

0,72

0,24

-0,0770

0,24

0,0995

0,0225

 

0,18

-0,0626

0,18

0,0780

0,0154

 

0,12

-0,0452

0,12

0,0527

0,0075

 

0,06

-0,0246

0,06

0,0266

0,0020

Из уравнения (21) следует, что при равенстве скоростей изменения u0β, f0α

f

условие равенства нулю левой части уравнения не выполняется, хотя по графикам оно должно быть выполнено. Дело в том, что значения ∂x40γ / ∂u, ∂x40γ / ∂f0a, взятые в табл. 1, относятся только к точке u = 1,08, f0a = 1,08, а изменение x40g - к концу интервала от 1,08 до 1,02. Если посмотреть значения функций чувствительности, то можно увидеть, что они увеличиваются по линейному закону, как и следовало ожидать по их уравнениям в табл. 3. В то же время по графикам, например 4 и 6, они должны быть постоянными. Здесь, по-видимому, проявляется вероятность получения бессмысленных результатов при дифференцированной экспериментальных данных [2]. Здесь это проявляется в том, что первичное уравнение недостаточно точно описывает график, близкий к прямой линии, которое было взято для описания уравнения прямой линии, чтобы сохранить подобие другим уравнениям. Поэтому для проверки существования решения уравнения (17) возьмем данные из графиков. Из них следует, что F0g (u40b, f10a) Крутизны графиков 1 и 4 можно считать равными и, как следствие, отношения  x40g к шагу изменения un0b = 1,08 - 1,02 = 0,06, f10a = 1,08 - 1,02 = 0,06

f

Выражение (17) принимает вид (опуская номера зон в индексах):

f

Запишем его в виде

f

По справочнику [3] решение уравнения (17) записывается

f

что приводит к выражению x = F (f1 - e- t), а для случая du / dt = df0a / dt при  t → ∞ x0γ = 0,001.

Проверим в случае неравенства функций чувствительности для графиков 3 и 6:

1. Для изменения df0a = du = 0,36, F = 0,0477

f

решение пишется в виде f

где             f

Для   du / dt = df0a / dt, u = 1 - e- 0,8t

f

2. Для изменения df0a = du = 0,30, F = - 0,0317

f

Для    f

f

и т.д. Результаты представлены в таблице 6.

Из данных графы 8 таблицы 6 следует, что результаты вычислений значений плотности на выходе четвертой зоны аппарата совпадают с данными графиков на рис.3. Проверка осуществлена для 24-х графиков. Такая высокая точность обусловлена тем, что здесь приняты равными скорость изменения подачи воды на орошение и скорость изменения поступления смеси в массе. Естественно, на практике эти скорости отличаются и уравнение (17) принимает вид:

f

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

Таблица 6. Значения плотности x вычисленные по уравнению  (22)

u0β

f

f0a

f

А=гр.2+

+гр.4

А f0a

F = гр.8

табл.5

δ = гр.6 - - г р.7

1

2

3

4

5

6

7

8

0,06

0,3133

0,06

-0,2967

0,0166

0,001

0,0010

0

0,12

0,3142

0,12

-0,2942

0,1200

0,0024

0,0024

0

0,18

-0,3494

0,18

-0,2917

0,0523

0,0094

0,0104

0,001

0,24

-0,3688

0,24

-0,2879

0,0809

0,0194

0,0194

0

0,30

-0,3883

0,30

-0,2827

0,1056

0,0317

0,0317

0

0,36

0,4036

0,36

-0,2761

0,1325

0,0477

0,0477

0

ВЫВОДЫ

  1. Реальные технологические процессы в промышленном производстве  характеризуются непрерывным изменением воздействующих параметров и состояния. Как динамические системы они представляют собой бесконечную последовательность переходных процессов, при которых нет постоянных начальных значений. Каждое предыдущее значение параметров и состояния является начальным для последующего их изменения, продолжительности переходных процессов соизмеримы, управляемый параметр следует рассматривать как сложную функцию от функций, зависящих от одного независимого аргумента - времени. Соотношение начальных значений и изменений параметров и состояния определяется двусторонним неравенством.
  2. Получено уравнение потенциальных изменений на примере процесса смешения двух потоков, на основании которого построено обобщенное дифференциальное уравнение в частных производных, описывающее динамику многопараметрического процесса с учетом чувствительности выходного параметра к изменению воздействующих и скоростей их изменений. Осуществлена проверка существования решения полученного уравнения и показано множество решений, зависящих от множества воздействующих параметров и множества их изменений.
  3. Достоверность использованных при построении обобщенного дифференциального уравнения положений и реальность собственно дифференциального уравнения доказаны путем преобразований его в классическое однопараметрическое линейное дифференциальное уравнение при уменьшении отношения продолжительности изменения воздействующих параметров к продолжительности изменения выходного параметра, уменьшении числа воздействующих параметров до единицы, уменьшении отношения скорости изменения выходного параметра к скоростям изменений воздействующих параметров.
  4. Приведены результаты использования построенного дифферен циального уравнения в частных производных для описания многоступенчатого многопараметрического процесса вымывания технологической жидкости при встречно-пересекающихся потоках.

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

  1. Воронков И.М. Курс теоретической механики. - М.: Наука, 1965, 596с.
  2. Шуп Т. Решение инженерных задач на ЭВМ. - М.: Мир, 1982, 295с.
  3. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. - М.: Наука, 1971, 576с.