СИСТЕМА ОПРЕДЕЛЕНИЯ ГИДРОДИНАМИЧЕСКИХ КОЭФФИЦИЕНТОВ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ДВИЖЕНИЯ СУДНА
SYSTEM FOR DETERMINING HYDRODYNAMIC COEFFICIENTS OF MATHEMATICAL MODEL OF SHIP MOVEMENT

Изобретение относится к области судовождения - автоматическому управлению движением судна, в частности к идентификации гидродинамических коэффициентов математической модели движения судна.
Известны системы определения параметров математической модели движения судна, которые основаны на использования заданной структуры математической модели движения судна, вида:
x
y
z
где V, V, ω, φ - продольная, боковая и угловая скорости судна, угол курса, δ - угол перекладки руля,
ij
11
12
13
21
22
23
a - неизвестные параметры математической модели движения судна (подлежащие уточнению гидродинамические коэффициенты: a, a, a, a, a, a),
1
2
3
i
ς, ς, ς - входные помехи с нулевыми средними и известными ковариационными матрицами  M[ς(t)]=0 и
[П.В Попов Идентификация уравнений динамики подвижного аппарата по натурным данным ходовых испытаний ОАО «Концерн НПО «Аврора» С.Петербург // Сборник трудов ХХХУ111 Всероссийской конференции издательство ИПУ РАН 2012, г. Москва]
Известна система экспериментального определения влияния вариации коэффициента математической модели объекта на его динамику (SU 205913 А1, 01.02.1967). Выходной сигнал с исследуемого объекта запоминают (задерживают), формируют невязку, а затем пропускают через электронную модель исследуемого объекта, в котором варьируют величину коэффициента (постоянной времени) модели. Недостатком такой системы является идентификация только одного параметра идентифицируемой системы.
Известна система экспериментального определения параметров математической модели движения судна (RU 2151713 C1, 27.06.2000, принятый нами в качестве прототипа). Известная система использует только часть измеряемых фазовых координат состояния судна. Используя измеренные параметры и сформированные с использованием электронной модели движения судна, получают невязки, которые минимизируют путем вариации искомых величин гидродинамических коэффициентов в электронной модели.
Недостатки:
1) использование только части измеряемых фазовых координат состояния судна,
2) низкая точность идентификации гидродинамических коэффициентов математической модели движения судна,
3) высокие требования к виду задаваемых специальных режимов плавания для идентификации параметров гидродинамических коэффициентов математической модели.
Техническим результатом предлагаемой системы определения гидродинамических коэффициентов математической модели движения судна является:
- создание высокоточной математической модели движения судна, параметры которой уточняются в процессе рейса судна,
- повышение качества автоматического управления движением,
- повышение безопасности проводки судна в узкостях,
- снижение загрузки рулевого привода при сильном волнении.
Поставленные цели достигаются с помощью:
- введения подсистемы формирования коэффициентов усиления в процессе идентификации гидродинамических коэффициентов (введением блоков 10. и 11.);
y
z
- использования полного набора измерителей фазового состояния судна (включая измерители линейных dV/dt и угловых dω/dt ускорений);
- проведения идентификации при облегченных идентификационных режимах движения судна,
- определения значений гидродинамических коэффициентов путем решения двух независимых систем алгебраических уравнений.
Предлагаемая система идентификации позволяет в процессе рейса судна (практически не нарушая заданный режим движения) уточнять гидродинамические коэффициенты математической модели движения судна.
Математическая модель принята вида:
x
y
z
где V, V, ω - продольная, боковая и угловая скорости движения судна,
y
dV/dt - боковое (линейное) ускорение судна,
z
dω/dt - угловое ускорение судна,
δ - кормовой руль,
ij
11
12
13
21
22
23
a - неизвестные, подлежащие уточнению (шесть оценок) гидродинамические коэффициенты: a, a, a, a, a, a путем решения систем алгебраических уравнений.
yi
yi
xi
zi
i
zi
Технический результат достигается тем, что система идентификации гидродинамических коэффициентов математической модели движения судна содержит задатчик идентификационных маневров управления движением судна; блок измерения, включающий датчики бокового ускорения - dV/dt, боковой скорости судна - V, продольной скорости - V, угловой скорости - ω, угла перекладки руля - δ, углового ускорения - dω/dt в момент времени i, в блоке формируются также сигналы: .,
yi
zi
dV/dt, dω/dt; блок формирования правой части уравнения математической модели боковой скорости:
блок формирования правой части уравнения математической модели угловой скорости:
блок формирования невязок канала боковой скорости:
блок формирования невязок канала угловой скорости:
блок формирования коррекции оценок гидродинамических коэффициентов в канале боковой скорости:
блок формирования коррекции оценок гидродинамических коэффициентов в канале угловой скорости:
11,i
12,i
13,i
21,i
22,i
23,i
11,i
12,i
13,i
22,i
23,i
33,i
блок памяти, содержащий оценки гидродинамических коэффициентов судна a, a, a, a, a, a и элементы ковариационной матрицы: P, P, P, P, P, P;
11,i
12,i
13,i
21,i
22,i
23,i
11,i+1
12,i+1
13,i+1
21,i+1
22,i+1
23,i+1
ij
блок сравнения, формирующий модули разностей гидродинамических коэффициентов судна, используя сигналы из блока памяти: a, a, a, a, a, a и из блоков формирования коррекции гидродинамических коэффициентов каналов боковой и угловой скоростей: a, a, a, a, a, a. В блоке сравнения производится сравнение модулей разности с заданными допустимыми значениями - ε:
Для решения поставленных задач вводят блок формирования коэффициентов усиления в процессе идентификации гидродинамических коэффициентов судна:
и блок формирования коррекции элементов ковариационной матрицы:
Обеспечить правильное функционирование всех рассмотренных выше блоков возможно, если выполнить все связи между ними в соответствии с фиг.1.
Если все шесть модулей разности, сформированные в блоке сравнения:
ij
y
z
11,i+1
12,i+1
13,i+1
22,i+1
23,i+1
33,i+1
11,i
12,i
13,i
22,i
23,i
33,i.
11,i+1
12,i+1
13,i+1
y
21,i+1
22,i+1
23,i+1
z
11,i
12,i
13,i
21,i
22,i
23,i
меньше допустимых значений ε, то полученные оценки гидродинамических коэффициентов, в блоке формирования коррекции оценок гидродинамических коэффициентов в канале боковой скорости - V и блоке формирования коррекции оценок гидродинамических коэффициентов в канале угловой скорости - ω близки к истинным значениям искомых параметров математической модели движения судна в данный момент времени плавания, если же хотя бы одно из шести условий, полученных в блоке сравнения, не удовлетворяется, то шесть элементов ковариационной матрицы P, P, P, P, P, P из блока формирования коррекции элементов ковариационной матрицы вводят в блок памяти взамен шести элементов ковариационной матрицы P, P, P, P, P, P. Три оценки гидродинамических коэффициентов a, a, a из блока формирования коррекции оценок гидродинамических коэффициентов в канале боковой скорости - V и три оценки гидродинамических коэффициентов a, a, a из блока формирования коррекции оценок гидродинамических коэффициентов в канале угловой скорости - ω вводят в блок памяти взамен a, a, a, a, a, a.
При поступлении новых измерений через интервал времени Δt рассмотренный выше процесс уточнения гидродинамических коэффициентов вновь повторяется. По истечении «n» таких циклов идентификации, когда выполнятся все условия сравнения  полученные значения оценок гидродинамических коэффициентов будут соответствовать искомым, т.е. определение гидродинамических коэффициентов математической модели движения судна успешно заканчивается.
Построение и работа системы определения параметров математической модели движения судна
y
y
x
z
z
y
z
Блок-схема системы идентификации приведена на фиг.1 и включает: 1 - задатчик идентификационных маневров управления движением судна; 2 - судно; 3 - измерительный блок (включающий датчики бокового ускорения dV/dt, боковой скорости хода V, продольной скорости V; угловой скорости ω, угла перекладки руля δ и углового ускорения dω/dt); 4 и 5 - блоки формирования правых частей математической модели боковой скорости судна - V и модели угловой скорости судна - ω; 6 и 7 - блоки формирования невязок каналов боковой и угловой скоростей судна; 8 и 9 - блоки коррекции гидродинамических коэффициентов в каналах боковой и угловой скоростей; 10 - блок формирования коэффициентов усиления в процессе идентификации гидродинамических коэффициентов; 11 - блок коррекции элементов ковариационной матрицы; 12 - блок памяти; 13 - блок сравнения.
Рассматриваемая система может быть реализована как с применением аналоговых счетно-решающих элементов, так и цифровых (последние предпочтительней). Измерительный блок следует оснастить датчиками с точностью не ниже 2%.
Рассмотрим процесс формирования уравнений для определения (корректировки) гидродинамических коэффициентов математической модели движения судна.
Процесс вычисления оценок гидродинамических коэффициентов начинается с задания:
- начальных значений гидродинамических коэффициентов (в блоке памяти 12):
11,i
12,i
13,i
21,i
22,i
23,i
a=0., a=0., a=0., a=0., a=0., a=0,
- начальных значений элементов ковариационной матрицы
11,i
22,i
33,i
12,i
13,i
23,i
4
P=P=P=10, P=P=P=0 (в блоке-памяти 12),
- идентификационных маневров судна, которые формируются в задатчике идентификационных маневров движения судна 1.
11,i
12,i
13,i
21,i
22,i
23,i
Формирование двух каналов получения оценок гидродинамических коэффициентов: а) a, a, a и б) a, a, a.
y
11,i
12,i
13,i
y
а. Сигналы  с блока измерения 3 поступают на первый вход блока формирования правой части математической модели боковой скорости судна V 4, на второй вход которого поступают из блока памяти 12 гидродинамические коэффициенты a, a, a. В блоке 4 формируется сигнал, описывающий правую часть уравнения боковой скорости V:
11,i
12,i
13,i
22,i
23,i
33,i
1i
2i
3i
С блока измерения 3 сигналы  поступают на первый вход блока 10, на второй вход которого из блока памяти 12 поступают элементы ковариационной матрицы P, P, P, P, P, P. В блоке 10 формируются коэффициенты усиления процесса идентификации гидродинамических коэффициентов K, K, K.
yi
1i
2i
3i
y
С выходов блоков 3, 4 и 10 сигналы: dV/dt,  и K, K, K поступают на вход блока невязок канала боковой скорости - V 6, где формируют невязки:
11,i
12,i
13,i
Невязки, сформированные по зависимости (2) из блока 6, поступают на первый вход блока коррекции гидродинамических коэффициентов 8, на второй вход которого поступают гидродинамические коэффициенты a, a, a из блока памяти 12. В блоке 8 формируют обновленные гидродинамические коэффициенты:
Гидродинамические коэффициенты, сформированные по зависимости (3) в блоке 8, вводятся в блок сравнения 13.
z
21,i
22,i
23,i
z
б. Сигналы  снимаемые с блока измерения 3, поступают на первый вход блока формирования правых частей математической модели углового движения - ω 5, на второй вход которого поступают из блока памяти 12 гидродинамические коэффициенты a, a, a. На выходе блока 5 формируется сигнал, описывающий правую часть уравнения углового движения - ω:
zi
1i
2i
3i
С выходов блоков 3, 5 и 10 сигналы: dω/dt,  и K, K, K поступают на вход блока невязок канала угловой скорости 7, где формируют невязки:
21,i
22,i
23,i
Невязки с выхода блока 7. поступают на вход блока коррекции гидродинамических коэффициентов 9, на второй вход которого из блока памяти 12 поступают гидродинамические коэффициенты a, a, a.
В блоке 9 формируют обновленные гидродинамические коэффициенты:
11,i+1
12,i+1
13,i+1
21,i+1
22,i+1
23,i+1
11,i
12,i
13,i
21,i
22,i
23,i
Гидродинамические коэффициенты, сформированные по зависимости (4) в блоке 9, вводятся в блок сравнения 13. В блоке сравнения 13 формируются модули разностей обновленных гидродинамических коэффициентов: a, a, a, a, a, a и гидродинамических коэффициентов: a, a, a, a, a, a (вводят из блока памяти 12):
11,i+1
12,i+1
13,i+1
21,i+1
22,i+1
23,i+1
Если все разности (5) и (5a) не выходят из заданной точностной трубки, то процесс идентификации считается законченным, а полученные коэффициенты: a, a, a, a, a, a соответствуют истинным значениям гидродинамических коэффициентов математической модели движения в данный момент времени.
11,i
12,i
13,i
21,i
22,i
23,i
11,i+1
12,i+1
13,i+1
21,i+1
22,i+1
23,i+1
Если хотя бы одна из зависимостей (5) или (5a) не удовлетворяется, то в блок памяти 12 вместо (старых) значений гидродинамических коэффициентов: a, a, a, a, a, a заносятся новые: a, a, a, a, a, a (из блока 13).
1i
2i
3i
11,i+1
12,i+1
13,i+1
22,i+1
23,i+1
33,i+1
С измерительного блока 3 сигналы  поступают на первый вход блока коррекции элементов ковариационной матрицы 11, на второй и третий входы которого поступают коэффициенты усиления невязок процесса идентификации гидродинамических коэффициентов K, K, K из блока 10 и элементы ковариационной матрицы P, P, P, P, P, P из блока 12. На выходе блока 11 формируют обновленные элементы ковариационной матрицы:
где  Элементы ковариационных матриц вводят в блок памяти 12.
Затем через интервал времени Δt весь рассмотренный выше цикл уточнения гидродинамических коэффициентов вновь повторяется. Через «n» циклов идентификации, когда выполнятся все условия (5) и (5а), полученные значения оценок гидродинамических коэффициентов будут соответствовать искомым, таким образом (в данный момент времени) идентификация гидродинамических коэффициентов математической модели движения судна успешно заканчивается.
Результаты моделирования
11
12
13
21
22
23
x
На фиг.2 и 3 представлены оценки и истинные значения гидродинамических коэффициентов математической модели судна «Сиберяков» a, a, a, и, a, a, a при скорости судна V=5 m/c с использованием идентификационного режима «зигзаг».
Из фигур следует, что процесс корректировки оценок гидродинамических коэффициентов занимает около восьми секунд (расхождение восстановленных оценок гидродинамических коэффициентов от истинных значений не более 5%).