Автоматизация научных исследований
морей и океанов
УДК 551.46:681.51
Л.А. Краснодубец*, В.И. Забурдаев**, В.В. Альчаков*
Управление морскими буями-профилемерами как метод повышения репрезентативности термохалинных измерений. Синтез и моделирование системы управления
Представленная работа является второй частью статьи тех же авторов, опубликованной в «Морском гидрофизическом журнале» в № 4 за 2012 г. В ней выполнен анализ нелинейных дифференциальных уравнений, описывающих процесс управления гидрофизическим зондом с изменяемой плавучестью. Показана возможность эффективного управления скоростью погружения (всплытия) зонда с целью уменьшения динамической погрешности измерения вследствие различной инерционности используемых первичных преобразователей.
Ключевые слова: буй-профилемер, профильные измерения, управление движением, адаптивный регулятор, моделирование движения.
Постановка и решение задачи выведения морского профилемера на заданную глубину с нулевой вертикальной скоростью. Математическая модель морского профилемера как объекта управления, построенного на основе ныряющего буя с регулируемой плавучестью, представлена в работе [1] нелинейным дифференциальным уравнением (20) с начальными условиями (22) и параметрами (31), характеризующим процесс погружения. Поставим задачу (об адаптивном управлении с эталонной моделью) - найти управление u(t, z, z) в форме обратной связи, которое обеспечит перевод профилемера из начального состояния (22) в [1], соответствующего его неподвижному положению на поверхности в точке гидрологической станции, в новое неподвижное состояние стационарного равновесия z(t) = z ; Z(t) = 0, соответствующее заданной глубине z . При этом время перехода в требуемое состояние tp также должно задаваться. Далее потребуем, чтобы процессы z(t) ® z ;
z(t) ® 0 в замкнутой системе управления
z(t) + a,z2 (t) = -b0u(t, z, z) + fp (z), t e [t0, tf ], tp £ tf, (1)
полученной на основе (20) в статье [1], проходили в малой окрестности эталонных процессов z*(t) ® z ; z*(t) ® 0, формируемых эталонной моделью
z*(t) + aiz*(t) + a0z(t) = b0z; a0,a,b0 = const>0, (2)
© Л.А. Краснодубец, В.И. Забурдаев, В.В. Альчаков, 2012
70 ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 5
определяющей желаемые динамические и статические свойства проектируемой системы автоматического управления движением профилемера. Выбор эталонной модели в виде (2) обусловлен только порядком уравнения (1). Степень приближения эталонных процессов к соответствующим процессам в управляемой системе будем оценивать значением целевой функции
G« Л[.-(0-z(t,u)]2,, * О, (3)
где G(u) представляет собой нормированное по массе значение энергии ускорения [2], вычисляемой в окрестности траектории движения эталонной модели. Чем меньше величина G(u), тем в большей степени процессы в системе приближаются к процессам в эталонной модели. Целевая функция (3) при такой постановке служит одновременно и критерием управления и целевым условием адаптации. При этом в роли физически реализуемого аргумента целевой функции выступает текущее значение функции управления.
Решение поставленной задачи выполним следующим образом. Пусть начальное состояние модели (2) в точности соответствует начальному состоянию объекта (1), а степень приближения управляемых процессов к эталонным процессам соответственно в системе и эталонной модели оценивается значением целевой функции (3) в ходе функционирования системы автоматического управления.
В соответствии с градиентным методом отыскания min G(u) уравнение для вычисления оптимального управления примет вид
du . dG(u)
— = -1—:— , (4)
dt du
где l = const характеризует скорость, с которой управление u(t) приближается к оптимальному значению uopt. Соотношение (4) соответствует уравне-
m тт dG(u) нию скоростного градиента [3]. Далее, вычисляя производную - при
du
помощи (3) с учетом (1) и подставляя полученный результат в (4), можно получить уравнение, определяющее искомое управление в дифференциальной форме
^ = -ЛЬ0 [z * -z(t, u)]. (5)
dt
Анализ уравнения (1) показывает, что изменение какого-либо из параметров объекта a0, a1, Ь0(или всех вместе) вызывает изменение ускорения &z&(t,u) , что в соответствии с (5) приводит к вычислению нового управления u(t) , обеспечивающего адаптацию системы управления к новым условиям функционирования. По этой причине закон управления (5) можно назвать законом прямого адаптивного управления. Структура регулятора, реализующего такой закон, формируется в ходе преобразования уравнения (5) к виду,
ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 5 71
удобному для технической реализации. Сначала с учетом (1) и (2) выражение (5) принимает вид
Ши л Г ,_ * ч , * .."I
— = -М0[а0(г - г)-[ г- г\, (6)
ш
* . * V,
где г и г - переменные эталонной модели, которые принимают участие в формировании управляющей функции. Заменой в (6) этих переменных соответствующими переменными, поступающими от объекта управления г = г . *
иг = г, формируются цепи обратной связи. Далее путем интегрирования при нулевых начальных условиях обеих частей модифицированного уравнения (6) получается окончательная форма (аналитическая структура) закона управления для адаптивного регулятора в виде
и = -1Ь0 |а01[(г - -а1г - = -к|а0|[(г - г)Ш -а1г - .¿\|. (7)
По выражению (7) можно построить структурную схему адаптивного регулятора, которая изображена на рис. 1.
Р и с. 1. Структурная схема адаптивного регулятора
Задача синтеза параметров адаптивного регулятора сводится к выбору параметра 1, а также вычислению параметров а0 и а1 эталонной модели, которые одновременно являются параметрами регулятора и определяются требуемыми показателями качества переходной характеристики проектируемой системы управления. Они могут быть вычислены различными методами [4].
Исследуем влияние параметра 1 на свойства замкнутой системы управления с регулятором, реализующим сконструированный закон управления (7). Прежде всего отметим, что по построению дифференциальный закон управления (6) при заданных начальных условиях описывает траекторию
72
0233-7584. Мор. гидрофиз. журн., 2012, № 5
движения управляющей функции u(t) к ее оптимальному значению uopt. Далее введем в рассмотрение функцию
u(t) = u(t) - uopt, uopt = const, (8)
характеризующую текущее отклонение управляющей функции u(t) относительно ее оптимального значения uopt , которое можно найти из условия минимума целевой функции (3) как решение уравнения
^ = b0[z* -z(u)] = 0. (9)
ou
Действительно, из уравнения (9) следует, что
z* = z(uopt). (10)
В таком случае в соответствии с (1) для стадии погружения (u < 0) выражение для оптимального управления определяется соотношением
uopt =- b1 z&(uopt) - b- a1z2 + b- fp , uopt < 0 , b0 b0 b0
которое с учетом (10) преобразуется к виду
1 ... 1 ,2 1
— z--axz + —
b b ь,
uopt =-T-z*-T"a1z2 +T- fp, uopt < (11)
К аналогичному виду, определяющему функцию управления, можно преобразовать уравнение (1). Имеем
«(*) = --1 г - -1 а/ + -1 /, «(*) < 0. (12)
Ьа Ьо Ьо
Принимая во внимание (11) и (12), перепишем (8) с учетом «(*) < 0 следующим образом:
«(0 = -[«(0 - « ] = -1 г + -1 а1Ч2 - -1 / - -1 - -1 а1Ч2 + -1 / =--Цг* - Ч) .(13)
Ьо Ьо Ьо Ьо Ьо Ьо Ьо
Далее, дифференцируя обе части (8) и подставляя полученный результат вместе с (13) в уравнение (5), получим дифференциальное уравнение для отклонения и (*) в виде
+ 10« = 0; ы(0 = «0; * ]. (14)
Ш
Соответствующее характеристическое уравнение записывается следующим образом: 5 + 1Ь02 = 0, а его решение определяет характеристическое число 5 = -1Ь02 .
ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 5
73
Таким образом, решение уравнения (14) с указанными начальными условиями можно представить в виде
и = и0 ехр(-1Ь020, I е [¿0, ^ ]. (15)
Из выражения (15) следует, что при 5 = -ЩЬ\ для ^ отклонение
управляющей функции и (^) от ее оптимального значения иор(. по экспоненте
стремится к нулю. Причем скорость этого стремления пропорциональна величине
Щ = -, (16)
Та
где Та - постоянная времени процесса адаптации, который характеризует движение управляющей функции от ее начального значения до оптимального. При этом необходимым и достаточным условием устойчивости этого процесса является неравенство 1 < 0 .
Если задать желаемое время 1ра переходного процесса адаптации, то
значение постоянной времени адаптации определится соотношением
Т »
а ~ '
п
где п = 5...7 .
При этом следует отметить, что для эффективной работы регулятора постоянная времени процесса адаптации Та должна быть, по крайней мере, на порядок меньше постоянной времени Т основного процесса управления, т.е.
Та < 0,1т. (17)
Учитывая последний результат и формулу (16), выражение, которое определяет абсолютное значение параметра 1, можно записать в виде
Щ- -П- №
0 ра
Коэффициент усиления, входящий в выражение (11) в работе [1], с учетом (22) в [1] вычисляется по формуле
к=щ »п (1«)
0 ра
Можно показать [5], что при неограниченном росте значения параметра к = ЯЬ0 (коэффициента усиления адаптивного регулятора) уравнение динамики проектируемой системы управления вырождается в уравнение эталонной модели (2), которая по построению является устойчивой. По этой причине замкнутая система с адаптивным регулятором при к ® ¥ остается устойчивой.
74
0233-7584. Мор. гидрофиз. журн., 2012, № 5
Полученный результат свидетельствует о том, что коэффициент усиления к = Xb0 можно выбирать без учета параметров X и b0 при настройке адаптивного регулятора в процессе его наладки для работы с конкретным объектом управления в составе замкнутой системы. На практике начальное значение коэффициента усиления к можно принять равным единице, а далее экспериментально довести его до необходимого уровня.
Расчет параметров адаптивного регулятора a0 и а1 проведен по формулам, которые получены при помощи канонической формы записи колебательного звена, выполняющего функцию эталонной модели, и имеют вид
ц = А, ц = 2V-, (20)
t t
где t- постоянная времени процесса управления; д - декремент затухания (параметр, характеризующий колебательные свойства эталонной модели).
Если принять во внимание известное приближенное соотношение tp = (5...7)t, то формулы (20) можно преобразовать к виду
5...7 5...7
ц » ——, Ц »2д-, (21)
tp tp
где tp - время желаемой переходной характеристики проектируемой системы управления (длительность переходного процесса от начального состояния системы управления до момента достижения профилемером заданной глубины); д = 0,907 соответствует желаемому переходному процессу без перерегулирования.
Таким образом, применяя формулы (21) и учитывая неравенство (17) вместе с формулой (18), можно вычислить все параметры адаптивного регулятора а0,ц и X по всего лишь одной заданной величине tp - времени желаемой переходной характеристики. При этом время процесса адаптации выбрано из условия tpa = 0,0001 tp .
Процесс выведения морского профилемера на заданную глубину можно рассматривать как первую фазу его движения с целью подготовки к выполнению профильных измерений. Вторую фазу управляемого движения морского профилемера будем рассматривать как процесс профильных измерений, выполняемых в ходе управляемого всплытия с заданной глубины и с заданной вертикальной скоростью, которая удовлетворяет неравенству (4) в [1].
Постановка и решение задачи управления движением морского про-филемера в процессе выполнения термохалинных измерений. Вторая фаза движения профилемера, непосредственно связанная с выполнением измерений, начинается из неподвижного состояния на заданной глубине и продолжается в режиме всплытия с заданной постоянной вертикальной скоростью до момента достижения границы разделения сред (морской поверхности). Соответствующий закон управления найдем следующим образом.
Сначала построим математическую модель, характеризующую процесс изменения скорости движения профилемера на стадии всплытия. Для этой цели в уравнении (23) из работы [1] выполним замену переменных
ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 5 75
.. dv
z = — , & = v, dt
где v - скорость всплытия. Полученное неоднородное нелинейное дифференциальное уравнение
— = av2 - bu + fp, u > 0 (22)
dt Jp
в соответствии с описанием стадии всплытия следует рассматривать с начальными условиями
to = 0; v(to) = 0, (23)
при этом функции u и fp должны быть известными.
Поставим задачу найти управление (управляющую функцию u) в форме обратной связи, которое обеспечит перевод профилемера из начального состояния (2), соответствующего его неподвижному положению на заданной глубине (вблизи проекции точки гидрологической станции), в состояние движения (всплытия) с постоянной скоростью вплоть до момента tf выхода
на поверхность. При этом время tp достижения заданной скорости v следует
задать. Далее потребуем, чтобы процессы v(t) ® v , dv(t) ® 0 в замкнутой
dt
системе управления
d- = av2 -bu(t,v) + fp , te[t0, tf ], tp £ tf (24)
проходили в малой окрестности эталонных процессов v*(t) ® v ; v*(t) ® 0, формируемых эталонной моделью
v*(t) + a0v(t) = J30v; bo = a0 = const > 0 , (25)
определяющей желаемые динамические и статические свойства проектируемой системы автоматического управления движением профилемера. Выбор эталонной модели в виде (25) обусловлен не только порядком уравнения (22), но и выполнением равенства b0 = a0, что придает выбранной модели астатические свойства. Степень приближения эталонных процессов к соответствующим процессам в управляемой системе будем оценивать значением целевой функции
G(u) ^M-^rf^0- (26)
где G(u) представляет собой нормированное по массе значение энергии ускорения [2], вычисляемой в окрестности траектории движения эталонной модели (25). Вывод закона управления, определяющего искомую управляющую
76
ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 5
функцию, выполним по аналогии с решением задачи вывода буя-профилеме-ра на заданную глубину с нулевой скоростью.
Принимая во внимание, что в рассматриваемом случае
^ = ) - ^, и)],
ди
а также учитывая (4) и (25), уравнение для закона управления в дифференциальной форме можно записать в виде
^ = ЛЬ[а0(у -V*)- V]. (27)
Далее, выполняя замену V* = V в (27), что соответствует замыканию объекта управления вместе с регулятором отрицательной обратной связью, а также интегрируя обе части полученного уравнения при нулевых начальных условиях, можно получить аналитическое выражение для искомой управляющей функции в следующем виде:
í
и (*) = ЛЬ[а0 |(у - у)сИ - V]. (28)
0
Соответствующая структурная схема адаптивного регулятора, обеспечивающего управляемое движение морского профилемера в процессе термохалин-ных измерений с заданной постоянной скоростью в режиме всплытия, изображена на рис. 2.
V
1
ЛЬ
и
< V
Р и с. 2. Структурная схема адаптивного регулятора скорости всплытия морского профилемера Параметры регулятора а0 и к рассчитываются по формулам
«0 ="
5...7
к=
700
(29)
которые получены с учетом того, что эталонная модель описывается дифференциальным уравнением первого порядка, а время адаптации выбрано из
соотношения = 0,0.
ра р
¡ББЫ 0233-7584. Мор. гидрофиз. журн., 2012, № 5
77
Моделирование адаптивной системы автоматического управления.
Синтезированный адаптивный регулятор исследовался в среде визуального моделирования БтыНпк [6] совместно с предложенными математическими моделями морского профилемера, параметры которых определены соотношениями (32) в статье [1].
На первом этапе моделировался процесс вывода профилемера на заданную глубину. При этом использовался закон управления в виде (7). Параметры регулятора вычислены с помощью формул (17), (18), (21) и имеют следующие значения: 1 = 0,0000725; к = 0,0713; а0 = 0,000049; а1 = 0,0127.
Объект управления и подключенный к нему адаптивный регулятор охвачены отрицательной обратной связью. Структурная схема моделирования замкнутой системы управления процессом погружения морского профилемера приведена на рис. 3. Результаты моделирования иллюстрируют рис. 4 и 5. На первом изображен график процесса погружения на заданную глубину 500 м. На втором - график изменения скорости погружения. Как видно из графиков, профилемер выходит на заданную глубину за требуемое время и с нулевой вертикальной скоростью.
и« V- 'М
Л
Профилемер ►
1 V« 1 2(1) I |
в в
Р и с. 3. Структурная схема моделирования управляемого процесса погружения морского профилемера
г, м 0
100
200
300
400
500
600
0
200
400
600
800
1000 г, с
V, м/с 1.4
1.2
0.8
0.6
0.4
0.2
200
400
600
800
1000 г, с
Р и с. 4. График процесса погружения профилемера на заданную глубину с нулевой вертикальной скоростью
Р и с. 5. График скорости погружения морского профилемера
78
0233-7584. Мор. гидрофиз. журн., 2012, № 5
0
0
На втором этапе моделировался процесс управляемого движения профи-лемера при выполнении профильных измерений в режиме всплытия с начальной глубины 500 м с заданной постоянной скоростью 0,2 м/с. При этом использовалась «скоростная» модель профилемера (22) и адаптивный закон управления в виде (28). Структурная схема моделирования приведена на рис. 6. Параметры регулятора, вычисленные с помощью формул (29) для 1 = 2 с и 1 = 0,0к , имеют следующие значения: к = 0,24; а0 = 2,33 .
н>-
о
1 ли 1 1
в
0.2
9.821
Р и с. 6. Структурная схема моделирования управляемого процесса всплытия морского профилемера с заданной постоянной скоростью при выполнении термохалинных измерений
Графики скорости всплытия морского профилемера приведены на рис. 7, а на рис. 8 приведен график изменения глубины при выполнении профильных измерений.
V, м/с 0.25
0.2 0.15
0.1 0.05 0
0 1 2 3 4 5 0 500 1000 1500 2000 2500
I, с I, с
а б
Р и с. 7. Графики скорости всплытия морского профилемера при выполнении измерений во время начальной фазы (а) и в установившемся режиме (б)
На третьем этапе выполнялось сравнительное моделирование управляемого и неуправляемого движений морского профилемера в режиме всплытия с одинаковой глубины (500 м), которое сопровождалось скачкообразным изменением вертикальной стратификации водной среды. При этом начальные состояния и объемы камер плавучести управляемого и неуправляемого про-
0233-7584. Мор. гидрофиз. журн., 2012, № 5 79
филемеров принимались одинаковыми. В качестве возмущающего воздействия рассматривалось гипотетическое скачкообразное изменение на 50% плотности морской воды (рис. 9). Изменения скорости и глубины при всплытии неуправляемого и управляемого профилемеров показаны соответственно на рис. 10 и 11.
2, М
г, с
Р и с. 8. График изменения глубины при выполнении профильных измерений в режиме всплытия морского профилемера
рх100,% 0.6
0.5
0.4
0.3
0.2 0.1 0
0 500 1000 1500 2000 2500
г, с
Р и с. 9. График гипотетического скачкообразного изменения плотности морской воды
80
0233-7584. Мор. гидрофиз. журн, 2012, № 5
V, м/с 0.4 —
0.3 0.2 0.1 0
2, М 0
500
500 1000 1500 2000 2500 0 500 1000 1500 2000 2500
г, с л г, с
а б
Р и с. 10. Графики изменения модуля скорости (а) и глубины (б) при неуправляемом всплытии профилемера
V, м/с
0.2 ■ 0.15 0.1 0.05 0
2, М 0
500
500 1000 1500 2000 2500 0 500 1000 1500 2000 2500
г,с * г,с
а б
Р и с. 11. Графики изменения модуля скорости (а) и глубины (б) при управляемом всплытии профилемера
Заключение. Выполненные исследования показывают возможность повышения репрезентативности термохалинных измерений буями-профилеме-рами за счет управления скоростью движения по траектории. Система автоматического управления позволит с высокой точностью выводить буй-про-филемер на заданную глубину с нулевой вертикальной скоростью за требуемое время. Это позволит выполнять синхронные измерения несколькими про-филемерами. Автоматическая стабилизация скорости движения буя-профиле-мера относительно задаваемых расчетных величин снижает влияние инерционности датчиков на результаты измерений. Адаптивные свойства регулятора, управляющего плавучестью буя-профилемера, позволяют устранять влияние возмущающих сил, вызванных неоднородностями водной среды, на параметры движения, которое можно формировать при помощи программного устройства в соответствии с выбранной методикой измерений.
0233-7584. Мор. гидрофиз. журн., 2012, № 5
81
СПИСОК ЛИТЕРАТУРЫ
1. Краснодубец Л.А., Забурдаев В.И., Альчаков В.В. Управление морскими буями-профилемерами как метод повышения репрезентативности термохалинных измерений. Модели движения // Морской гидрофизический журнал. - 2012. - № 4. - С. 69 - 79.
2. Крутько П.Д. Новые технологии аналитического проектирования алгоритмического обеспечения систем управления движением // Управление, автоматизация и окружающая среда: Матер. междунар. науч.-техн. конф., г. Севастополь, 8 - 13 сентября 2008 г. - Севастополь: Изд-во СевНТУ, 2008. - С. 4 - 24.
3. Антонов В.Н., Терехов В.А., Тюкин И.Ю. Адаптивное управление в технических системах. - СПб.: Изд-во С.-Петербургского ун-та, 2001. - 244 с.
4. Краснодубец Л.А. Терминальное управление в морских наблюдательных системах с подвижными платформами сбора данных // Изв. РАН. Серия ТиСУ. - 2008. - № 2. -С. 141 - 153.
5. Краснодубец Л.А. Аналитическое конструирование адаптивных регуляторов на основе концепций обратных задач динамики и локальной оптимизации // Вестник СевНТУ. -Севастополь. - 2010. - Вып. 108. - С. 5 - 9.
6. Гультяев А. MATLAB. Имитационное моделирование в среде Windows: практическое пособие. - СПб.: Корона, 1999. - 288 с.
*Севастопольский национальный технический университет E-mail: lakrasno@gmail.com alchakov@gmail.com **Морской гидрофизический институт НАН Украины, Севастополь
E-mail: oaoimhi@inbox.ru
Материал поступил в редакцию 06.05.11 После доработки 18.08.11
АНОТАЦ1Я Представлена робота е другою частиною статл тих саме авторiв, виданоТ в «Мор-ському гiдрофiзичному журналЬ> в № 4 за 2012 р. У нш виконано аналiз нелшшних диференщ-альних рiвнянь, яю описують процес управлшня гiдрофiзичним зондом iз змшною плавучютю. Показана можливють ефективного управлшня швидюстю занурення (спливання) зонда з метою зменшення динамiчноТ погршносп вишрювання внаслщок рiзноТ шерцшносп використо-вуваних первинних перетворювачiв.
Ключовi слова: буй-профiлемiр, профшьш вишрювання, управлшня рухом, адаптивний регулятор, моделювання руху.
ABSTRACT The present work is the second part of the paper by the same authors published in «Marine Hydrophysical Journal», № 4, 2012. Non-linear differential equations describing the process of control of a hydrophysical probe with variable buoyancy are analyzed. Possibility of efficient control of the probe submergence (emergence) speed aimed at reducing the measurement dynamic error resulting from different inertance of the applied primary converters is shown.
Keywords: buoy-profilemeter, profile measurements, control of movement, adaptive regulator, modeling of movement.
82
ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 5