ОТУ_Поляков_Часть2_2009. Теория автоматического управления для чайников Часть ii. Управление при случайных
 Скачать 1.41 Mb.
|
|
4. Морское волнение 4.1. Что такое морское волнение? Морское волнение – это колебания поверхности воды, вызванные ветром (а также прили- вами, отливами и другими причинами). Простейшая модель морского волнения – гармониче- ские колебания поверхности, когда волновая ордината (вертикальная координата точек по- верхности) изменяется по закону синуса. Такое волнение называют регулярным (синяя штри- ховая линия на рисунке). Однако, в самом деле волновая ордината меняется по более сложному закону, его можно (опять-таки приближенно) представить как сумму большого количества гармоник (синусоид) с разными амплитудами и фазами (теоретически – это сумма бесконечного числа гармоник). Это так называемое нерегулярное волнение, которое чаще всего описывается как случайный про- цесс. Можно предположить, случайные волновые ординаты распределены по нормальному за- кону, и это действительно подтверждается экспериментами. При этом плотность распределения высот волн подчиняется закону Рэлея. В нашей стране для характеристики волнения используется высота волны 3%-ной обеспе- ченности % 3 h – это высота волны, вероятность превышения которой составляет 3%. Значения % 3 h определяются по таблице в зависимости от интенсивности волнения (в баллах): Волнение, баллы % 3 h , м Волнение, баллы % 3 h , м Волнение, баллы % 3 h , м 1 0 – 0,25 4 1,25 – 2,0 7 6,0 – 8,5 2 0,25 – 0,75 5 2,0 – 3,5 8 8,5 – 11 3 0,75 – 1,25 6 3,5 – 6,0 9 более 11 По значению % 3 h можно определить дисперсию волновой ординаты r D . Известно, что вы- сота волны h – это случайная величина, распределенная по закону Рэлея: 0 , 8 exp 4 ) ( 2 ≥ ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − = h D h D h h f r r Вероятность превышения некоторого уровня % 3 h можно рассчитать как интеграл от плотности вероятности на интервале ) ; [ % 3 ∞ h : 0 h ) (h f % 3 h 0 t 1 r © К.Ю. Поляков, 2009 32 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − = ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − = ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − = > ∞ ∞ ∫ r h r h r r D h D h dh D h D h h h p 8 exp 8 exp 8 exp 4 ) ( 2 % 3 2 2 % 3 % 3 % 3 Эта вероятность должна быть равна 3% или 0,03; поэтому 2 % 3 2 % 3 2 % 3 0356 , 0 03 , 0 ln 8 03 , 0 8 exp h h D D h r r ≈ − = ⇒ = ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − Значение % 3 h определяет также средний период волнения T , среднюю частоту волнения ω и частоту максимума спектра m ω : % 3 1 , 3 h T ≈ , T π ω 2 = , ω ω 71 , 0 ≈ m Существуют два типа моделей нере- гулярного волнения, двухмерная и трех- мерная. В двухмерной модели предполага- ется, что гребни волн имеют бесконечную длину и перемещаются параллельно друг другу в одном направлении (рисунок слева). Каждый, кто видел реаль- ное волнение, знает, что на практике это не так. Более точная (но и более сложная) трехмерная модель учитывает сложение множества двухмер- ных волн, идущих в разных направлениях (рисунок справа). Строго говоря, волнение – это нестационарный процесс. Оно начинается с ряби, затем ве- тер за счет завихрений потоков воздуха разгоняет волны. Но если ветер с постоянными свойст- вами действует достаточно длительно время (несколько часов) на большой акватории (десятки километров), можно говорить о развитом (установившемся) волнении, которое считают ста- ционарным эргодическим процессом. Такой подход позволяет использовать методы анализа случайных процессов на основе корреляционных функций и спектральных плотностей. Под термином «спектр морского волнения» обычно понимают спектральную плотность волновой ординаты. С экспериментальными данными лучше всего согласуются экспоненци- альные спектры вида ( ) n k r B A S − − − = ω ω ω exp ) ( , где параметры A , B характеризуют интенсивность волнения, а k и n зависят от его особенно- стей. Вообще говоря, спектры волнения в разных районах от- личаются. Если информации недостаточно, рекомендуют ис- пользовать типовой спектр Пирсона и Мошковица ( 5 = k и 4 = n ), рекомендованный 12-й Международной конференци- ей опытовых бассейнов (МКОБ): ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ − ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ = 4 5 25 , 1 exp 06 , 7 ) ( ω ω ω ω ω π ω m m r r D S Известно множество других выражений для спектров. Часто используют, например, спектр Неймана: 0 ω , рад/c ) ( ω r S , м 2 /Гц 0,5 1 2 3 спектр МКОБ спектр Неймана © К.Ю. Поляков, 2009 33 ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − = 8 , 0 % 3 2 6 575 , 4 exp 2 , 1 2 ) ( h S r ω ω π ω На рисунке справа показаны графики спектра МКОБ (синяя линия) и спектра Неймана (зеленая линия) для волнения 4 балла ( 2 % 3 = h м). Экспоненциальные спектры хорошо согласуются с экспериментальными данными, однако неудобны для моделирования и расчета систем управления. Вместо них используются дробно- рациональные спектры 1: 2 2 2 2 2 2 4 2 2 2 ) ( ) ( 2 ) ( 2 ) ( β α ω β α ω ω β α α ω + + − + + + = r r D S , 2: 2 2 2 2 2 2 4 2 2 ) ( ) ( 2 ) ( 4 ) ( β α ω β α ω β α α ω + + − + + = r r D S , 3: 2 2 2 2 2 2 4 2 ) ( ) ( 2 4 ) ( β α ω β α ω αω ω + + − + = r r D S Здесь r D – дисперсия волновой ординаты, α – коэф- фициент затухания и β – угловая частота корреляционных функций. Для развитого морского волнения m ω β 02 , 1 = и β α 21 , 0 = На рисунке справа показаны спектр МКОБ и дробно-рациональные спектры 1 и 3 (кривая для спектра 2 проходит между ними). Графики дробно-рациональных спектров проходят выше экспоненциальных на низких частотах, которые больше всего влияют на поведение инерцион- ных морских объектов (судов). Поэтому считается, что использование дробно-рациональных спектров приводит к несколько завышенным оценкам параметров движения (качки, рыскания). Заметим, что спектр 3 равен нулю на нулевой частоте, так же, как и экспоненциальные спектры. Однако, если учитывать трехмерность волнения, по современным представлениям спектр на нулевой частоте в большинстве случаев отличен от нуля. Это оправдывает использо- вание дробно-рациональных спектров 1 и 2. 4.2. Кажущиеся спектры До этого мы рассматривали спектры волнения без учета движения судна. Очевидно, что двигаться против волнения труднее, чем «по волне», а движение судна «лагом» к волне (при бортовой волне) вызывает сильную качку. Это означает, что эффект действия волнения на суд- но зависит от скорости и направления его движения по отношению к основному направлению распространения волн. Математически это выражается в изменении спектральной плотности возмущения. В физике хорошо известен эффект Доплера, который состоит в том, что при движении датчика измеренная им частота волны изменяется в зависимости от его собственной скорости по закону ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ + ⋅ = u v k 1 ω ω , где ω – истинная частота волны, k ω – кажущаяся (измеренная) час- тота, v – скорость движения датчика в направлении источника волн, а u – скорость распро- странения самих волн. При движении судна против волнения ( 0 > v ) кажущаяся частота волн будет больше, чем истинная, а при движении «по волнению» ( 0 < v ) – меньше истинной. 0 ω , рад/c ) ( ω r S , м 2 /Гц 0,5 1 2 3 спектр МКОБ спектр 1 спектр 3 © К.Ю. Поляков, 2009 34 Из гидродинамики известно, что скорость распространения волны с частотой ω (на глу- бокой воде) равна ω g u = . Кроме того, если судно движется со скоростью V (в м/с) по углом ξ к направлению распространения волн (считается, что 0 = ξ соответствует движению против вол- ны), оно приближается к источнику со скоростью ξ cos V v = . Таким образом, формула преоб- разования истинной частоты в кажущуюся приобретает вид ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ + ⋅ = ω ξ ω ω g V k cos 1 Если судно «убегает» от волн, то при увеличении скорости кажущаяся частота, вычислен- ная по этой формуле, оказывается отрицательной – это означает, что судно обгоняет волны. По- этому для того, чтобы работать только с положительными частотами, нужно взять модуль вы- ражения в правой части. Окончательно получаем ( ) ω ω ω a k + ⋅ = 1 , (8) где g V a ξ cos = – фактор относительного движения. Итак, составляющую волнения с частотой ω судно воспринимает с кажущейся частотой k ω Чтобы построить кажущийся спектр, нужно научиться решать обратную задачу – по за- данной кажущейся частоте k ω определить частоту (или частоты!) исходного спектра, которые судно воспринимает как k ω . Используя (8), формально получаем два квадратных уравнения 0 2 = + + k a ω ω ω и 0 2 = − + k a ω ω ω , (9) Нас интересуют все вещественные и положительные решения этих уравнений (фактиче- ски таких «подходящих» частот может быть от одной до трех). Предположим, что судно идет под острым углом к волне, так что 0 cos > ξ и 0 > a . В этом случае первое уравнение в (9) заведомо не имеет подходящих решений (поскольку 0 2 ≥ + ω ω a ), а второе имеет только одно положительное решение. На рисунке слева изображен график функции ω ω ω + = 2 a k , который ясно показывает, что одной кажущейся частоте 0 1 > k ω соответствует только одна частота 1 ω исходного спектра и наоборот. Теперь предположим, что судно идет под тупым углом к волне («убегает» от волн), так что 0 cos < ξ и 0 < a . На рисунке справа показан график функции ω ω ω + = 2 a k для этого слу- чая. Используя знания школьной математики, легко увидеть, что парабола пересекает ось 0 = k ω при a 1 0 − = = ω ω , а ее вершина находится в точке ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ 4 ; 2 0 0 ω ω . Все волны с частотами больше 0 ω судно опережает, поэтому кажущаяся частота становится отрицательной. 0 ω k ω 0 < a 0 ω 2 / 0 ω 4 / 0 ω 1 ω 2 ω 3 ω 1 k ω 1 k ω − 0 ω k ω 0 > a 1 ω 1 k ω © К.Ю. Поляков, 2009 35 Если интересующая нас кажущаяся частота 1 k ω меньше максимальной (равной 4 / 0 ω ), то энергия волн на трех частотах ( 1 ω , 2 ω и 3 ω ) складывается на частоте 1 k ω кажущегося спектра. Если 4 / 0 1 ω ω > k , остается только оно решение – 3 ω Нужно понимать, что мощность волнения (средний квадрат волновой ординаты) не зави- сит от скорости движения судна. Для движущегося объекта эта мощность просто перераспреде- ляется по частотам. Поэтому при расчете кажущегося спектра важно сохранить мощность. Это значит, что средние квадраты волновой ординаты для исходного и кажущегося спектров долж- ны быть равны: ∫ ∫ ∞ ∞ = 0 0 ) ( 1 ) ( 1 k k rk r d S d S ω ω π ω ω π Чтобы добиться этого, приравнивают «элементарные» мощности: k k rk r d S d S ω ω ω ω ) ( ) ( = , что дает k r k rk d d S S ω ω ω ω ) ( ) ( = Выражение k d d ω ω / в правой части последнего равенства можно найти как обратное значение для производной (при 0 > a ) ( ) 1 2 2 + = + = ω ω ω ω ω ω a a d d d d k , так что 1 2 1 + = ω ω ω a d d k При 0 < a нужно взять эту производную по модулю (чтобы все составляющие спектра склады- вались, а не вычитались). В самом сложном случае, когда на частоте 1 k ω кажущегося спектра складываются три составляющие исходного спектра (соответствующие частотам 1 ω , 2 ω и 3 ω ), при численном пересчете спектра имеем: 3 3 2 2 1 1 1 2 1 ) ( 2 1 ) ( 2 1 ) ( ) ( ω ω ω ω ω ω ω a S a S a S S r r r k rk + + + + + = На рисунке показаны спектры, полученные преобра- зованием спектра МКОБ (Пирсона и Мошковица) при разных углах встречи с волной (волнение 4 балла, 2 % 3 = h м, скорость судна 4 м/с). При встречном волнении (зеленая линия) максимум спектра смещается в область высоких частот, спектр «размазывается» по оси ω . Учитывая, что инерционный объект подавляет высокочастотные возмущения, можно считать, что влияние волнения на судно уменьшается. При попутном волнении (красная линия) спектр имеет разрыв (скачок), такие разрывы подтверждаются экспериментально. Вся энергия волн сосредоточена в уз- кой полосе на низких частотах. При этом возмущение оказывает очень сильное влияние на суд- но, и это влияние сложно (а иногда и невозможно) скомпенсировать с помощью управления. Поэтому движение на попутном волнении – это один из самых сложных режимов движения в точки зрения стабилизации курса. 0 ω , рад/c ) ( ω r S , м 2 /Гц 3 1 2 0 = ξ o 45 = ξ o 135 = ξ 3 2 1 © К.Ю. Поляков, 2009 36 Пересчет дробно-рациональных спектров выполня- ется достаточно просто. Вычисляют частоту максимума кажущегося спектра 2 m m mk a ω ω ω + = , и в приведенных выражениях для спектров используют mk k ω β 02 , 1 = и k k β α 21 , 0 = вместо α и β . Дисперсия r D сохраняется та же самая. На рисунке слева показаны дробно- рациональные спектры (типа 1) для разных углов встречи с волной (волнение 4 балла, 2 % 3 = h м, скорость судна 4 м/с). Видно, что общий характер смещения спектров та- кой же, как и для пересчитанного спектра МКОБ. 4.3. Моделирование действия морского волнения на судно Спектры, рассмотренные выше, относились к изменению волновой ординаты, тогда как для разработчика систем управления важно выяснить влияние волнения на динамику судна (угол рыскания, качку и т.п.). Как связано изменение волновой ординаты и рыскание (качка) судна? К сожалению, ответ на этот вопрос не прост и требует довольно полной информации о судне, например, об обводах корпуса. В конечном счете, (для линейной модели) требуется по- строить частотную характеристику ) ( ω j H , которая связывает волновую ординату и требуемую характеристику движения судна, например, угол рыскания или угол крена. Рассмотрим более подробно задачу стабилизации судна по курсу. Простейшая линейная модель судна – это модель первого порядка (модель Номото), которая связывает угловую ско- рость вращения вокруг вертикальной оси (оси Z) и угол перекладки руля: δ ω ω S z S z T K T + − = 1 & , где z ω – угловая скорость (в рад/с), δ – угол поворота руля (в радианах), S T – постоянная времени (в секундах) и K – безразмерный коэффициент. Угловая скорость (в такой простейшей модели) равна производной от угла рыскания ϕ (так называется угол отклонения от заданного курса). Поэтому к модели нуж- но добавить еще одно уравнение: z ω ϕ = & Обычно угол рыскания считается положительным при вращении против часовой стрелки. Волновое возмущение приводит к тому, что появляется дополнительное вращение, откло- няющее судно от заданного курса, то есть, уравнение Номото с учетом возмущения принимает вид w T K T S z S z + + − = δ ω ω 1 & , (10) где w – возмущение, вызванное морским волнением. Будем считать, что все случайные процес- сы в системе – центрированные, то есть, их математические ожидания равны нулю. Спектр ) ( ω w S возмущения w рассчитывается по общей формуле ) ( ) ( ) ( ) ( ω ω ω ω r w S j H j H S − = , 0 ω , рад/c ) ( ω r S , м 2 /Гц 1 1 2 0 = ξ o 45 = ξ o 135 = ξ 3 0,5 заданный курс ϕ z ω © К.Ю. Поляков, 2009 37 где ) ( ω r S – спектральная плотность волновой ординаты, а ) ( ω j H – частотная характеристика, описывающая реакцию судна на волнение. Для построения ) ( ω j H нужно использовать уравне- ния движения конкретного судна. Спектральную плотность ) ( ω w S часто удается достаточно успешно аппроксимировать дробно-рациональным спектром вида 4 0 2 2 2 0 4 2 2 ) 1 ( 2 ) ( ˆ ω ω λ ω ω ω ω + − + = w w K S , где 0 ω – доминирующая частота волн, λ – коэффициент демпфирования ( 1 0 < < λ ), w w K σ λω 0 2 = и w σ – коэффициент, определяющий интенсивность волнения. Формирующий фильтр для такой спектральной плотности имеет вид 0 0 2 2 ) ( ω λω + + = s s s K s F w Легко проверить, что ) ( ) ( ) ( ω ω ω j F j F S w − = . Таким образом, для моделирования волне- ния нужно пропустить белый шум через формирующий фильтр с передаточной функцией ) (s F и полученный сигнал ) (t w подключить к модели Номото (10). Соответствующая блок-схема системы выглядит так: Эта схема описывает разомкнутую систему (без обратной связи). Система управления с обратной связью включает также регулятор, привод руля, измерительную систему (датчики): Точность стабилизации на курсе определяется среднеквадратическим отклонением (СКВО) ϕ σ угла рыскания (или его дисперсией). Кроме того, требуется, чтобы мощность управления (СКВО угла перекладки руля δ σ ) не была слишком большой. Как правило, эти за- дачи противоречивы. Это значит, что для уменьшения рыскания нужно увеличивать мощность управления, а при уменьшении активности руля увеличивается рыскание. Кроме внешних возмущений, нужно обязательно учитывать шум (помехи) измерений (сигнал m на схеме). Дело в том, что неправильно спроектированный регулятор может хорошо работать в идеальных условиях, но приводить к «раскачиванию» системы при малейших ошиб- ках измерений. w δ z ω 1 + s T K S s 1 ) (s F w белый шум привод регулятор u датчики + – ε 0 ϕ ошибка заданный курс шум измерений m ϕ w δ z ω 1 + s T K S s 1 ) (s F w белый шум ϕ © К.Ю. Поляков, 2009 38 Шум измерений обычно моделируется как белый шум, имеющий равномерный спектр на всех частотах. Поскольку объект (апериодическое звено) является фильтром низких частот, сигнал ϕ (изменение курса) не может быть высокочастотным. Поэтому регулятор должен по- давлять высокочастотные сигналы, которые явно вызваны ошибками измерений. |