Построение математической модели корабельной платформы для подавления эффектов качки
Bui Van Tam
Новосибирский государственный технический университет, Новосибирск, Россия
Аннотация: Данная статья посвящена стабилизации верхней части механической платформы, находящейся на палубе корабля. Цель стабилизации состоит в обеспечении стабильного горизонтального положения различных измерительных приборы. Система стабилизации должна подавлять движения в двух степенях свободы, которые порождают наклоны платформы в двух осях: вдоль и поперек корабля. Статья анализирует известные методы стабилизации, осуществляет выбор метода решения. Также в статье даются обоснования для вывода математической модели системы и предлагается модель в виде структурной схемы. Предложен также метод проектирования регулятора и предварительные результаты моделирования.
Ключевые слова: стабилизация, автоматика, управление, регулятор, оптимизация, динамическая ошибка, статическая ошибка, моделирование
Введение
Качка любого морского судна отрицательно влияет на работу многих приборов. Особенно высокоточные приборы желательно
использовать в режиме без качки, или с существенным ослаблением её влияния. Для подавления качки оборудования, требуется разработка стабилизированной платформы, размещаемой на корабле. Если результаты измерений углов наклона такой платформы использовать в контуре обратной связи, можно обеспечить поддержание исходного
невозмущенного состояния рабочей поверхности платформы, чтобы ослабить или полностью устранить влияние качки на точность работы приборного оборудования. Для решения этой задачи желательно знать математическую модель движения такой платформы, т. е. связь между управляющими сигналами и перемещениями на выходе платформы, а также связь между этими перемещениями и сигналами на выходе датчиков углов наклона. Данная статья посвящена выбору технического решения для стабилизации верхней плоскости платформы, находящейся на корабле в условиях морской качки.
1. Постановка задачи
Требуется разработать эскизный проект для решения задачи подавления влияния качки на рабочую поверхность активной платформы. С этой целью необходимо выбрать метод подавления качки, вид датчиков и исполнительных устройств, построить математическую модель системы и рассчитать регулятор для этой системы.
2. Известные методы решения задачи
Существует несколько вариантов стабилизируемой платформы для использования её на корабле.
Первый вариант. Для стабилизации платформы, размещенной на корабле, важно изолировать эту платформу от борта с помощью системы, которая может изменять длину, когда плоскость борта меняет угол под действием внешних сил. Это означает, что система может компенсировать углы отклонения этой плоскости при изменении положения борта. Система выполняет эту функцию, как правило, используя гидроцилиндры в качестве привода, показанного на Рис. 1.
Второй вариант. Стабилизационная платформа изолирована от палубы подушками, состоящими из пружины и гидравлической плунжерной системы. Система похожа на подвеску, принимает на себя, но не передаёт на верхнюю столешню девиации силы сопротивления опоры от палубы. Системы, работающие по этому принципу, показаны на Рис. 2.
Обсуждение. Безопорный подвес платформы, по-видимому, является более эффективным, и одновременно более энергоёмким методом стабилизации. Такой подвес целесообразен в дорогостоящих стационарных установках, например, в лазерных системах, которые требуется полностью изолировать от механических воздействий от грунта. Магнитный подвес может обеспечить наилучшую изоляцию от высокочастотных вибраций, поэтому активное подавление вибраций требуется в меньшей полосе частот, чем при подвесе на гидравлических опорах.
Формально можно предложить также подвеску стола на воздушной подушке, это также будет энергозатратный вариант, который
обеспечит лучшее подавление высокочастотных колебаний основания. Можно также предложить, например, стабилизацию на струях воды, что также является энергозатратным вариантом. Для стабилизации оптической плиты в опыте Майкельсона-Морли было
использовано дорогостоящее решение, которое состояло в том, что гранитное основание интерферометра плавало в ванне со ртутью. Этот вариант также неприемлем для использования на корабле.
Рис. 1. Устойчивая модель стабилизированной платформы с использованием конструкции гидроцилиндра
Рис. 2. Схема и внешний вид стабилизированной платформы с использованием электромагнитного механизма поддержки
Исходя из энергетических и прочих затрат следует предпочесть для корабля активную стабилизацию на гидравлических цилиндрах или иных аналогичных устройствах, т. е. по первому варианту.
3. Предложение по техническому
РЕШЕНИЮ ЗАДАЧИ
3. 1. ПРОСТРАНСТВЕННАЯ МОДЕЛЬ СТАБИЛИЗИРОВАННОЙ ПЛАТФОРМЫ
Предлагается создание стабилизированной платформы на четырех опорах. Некоторые из опор являются гидравлическим приводом, их длина может изменяться за счет управления сервоклапаном.
Простейший вариант стабилизированной платформы содержит, как минимум, две гидравлических опоры, как показано на Рис. 3.
Плоскость ОАВ - стабилизированная часть платформы. Эта плоскость предназначены для размещения на ней приборов, которые должны
быть защищены от корабельной качки. С этой стабилизированной частью платформы связана система координат ОXYZ. Точка О верхней части платформы закреплена подвижно на шарнире.
Плоскость О'А'В' связана с палубой корабля. Начальное расстояние между двумя плоскостями, измеренное на линии АВ, где размещены изменяемые опоры, обозначим к. Это расстояние, равное длине неизменяемых опор на противоположной стороне подвижной части. Стабилизированная платформа связана с палубой в точках О, А, В через вертикальную стойку ОО' и два гидравлических опор подножки АА' и ВВ' (длина двух гидравлических опор равны 1\, 12).
Плоскость ОАВ может вращаться вокруг оси ОХ и ОУ (но не может вращаться вокруг оси О2), сочленения А, В, А', В' являются шаровыми опорами, которые свободно вращаются вокруг осей ОХ, ОУ, О2.
Y
'A
Рис. 3. Модель стабилизированной платформы с двумя гидравлическими опорами
В системе координат ОХУ2 с началом О, ось ОХ совпадает с продольной осью корабля, ось ОУ параллельна плоскости палубы и влево, ось, О2 повернётся вниз (перпендикулярно плоскости ОАВ).
Верхняя подвижная часть платформы и установленный на ней прибор должны оставаться горизонтальными. В любой точке верхней части, в частности, вертикальное ускорение должно быть одинаковым. Пусть точка М - это пересечение медиан треугольника ОАВ. В этом случае если к этой точке приложить вертикально направленное усилие (вес установленного прибора), то в условиях равновесия это усилие распределяется на все
три опоры поровну. В точках О, А, В нагрузка на соответствующие опоры должны быть одинаковыми.
Пусть 01, 02 - углы поворота плоскости ОАВ вокруг оси ОУ и ОХ в системе координат ОХУ2.
Задача управления платформой может быть сведена к тому, чтобы найти необходимую длину изменяющихся гидравлических опор (и поддерживать её в найденном значении). Эта длина должна соответствовать углам поворота 01, 02 (при отклонении стабилизированной платформы от горизонтальной плоскости). Поэтому требуется найти функцию:
fi i = f AA) )l 2 = f AA)
(1)
То есть, когда стабилизированная платформа наклонена по сравнению со стандартной горизонтальной плоскостью с соответствующими углами в1, в2, требуется воздействовать на электрогидравлические клапаны, чтобы изменить длины (11, 12) опор, чтобы стабилизированная часть платформы оставалась параллельной горизонтальной плоскости.
3. 2. Вычисление длины опор
Исходное состояние стабилизированной части платформы - расположение параллельно нижней плоскости (основанию) и параллельно глобальной горизонтальной плоскости. Плоскость ОАВ связана с плоскостью палубы. Координаты точек А и В соответственно
А0(ХА0' уА0' 2 АО) , В0(ХВ0' уВ0' 2Во) Здесь:
ХА0 = ХВ0 = а'
Ъ
Уао = -ув =
2А0 = 2В0 = 0 (2)
В первом случае стабилизированная платформа наклонена на угол 01.
Затем, чтобы стабилизированная часть платформы вращалась для поддержания состояние равновесия, нужно повернуть плоскость ОАВ вокруг оси ОУ на угол 01. Тогда
A ^ A
B ^ B
соответственно
A1(XA1, yA1, ZA1) , b1(xb1, yb1, zb1)
Здесь:
— — cos 0-;
_ _b
Улг — -Увг — 2'
zai — zbi — a• sin
(3)
Во втором случае стабилизированная платформа продолжит наклонена на угол 02 (стабилизированная платформа вращается вокруг оси ОХ) показано на Рис. 5 а, Ь, с.
Рис. 4.а Взаимосвязь длин I и углов 0 когда стабилизированная часть платформы наклонена на угол 01
Ув
-i Э-
b
O
е-
УА Y
B
A
B-i-^
ZB1 - ZA1
V)-Ai
Z
h
Рис. 4.Ь. Взаимосвязь I и 0, когда стабилизированная часть платформы наклонена на угол 01
Рис. 5.a. Взаимосвязь l и в, когда стабилизированная часть платформы наклонена на угол в\, вг
Когда стабилизированная часть платформы в2. Чтобы привести плоскость OAB в продолжает наклоняться вокруг оси OX на угол равновесное состояние, необходимо повернуть
плоскость ОАВ вокруг оси ОХ на угол 02 в противоположном направлении, тогда А ^ А,
В ^ щ с соответствующими координатами
А2(ХА2> уА2, 2А2) , В2(ХВ2, уВ2' 2В2) • Здесь:
ХА2 = ХА1; УА2 = |c0s °2 - ZA1 sin °2; zA2 = zA1 C0S °2 + Ь sin ^2 XB2 = *4i; Ув2 =-^c0s^2 - zA1sin^2; zB2 = zA1 C0S ^2 - ^sin^2
(4)
Рис. 5.Ь Взаимосвязь I и 0 когда стабилизированная часть платформы наклонена на угол 01, 02 С другой стороны, имеется соотношение:
2а 1 = 2 в 1 = я-втб! (5)
Из выражений (3), (4) и (5) получим новые координаты точек А2 и В2 по оси ОХ в виде:
xA2 = a. cos ду, yA2 = ^ cos - a. sin sin 02; zA2 = a.sin cos д2 + Ь sin 02 xB2 = a.cos^1; yB2 =- — cos^2 - a. sin sin 02, zB2 = a. sin cos 02 - bsin02
Таким образом, когда плоскость ОАВ вращается вокруг ОУ и ОХ соответственно на угол 01 и 02 тогда точка А движется из точки
А0(ХА0' Уа0' 2А0) до А2(ХА2' УА2' 2А2) , точка В движется из точки Во(Хщо,Уво, 2во) до
(6)
В2( Хщ 2, У В 2, 2) соответственно изменения
длины цилиндров (1 01 = 1 02 = Ь ).
1 01 ^ 1
1
02 ^ 1 2
ABy
В
iT X >—
ABX
В2
Ув = Ув1
УВ2
X
AAy
XA2 = XAI *В2 = Xai
tLJ>
Хв2
*A2
Э-
(AA2)2 = (AAx)2 + (AAy)2 + (AAZ)2 (BB2)2 = (ABx)2 + (ABy)2 + (ABz)2
A2
O
УА2
*
b/2
A
Ai
УА = УА1
AA,
I Y
Рис. 5 с. Взаимосвязь I и 0, когда стабилизированная часть платформы наклонена на угол 01, Из предложенных схем мы можем рассчитать длину цилиндров в виде:
1 2 = (ЛЛ2)2 = (ЛАх)2 + (ЛАу)2 + (ЛЛ2)2
Здесь:
12 = (BB2 )2 = (ABx )2 + (ABy )2 + (ABZ )2
(7)
AAX = a - xA2 = a(1 - cos A)
AAy = b - yA2 = b - b cosP2 + a. sin A1 sinP2
AAz = h - Za2 = h - a. sin A1 cos P2 - b sin P2
(8)
ABX = a - xB2 = a(1 - cos A1)
ABy = - Ь - yB2 = - Ь + Ь cosP2 + a.sinP1sinP2
ABZ = h - zB2 = h - a.sinp cosP2 + Ь sinP2
(9)
В результате подстановки выражений (8) и (9) в выражение (7), получим:
а
b
12 = (a - acos61)2 + (Ь - Ьcos62 + a.sin61 sin62)2-
b 2
+ (h - a.sin 61 cos62 -—sin62)
b b
12 = (a - a cos 6^ + (---- — cos 02 + a.sin61 sin62)
2 2
Пусть:
b2
+ (h -a.sin61 cos62 + — sin62 )
b b
/1(6162) = Sqrt[(a - a cos61) + (— - — cos62 + a.sin 61 sin 62) +
b2
+ (h - a.sin 61 cos 62 -—sin62) ]
/2(6162) = Sqrt[(a - a cos61)2 + (-b + b cos62 + a.sin 61 sin 62)2 +
b2
+ (h - a.sin 61 cos 62 + —sin62 ) ]
(11)
Поскольку длины цилиндров является положительными числами 11, 12, тогда функции / (0^2), / (002) также имеет положительные значения, система уравнений (10) переписывается следующим образом:
11 = /1(6162) 1 2 = /2(6162)
(12)
4. Система управления гидравлической СИСТЕМОЙ С СЕРВОКЛАПАНОМ
Функциональная схема гидросистемы с сервоклапаном представлена на Рис. 6.
Предположим, что поршни двигаются по оси координат ОZ. Начало координат является точкой нормального равновесия. Когда борт корабля отклоняется на угол, положение плоскости также смещается на соответству-
ющую величину. Гидравлический привод обеспечивает движений для компенсации смещения, вызванного движущейся плоскостью борт. Построим уравнение движения поршней под действием сил, действующих от воздействия разности давления, поступающего в два отсека цилиндра, от давления массы объекта, помещенного на плоскость, а также воздействия волн и ветра (вертикальный компонент ускорения)
Принципиальная схема гидравлической системы, выполняющей функции стабилизации стабилизирующей платформы, может быть построена в соответствии с приведенной выше структурной схемой, показанной на Рис. 7.
Система, используемая для управления гидроцилиндром с сервоклапаном, показана на Рис. 8.
Рис 6. Модель управления гидроцилиндрами
Рис 7. Принципиальная схема гидравлической системы: 1 - система теплообмена; 2 - фильтр 3 - датчик температуры гидравлического масла; 4 - гидравлический масляный насос; 5 - фильтр; 6 - перепускной клапан; 7 -односторонний клапан; 8 - электрогидравлический сервоклапан; 9 - элемент измерения смещения цилиндра; 10 -исполнительная структура (гидравлический цилиндр); 11 - датчики давления; 12 - вычислительный прибор создает правила управления
Рис. 8: Структурная схема гидравлической системы с сервоклапаном: 1 - цилиндр; 2 - сервоклапан; 3 - усилитель и преобразователь; 4 - датчик положения; 5 - компьютер
5. Математическое описание
ГИДРАВЛИЧЕСКОЙ СИСТЕМЫ С
СЕРВОКЛАПАНОМ
Он состоит из сервоклапана, симметричного гидравлического цилиндра и других принадлежностей (Рис. 8). В секции сервоклапана смещение золотникового клапана ху связано с текущим ^ посредством функции первого порядка. Это сделано в публикации [3]
Уравнение баланса напряжения:
ТцХу Xv + kv iv ,
(13)
где zv - положение золотника сервоклапана [да], iv - входной ток моментного двигателя [тА], ту и к„ - постоянная времени и коэффициент усиления сервоклапана соответственно. Поток нагрузки Q1 связан со смещением золотникового клапана х^, и давлением нагрузки Рн . Это сделано в публикации [4], [5]
Qi = Cd шх !- (Ps - sgn(x) Ph) , (14)
где Ql - скорость потока нагрузки [т3/с], С -расход коэффициент, т - градиент площади
золотника сервоклапанаюмщ, Ps - давление подачи насоса [Ж/да2], Ph - перепад давления между PA и PB [N/да2], и р плотность гидравлического масла [кг/да3]. Применяя закон непрерывности каждой камеры привода, скорость потока нагрузки. Уравнение неразрывности дано как [3].
В уравнении (14) поток нагрузки Qi является линейным по отношению к нулевому отверстию и нулевому давлению нагрузки:
Qi = KqZv - Kc Ph
(15)
где К, = / дxv и К = -^1 / дРн, Са -коэффициент расхода, ю - градиент площади, К, и Кс - коэффициент усиления потока, отнесенный к углу открытия и усиление потока, отнесенное к давлению. Пренебрегая эффектами внешней утечки, динамику привода можно описать:
V
а=Ahzh+-^ph+chlph
4EU
(16)
где Ah - площадь поршня, Vh - объем камеры поршня, Eh - эффективный объемный модуль, - общий коэффициент утечки. Согласно принципу баланса сил, динамика штока поршня может быть описана следующим образом:
AhPh = mhZh+BhZh+Fh
(17)
где Bh - совокупный коэффициент смоделированных сил демпфирования и вязкого трения, Fh - внешняя нагрузка цилиндра, которая также является движущей силой, действующей на поверхность управления.
Сочетая (13)-(17), математическая модель динамики систему управления цилиндрами (СУЦ) может быть записана в виде структурной схемы, показанной на Рис. 9.
и Kv xr
h 1
Ка
кс
к,
Л
4К
S I С,„
Ра —►
Ah
Ah
Fh
-Ю
1 Zh 1
щР+Щ, 5
Zb
Рис. 9. Структурная схема модели СУЦ
В этой части проектирование контроллера по методу, изложенному в работе [2]. Путь х1 = 1к\х2 = ¿к;х3 = Рк;х4 = ху, как переменная состояния, систему легко описать с нелинейной моделью в пространстве состояний четвертого порядка. Комбинируя уравнение (3)-(5),
получаем систему уравнений гидравлической системы с сервоклапаном
А„
- _ ih _Bh
Х2 — Х3 Х2 '
mh mh mh
F,
h
X3 — y^
- sgn(x4)x3 -
4EhChi
Vh
Xo -
4EhAh Vh
(18)
X 4 — X4 + "
у — xi
Здесь: С —
IT
Это система нелинейных дифференциальных уравнений.
6. ПРОЕКТИРОВАНИЕ КОНТРОЛЛЕРА ЭЛЕКТРОГИДРАВЛИЧЕСКОГО СИСТЕМЫ НА ОСНОВЕ МЕТОДА BACKSTEPPING
В этом части синтезирование контроллера на основе по методу, изложенному [2].
Сочетая (13), (15), (16) и (17) мы получим линейную модель электрогидравлической системы, записанную в пространстве состояний
Х1 — х2
х = —X *2 mhХз
V
™-h
Х3 — Vh ХА
4Eh(kc+Chl) , Vh '
4EhAh
Vh
X2 (19)
1
k,,
Хл — Хл + и
У — х-
Обозначим e1 = x1 - Xld (для i = 1, 2, 3 и 4) ошибки слежения, т. е. разницы между каждой фактической переменной состояния и ее желаемым значением, где xld - желаемое положение, которое необходимо достичь.
Синтез контроллера гидросистемы проходит в четыре шага:
Первый шаг. Путь У1 = 1е2 как функция
Ляпунова для подсистемы х1 = х2. Цели этого шага - убедиться, что фактическая траектория соответствует желаемой траекторией. Это означает, что ошибка е1 сходится нулю.
Выполнение дифференциального уравнения У1 = ~е1, мы получим У1 = е1(х1 — х1а) =
е1(е2 + Х2й — хы) (20)
Сейчас мы выбираем хм как виртуальный управляющий сигнал, определяемый выражением (20):
х2й = Ф(е1,х1б) = хы — Ке1 (к1 > 0). Подставляя x2d в уравнение (20), получаем: У1 = —к1е1! + е1е2 (21) Если e2 = 0, то У1 = —к1е^ <0 и е1 будет сходиться к нулю по экспоненциальному закону
Второй шаг. Путь У2
е2+-е22
2 1 2 2
как
функция Ляпунова для подсистемы
Хп —
Ah mh'
Bh_ mh '
X
2
к
v
U
v
2
3
Т
т
V
V
3
2
Выполнение дифференциального уравнения 1 2 1 2 •
У2 = -е2 + -мы получим У2 = У1 + е2ё2 =
V1 + e2(x2-x2d) = V1 + e2(:^
X ——х —
\т* т*
*2а) = У1 + в2 (~~ез +^~хза - ^Тх2 - х2й) / т* т* '
(22)
Сейчас мы выбираем х3а как виртуальный управляющий сигнал, определяемый выражением (23).
хза = Ф(е1,е2,хы,хы) =^х2+^х2а -
^к2е2-^е1 (23) * *
А* 2 2 А* 1 У '
Подставляя х3а в уравнение (22), получаем:
% = -к1в2 - к2в2 +^*в2вз.
Если в3 = 0, то У2 = —к1е1 - к2е2 <0 и е2 будет сходиться к нулю по экспоненциальному закону
Третий шаг. Путь У3 =\е\_ +1е2 +1ез как функция Ляпунова для подсистемы х3 =
4Б*кд„ 4Е*(кс+С*1) ^ 4Е*А*
" Хл.
Х3 vh *2
Выполнение дифференциального уравнения У3 =-е2 +-е2 + -е2 мы получим; У3 = У2 +
ез{±1- хза) = % + вз(^е4 + ^х4а -V* V*
4М^ХЗ-44М**Х2-±.3<1 ) (24)
V* V*
Сейчас мы выбираем х4а как виртуальный
управляющий сигнал, определяемый
выражением (25)
x4d = $(e1,e2,e3,x1d,x1d,x1d) = —+-х3 +
kq
Ah . Vh . Vh ,
Т^х2 + TT^x3d - тг^кзвз -
AhVh
4Ehkq 3" 4Ehkq
e2 (кз>0) (25)
4т*Е*Кч
Подставляя х4а в уравнение (24), получаем: Уз = -к-е- - к2е22-кзе2 +^езе4
Если в4 = 0, то У3 = -к1е^ - к2е^-к3е^ < 0 и е3 будет сходиться к нулю по экспоненциальному закону
11
Четвертый шаг. Путь У4 = - е-2 + -е2 +
11
-е 32 +-е2 как функция Ляпунова для
1 .
подсистемы х4 =--х4 +—и
Выполнение дифференциального уравнения % = %+ е4(х4 - х4о) = Уз+^4 (-Тх4 +
\ Тр
кТи-хы) (26)
■Ср /
Сейчас мы выберем и является управляющим сигналом, описывается выражением (27)
_ ^ I ьр .
и = —Х4 + —X4d
^^-'Л^ез (27) кр 4 4 Vhkp 3 v '
Подставляя и в уравнение (26), получаем:
У^ = -^е2 - к2е2-к3е2-к4е2 (28) Из выражения (28) находим, что
\?4 = -к^^ - к2в^-кзе
значений к, > 0 (/'=1, 2, 3, 4)
У4 = -к- к2е%-к3е3-к4е4 < 0 для всех
7. МОДЕЛИРОВАНИЕ И РЕЗУЛЬТАТЫ
Параметры гидросистемы приведены в Таблице 1. Это сделано в публикациях [6], [7].
Таблица 1
q
р
р
р
Ном Параметры Обо Количество
1 Масса поршня и нагрузки m 200 (Kg)
2 Площадь поршня Ah 1.47-10"3 (m2)
3 Постоянная времени сервоклапана 0.001
4 Коэффициент усиления сервоклапана kv 3.04 10-3 (m/A)
5 Общий коэффициент утечки Chi 110-n[(m3/i)/Pa]
6 Совокупный коэффициент смоделированных сил демпфирования и вязкого трения Bh 1 -104(N-s/m)
7 Объем камеры цилиндра Vh 2.47*10-4(m3)
8 Эффективный объемный модуль Eh 8108 Pa
9 коэффициент усиления потока/ открытия Kq 2.7(m2/s)
10 усиление потока / давления Kc 1.75*10-11 [(m3/s)/Pa]
11 Длина стабилизированной платформы a 1.4 m
12 Ширина стабилизированной платформы b 0.8 m
13 Начальное расстояние между двумя плоскостью (Длина цилиндров) h 0.66 m
Подбор параметров: K1=5, K=2, K3=4, K4=10.
С выбранными параметрами систему можно моделировать в программном Simulink. Результаты моделирования в настоящее время слишком предварительные, они лишь демонстрируют возможность подавления гармонического возмущения по каждой из координат в отдельности. В дальнейшем планируется моделировать всю систему в целом.
Заключение
В данной работе осуществлено предварительное теоретическое решение стабилизации верхней части платформы, помещенной на палубу корабля. Цель стабилизации состоит в том, чтобы верхняя часть стола оставалась параллельной горизонтальной плоскости в то время, когда палуба корабля подвержена качке.
ЛИТЕРАТУРА
[1] H. E. Merritt, Hydraulic Control Systems. New York: Wiley, 1967
[2] Claude Kaddissi, Jean-Pierre Kenne, and Maarouf Saad, Identification and Real-Time Control of an Electrohydraulic Servo System Based on Nonlinear Backstepping, IEEE/ASME TRANSACTIONS ON MECHATRONICS, VOL. 12, NO. 1, FEBRUARY 2007
[3] Weirui Wang, Shuangxia Pan, Fang Wang, A Comparison Study of Two Disturbance Rejection Control Strategies for Hydraulic Position Servo Systems, the 30th Annual Conference of the IEEE
Industrial Electronics Society, November 2 - 6,2004, Busan, Korea.
[4] Yongsheng Yang, Electro-Hydraulic Servo System with Predictive Control, ICIEA 2009.
[5] K. X. Wang and B. X. Li, Solid rocket propellent and burning, Beijing:
National defense industry Press, 1983, pp. 1-10.
[6] Jianyong Yao, and Wenxiang Deng, Active Disturbance Rejection Adaptive Control of Hydraulic Servo Systems, DOI 10.1109/TIE.2017.2694382, IEEE
[7] Cun Shi, Adaptive decoupling control of dissimilar redundant actuatuon system for large civil aircraft, aerospace Science and technology 47 (2015)
Буй Ван Там - аспирант кафедры автоматики НГТУ. В 2006 году окончил военно-морскую академию в Нячанге. Он получил степень M.S. (магистра) в области кибернетики и автоматизации в военно-технической академии, Ханой, Вьетнам. В 2014 г. С 2015 по 2019 гг. Участвовал в исследованиях в Военно-морской академии Вьетнама. С 2009 года он преподает в Военно-морской академии Вьетнама. Его исследовательские интересы включают математически описательные объекты; обследование качества системы автоматического управления; синтезированный контроллер для объектов на кораблях. E-mail: tamtlpt@gmail .com
Статья получена 20.04.2021.
Construction of a Mathematical Model of a Ship Platform to Suppress Rolling Effects
Bui Van Tam
Novosibirsk State Technical University, Novosibirsk, Russia
Abstract: This article is devoted to the stabilization of the upper part of the mechanical platform located on the deck of the ship. The purpose of the stabilization is to ensure a stable horizontal position of the various measuring instruments. The stabilization system should suppress movements in two degrees of freedom, which generate tilts of the platform in two axes: along and across the ship. The article analyzes the known stabilization methods, makes the choice of the solution method. Also, the article provides justifications for the derivation of a mathematical model of the system and proposes a model in the form of a structural diagram. A method for designing a regulator and preliminary simulation results are also proposed. Key words: stabilization, automation, control, regulator, optimization, dynamic error, static error, modeling.
References
[1] H. E. Merritt, Hydraulic Control Systems. New York: Wiley, 1967
[2] Claude Kaddissi, Jean-Pierre Kenne, and Maarouf Saad, Identification and Real-Time Control of an Electrohydraulic Servo System Based on Nonlinear Backstepping, IEEE/ASME TRANSACTIONS ON MECHATRONICS, VOL. 12, NO. 1, FEBRUARY 2007
[3] Weirui Wang, Shuangxia Pan, Fang Wang, A Comparison Study of Two Disturbance Rejection Control Strategies for Hydraulic Position Servo Systems, the 30th Annual Conference of the IEEE
Industrial Electronics Society, November 2 - 6,2004, Busan, Korea.
[4] Yongsheng Yang, Electro-Hydraulic Servo System with Predictive Control, ICIEA 2009.
[5] K. X. Wang and B. X. Li, Solid rocket propellent and burning, Beijing:
National defense industry Press, 1983, pp. 1-10.
[6] Jianyong Yao, and Wenxiang Deng, Active Disturbance Rejection Adaptive Control of Hydraulic Servo Systems, DOI 10.1109/TIE.2017.2694382, IEEE
[7] Cun Shi, Adaptive decoupling control of dissimilar redundant actuatuon system for large civil aircraft, aerospace Science and technology 47 (2015)
^^^^ Bui Van Tam is PhD-student of g ^ NSTU, Department of Automation.
He graduated university in 2006 at naval academy, Mia Trang. He received M.S. degrees in ^^^ cybernetics and automation engineering from military technical academy, Ha Noi, Viet Nam.
ПРИМЕЧАНИЕ РЕДАКЦИИ
Автор статьи исходит из задачи стабилизации верхней части платформы относительно неподвижной системы координат. Это крайне сложно, и, вероятно, не является лучшей постановкой задачи. Если корабль качается на волнах, то невозможно поддерживать платформу неподвижно относительно дна океана. Ошибку управления в этом случае невозможно измерить. Также в статье не сказано о вертикальной качке, поскольку все рассуждения делаются в терминах углов отклонения от горизонтальной плоскости. Предлагается взглянуть на задачу наилучшего демпфирования всех видов качки. Качку следует рассматривать не в двух координатах, а в шести, как минимум. При этом необходимо учесть, что плавное поступательное движение на даёт качки.
Предлагается рассмотреть для одной точки (например, для центра стола) три вида ускорений: по координатам X, Y, Z. Если корабль движется с ускорением, то платформа также должна двигаться с ускорением, ведь она должна оставаться на корабле. Но самые большие амплитуды ускорения можно подавить активной стабилизацией. Например, если бы платформа была подвешена на магнитной подвеске, то ускорение основания передавалось бы платформе не мгновенно, а с некоторым запаздыванием и демпфированием. Удары по основанию также не передавались бы на платформу. Техническая постановка задачи состоит в том, чтобы ускорения основания платформы передавались бы верхней части, стабилизированной, не полностью, а с некоторой фильтрацией, но при этом интеграл по времени от ускорения, который является скоростью движения, а также второй интеграл по времени от ускорения, который является перемещением, в среднем должен передаваться достаточно точно.
Ожидаемый эффект можно
продемонстрировать следующим примером. Если палуба поднимается и опускается на несколько метров вверх и вниз, не может стоять задача стабилизации платформы на месте, но может стоять задача более плавного подъема и
In 2014. From 2015 to 2019, he was participating in research in Naval Academy, Viet Nam. Since 2009, he has been a teacher in Naval Academy, Viet Nam. His research interests include mathematically descriptive objects; automatic control system quality survey; synthesized controller for objects on ships. E-mail: tamtlpt@gmail .com
The paper has been received on 20/04/2021.
опускания, как если бы тяжелая платформа была подвешена на пружинах.
Упомянутые три координаты относятся только к центру стола. Требуется добавить наклон по продольной и по поперечной осям корабля, это даёт ещё две координаты. И, наконец, имеется обращение корабля вокруг своей оси (или, если угодно, обращение стола вокруг его центра), это даёт шестую координату.
Можно рассуждать и иначе. Плоскость можно задать тремя точками, каждая точка имеет три степени свободы, это даёт девять координат; но фиксированные точки на плоскости образуют треугольник, стороны которого фиксированы, это даёт три связи, что убирает три степени свободы. Девять минус три опять даёт шесть. Это рассуждение эквивалентно, оно приводит к тому же результату.
Таким образом, система обратной связи должна максимально сильно по возможности подавить ускорение, снизить его, но не полностью устранить. Ведь полностью устранить ускорение можно только если платформа окажется неподвижной относительно дна океана. Если платформа, допустим, будет подниматься очень медленно или опускаться очень медленно, и при этом не будет совершать качания, тогда задача будет решена. Любой пружинный или иной подвес не устраняет колебания полностью, он только их делает меньше, ослабляет, но не уничтожает.
Может быть сделана двухступенчатая система. Если грубая система устранит колебания, которые были большие по амплитуде и по частоте, но останутся небольшие остаточные колебания, тогда точная система сможет устранить остаточные колебания.
Поэтому для начала можно поставить задачу хотя бы как-то ослабить колебания. Именно обращать внимание следует на ускорения, а не перемещения.
Следовательно, будут нужны датчики ускорения. Такие существуют, их можно купить, поставить и снять их сигналы. Это во многом очень похоже на управления балансирующим роботом.