Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Березин С.Я. Системы автоматического управления движением судов по курсу

.pdf
Скачиваний:
17
Добавлен:
25.10.2023
Размер:
15.31 Mб
Скачать

Исходя из этого, исследование устойчивости систем автомати­ ческого управления движением судов по курсу целесообразно раз­ бить на два этапа — исследование устойчивости следящей системы управления рулем и исследование устойчивости всей системы.

Исследование устойчивости следящих систем управления рулем.

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

Кроме этих нелинейностей следящие системы управления рулем имеют нелинейность, обусловленную моментом сопротивле­ ния на баллере руля, и нелинейность типа люфта в механических передачах.

Рис. 11.1. Структурная схема следящей системы с учетом нели­ нейности.

ЛЭ] и ЛЭ2 — линейные элементы; НЭ — нелинейный элемент.

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

Точное исследование нелинейных систем автоматического упра­ вления является трудной задачей. Существующие точные методы (метод точечного преобразования, метод Ляпунова) пригодны только для решения некоторых простейших задач [38]. Поэтому для исследования следящих систем управления рулем целесооб­ разно воспользоваться методом гармонического баланса, предло­ женным Н. М. Крыловым и Н. Н. Боголюбовым и развитым Л. С. Гольдфарбом [18] и Е. П. Поповым [26].

Метод гармонического баланса дает незначительную погреш­ ность при исследовании систем автоматического управления дви­ жением судов по курсу благодаря тому, что объект регулирова­ ния (судно) представляет собой хороший низкочастотный фильтр

[26].

Согласно методу гармонической линеаризации [18, 26] следя­ щая система управления рулем состоит из стационарных линейных и безынерционных нелинейных звеньев.

Структурная схема следящей системы управления рулем с уче­ том нелинейности приведена на рис. II. 1. Для исследования сле­ дящих систем управления рулем воспользуемся работами [18, 26]

49

и рассмотрим несколько другой вариант метода Крылова и Бого­ любова, который удобнее для практических исследований автору­

левых.

Предположим, что управляющее воздействие ср(0 отсутствует, а система находится на границе устойчивости и в ней на входе нелинейного элемента возникли незатухающие колебания неко­ торой частоты со и амплитуды а. Тогда уравнения первого прибли­ жения, согласно [18, 26], запишутся следующим образом:

X i =

Л-4&0. с»

 

хя =

х1\Рл1(/'ш);

( П Л )

x3

= x2W„(a, со);

 

х4

= х ^ л2{}а).

 

Откуда уравнение свободных колебаний системы будет иметь вид

l + WH(a, со)Гл1( Н Г л2(/со)/е0.с- 0 .

(II.2)

Обозначим второе слагаемое в левой части уравнения (II.2) через W3K3{a, /со) и назовем его эквивалентной комплексной пере­ даточной функцией нелинейной следящей системы. Обратную ве­ личину обозначим через УЭКв(а, /со) и назовем обратной эквива­ лентной комплексной передаточной функцией системы:

Y 3(a, /со) =

1 . = Кл1 (/со) Y л2 {ja)YH(а, со).

(II.3)

W3 (a, /со)

 

При a = const можно

построить эквивалентную обратную ам­

плитудно-фазовую характеристику разомкнутой системы УЭкв(а, /со). Семейство таких характеристик, построенных для различных зна­ чений а, изображено на рис. II.2. Каждому значению а будет со­ ответствовать определенная эквивалентная обратная амплитудно­

фазовая характеристика.

При а = ап кривая пройдет через точку

(-1-JO).

устойчивости

периодического

решения

Для исследования

с амплитудой а = ап дадим малое приращение Аа

амплитуде

ап.

При а = ап.+ Аа характеристика займет

положение

I или

II

(см.

рис. II.2). Из линейной теории известно, что кривая I, охватываю­ щая точку (—/; /О), соответствует затухающим колебаниям пере­ ходного процесса, а кривая II — расходящимся колебаниям. По­ этому, если при увеличении амплитуды (Да>0) эквивалентная обратная амплитудно-фазовая характеристика будет охватывать точку (—1; /О), а при уменьшении амплитуды (Аа<0) не будет ее охватывать, то переходный процесс в системе будет таким, что колебания с амплитудой, большей, чем ап, затухают, а колебания с амплитудой, меньшей, чем а„, расходятся. Следовательно, пере­ ходный процесс с обеих сторон сходится к исследуемому периоди­ ческому процессу с амплитудой ап. Это означает устойчивость про­ цесса, т. е. в системе имеют место автоколебания.

50

Если лее при увеличении амплитуды (Да>0) характеристика (кривая II) не охватывает точку (—1\ jO), а при Аа < а — охва­ тывает (кривая /), то переходный процесс в обе стороны расхо­ дится, т. е. исследуемое периодическое решение неустойчиво (си­ стема устойчива в малом и неустойчива в большом).

Значения а и со, соответствующие предельному циклу, можно найти при модуле обратной амплитудно-фазовой характеристики, равном единице, и при ее фазе ф=180°, т. е. определяются точкой

(-1-JO).

Таким образом, по расположению эквивалентной обратной амплитудно-фазовой характеристики на комплексной плоскости можно судить об устойчивости автоколебаний системы.*

Рис. 11.2. Эквивалентные обратные

Рис. П.З. Характеристика

нелинейно-

амплитудно-фазовые характеристики.

сти вида насыщения.

Обратная передаточная функция замкнутой следящей системы

управления рулем имеет вид (см. рис. II.1)

 

Х с(а, p) = k0,c+ Y Jll(p)YHlYM{p)YM(a)Yia(a),

(П.4)

где yHi(a), Упг(а), Унз(а)— обратные передаточные функции нели­ нейных элементов соответственно от нелинейности вида насыще­ ния, от момента сопротивления на баллере руля, от люфта в меха­ нической передаче.

В следящей системе управления рулем в функцию УЛ1(Р) вхо­ дит коэффициент усиления чувствительного элемента и передаточ­ ные функции усилителей:

Улг(Р)

*i (Р)

к (р)

*2 (Р)

 

 

 

где k(p) — многочлен от р со свободным членом;

ki — коэффициент усиления

звеньев, входящих в линейный

элемент.

 

 

* Условие устойчивости является лишь необходимым, но не достаточным.

51

Аналогично для обратной передаточной функции исполнитель­ ного устройства

Yni{P) = ^ p r = ~rPfh(p),

(П. 5)

*4 (Р) k2

 

где k2— коэффициент усиления по скорости исполнительного дви­ гателя; k2(p) — многочлен от р со свободным членом.

Рассмотрим основные нелинейности, присущие следящим си­ стемам управления рулем.

Н е л и н е й н о с т ь в и д а н а с ы щ е н и я . В усилительных уст­ ройствах следящего рулевого привода имеет место нелинейность такого вида. Наиболее часто эта нелинейность проявляется в око­ нечном усилителе (электромашинном или гидравлическом), в ко­ тором насыщение наступает, как правило, при значительных от­ клонениях регулируемой величины. Предварительные усилители конструируются обычно с большим запасом по напряжению, и на­ сыщение их наступает гораздо позднее.

Эквивалентная передаточная функция нелинейности вида на­

сыщения известна [9] и имеет вид

 

 

q (а) = — ky

arcsin----

(II.6)

 

 

akу

 

где ky

(рис. II.3);

а —- амплитуда колебаний на

входе усили­

 

в

 

 

теля.

При а<В q(a)—ky — линейная характеристика.

Для удобства дальнейших исследований перейдем от выраже­ ния (II.1) к уравнению

*7(а)

2

arcsin-----1-----

(II-?)

k y

Л

А

А

 

где А = — , п построим графическую зависимость

НА).

В

 

 

 

Я (о)

На рис. П.З приведена зависимость

= f(A) = Y nl(A).

Я (а)

Пользуясь этим рисунком, можно определить значение Y-al(A) при заданных значениях ky и В для выбранной амплитуды гармониче­ ских колебаний на входе нелинейного элемента с нелинейностью вида насыщения.

Из рис. П.З нетрудно установить, что абсолютная величина

q (а) = — ^ — всегда меньше ky и только при У„1(А) = 1 ky=\q(a)\.

*HI (Л)

Поэтому, если обозначить через k0 коэффициент усиления следя­ щей системы при отсутствии насыщения k0 = k0iky и через ko(а) — при наличии насыщения, то \ko(a) \<iko при всех значениях ампли­ туд, если рассматриваемая нелинейность имеет место.

52

ко. с
При заданных значениях частоты и амплитуды на входе си­ стемы можно определить величину амплитуды а на входе нели­ нейного элемента. Из рис. II.1, где амплитуда а на входе нелиней­ ного элемента обозначена х2, видно, что
У Л 2 (/< * > )

Таким образом, наличие в основном контуре следящей системы управления рулем нелинейностей вида насыщения будет эквива­ лентно уменьшению коэффициента усиления системы и влияние их на устойчивость будет благоприятно. Качество системы при этом будет ухудшаться. Поэтому при исследовании следящей си­ стемы управления рулем необходимо учитывать возможное ухуд­ шение качества регулирования. Так, например, если известно при­ мерное значение амплитуды на нелинейном элементе, то, поль­ зуясь кривой на рис. II.3, можно определить для данной нелиней­ ности пределы уменьшения коэффициента усиления системы.

Для того чтобы установить амплитуду и частоту колебаний на входе нелинейного элемента, необходимо воспользоваться ампли­ тудно-частотной характеристикой системы управления движением судна по курсу и результатами натурных испытаний авторулевых. По амплитудно-частотной характеристике определяется рабочая частота системы, а по результатам натурных испытаний — значе­ ние амплитуды рыскания. Используя опыт эксплуатации автору­ левого «Самшит», можно считать, что в режиме стабилизации при наиболее неблагоприятных условиях величина рыскания колеб­ лется в пределах от 2,5 до 5,5° курса.

При следящем управлении курсом амплитуда управляющего воздействия должна быть задана.

Если предположить, что в системе автоматического управления

движением судна по курсу под влиянием управляющего или воз­ мущающего воздействия установились колебания с некоторой частотой ыо и амплитудой аоюах на входе системы,* то условное уравнение движения следящей системы управления рулем примет вид

1+ У Л 1 ( / « > ) У Н1 ( д ) р = ф (II.8>

*1 _

_________ Ц^'л! (/<■))_________ _

а

1+*о.сГЛ1 (/со) W (а) Г л* (/«*>)

 

(П.9).

 

Ул! (/ш) + k0

 

h i (а) Улг ( / ш )

Изображенная на рис. П.З зависимость kyYHi(A) =f(A) явля­ ется прямой линией, проведенной под углом ф к оси абсцисс. Поэтому для рассматриваемой нелинейности не представляет труда выразить Уш(а) через а, т. е.

KHi(a) =

a tg ф

(НЛО).

 

kvb

 

* Полагаем, что в системе свободные колебания отсутствуют (нормальный режим следящей системы).

53-

Подставляя это выражение в уравнение (П.9), получим

*2

1

(II.11)

а

k0. с ■ky * Ь

 

Кл1(/ш) +

Ут (/co)atg<p

 

Обратные передаточные функции УЛ1(/со) и Ул2(/со) обычно известны из статического расчета следящей системы управления рулем.

Из выражения (11.11) можно записать

«О шах — a

а

 

k0. Q_kyb

 

УЛ2 О®) « tg ф

 

откуда

 

 

шах_________f e 0 , ckyb_______________

( 11. 12)

 

 

Y*i (/со)

УЛ1 (/®) УЛ2 (/to) tg ф

 

 

При графическом построении функций УЛ1(/со) и УЛ2(/со) на

комплексной плоскости (рис. 11.4) получим

 

д __Цр ш ах__________Й0, ck yb______

 

 

“ I оА |

| ол 11 OB I tg ф

 

Полученное выражение для амплитуды колебаний на входе

нелинейности нельзя

считать окончательным, так как

при выводе

 

 

его не учитывались имеющиеся в си­

 

 

стеме корректирующие

устройства.

 

 

Поэтому величина а в процессе

 

 

дальнейшего расчета должна уточ­

 

 

няться.

 

о б у с л о в ­

 

 

Не л и н е й н о с т ь ,

 

 

л е н н а я

м о м е н т о м с о п р о ­

 

т и в л е н и я

на

б а л л е р е

 

руля.

Величина и характер мо­

 

мента на баллере руля, а значит, и

 

на исполнительном рулевом устрой­

 

стве зависят от разнообразных фак­

Рис. 11.4. К определению ампли­

торов:

от

влияния

конструкции

туды колебаний на входе нели­

руля,

корпуса судна и рулевого уст­

нейного элемента.

ройства, угла и скорости переклад­

 

ки руля, скорости судна [25].

В существующей литературе [25]

показано, что для балансир­

ных рулей, применяемых в настоящее время на судах, зависимость

МС= /(Р) с достаточной степенью точности может

быть представ­

лена в виде отрезков прямых (рис.

II.5). Уравнение Мс для кри­

вых, представленных на рис. II.5, а, запишется в виде

МС= М 0 при P < d

и при — < 0 ;

(11.13)

 

dt

Mc = Mo+ (P— d)tga

при d < p < d + p

И

54

Если входная величина (отклонение

руля)

изменяется по

си­

нусоидальному закону p= asinco^

(рис.

II.5,

б), то в точках

В,

D имеем

 

 

 

 

arcsin — = arcsin —-— ;

 

 

a

d л . ь

(11.14)

U« ■=arcsin а —

,

 

 

 

 

где а — амплитуда гармонического сигнала.

а)

В)

Проведя гармоническую линеаризацию характеристики момента сопротивления на исполнительном рулевом двигателе и исполь­ зуя [9], получим

Мс = д ( а ) + ^ р ,

 

где

2Мп 7 ^ 2 1 /! + - ^ -

 

Ч (а) = ла

 

 

4

4

 

d

 

 

 

 

I —

 

 

d2

 

а

a tg a— у

(11.15)

 

1— — d t g a

Если угол перекладки руля больше d, то

 

 

 

 

 

< 1 1 л 6 >

?,(<■ )= —

l-$ - +

- £ - i ) t t a (a > d ).

(11.17)

па

\ 2

2а

]

 

Если угол перекладки руля меньше d, то

Qi (а) = 0.

55

Для удобства определения значений q(a) и 91(a) построим их графики в функции от амплитуды. Разобьем q(a) на два слагае­

мых

q(a) = q(a) + q*(a),

(II.18)

где

 

9(a) = — ° ;

( И - 19)

ла

 

Рис. II.6.

График нелинейной

функ-

Рис. 11.7. График

нелинейных функций

 

 

 

 

 

 

 

Кривая / —

 

( я )

(

f

( а \

 

 

 

ции Л)о

 

 

 

------- =

 

— .

 

 

 

 

 

Крив»

I I -

tc а

 

1

Vd )

 

 

 

 

 

 

 

 

^

=

 

( - ) .

 

 

 

 

 

 

 

 

 

 

tg а

 

 

l <1)

 

Ha

рис. 11.6 приведена зависимость

^-^ = /(a),

пользуясь кото-

рой, при

известном

значении

 

М„

найти

9 (a).

 

 

М0 нетрудно

 

 

 

 

 

Т а б л и ц а

II .1

 

 

 

 

 

 

 

 

 

 

 

 

Результаты

вычисления Ч * (а ) и

9i (а)

 

 

 

 

 

 

 

 

 

 

 

 

tg я

tg я

 

 

 

 

 

 

 

а

 

1

 

1,5

2.0

2,5

3,0

3,5

 

 

4.0

5,0

 

 

 

 

 

 

Я* (а)

0

 

0,122

0,193

0,252

0,292

0,322

 

0,343

0,375

tga

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<7i

(а)

0

 

0,0356

0,0795

0,114

0,141

0,162

 

0,179

0,204

tga

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

56

На рис. II.7 приведена зависимость Я* (а) = л — (кривая /),

вычисленная по выражению

 

 

tg а

' \ d

 

(11.20). На том же рисунке построена

зависимость Ях (а)

«/ Д \ _

2 ( 1 I

< Р ___d_

(кривая //). Ре-

t g a

\ d /

я \ 2

2a2

a

 

 

зультаты вычисления обеих кривых сведены в табл. II.1.

рулем

При проектировании

следящей

системы

управления

обычно известна зависимость Mc = f(P), т. е известны М0, d

и tga,

поэтому, пользуясь кривыми I и II

(см. рис.

II.7), легко

найти

для заданного значения

амплитуды

величину

q(a) и qi(a).

 

tga

Рис.

11.8. — Амплитудно-фазовая

характери­

стика нелинейности от момента

сопротивле­

 

ния на баллере руля.

 

На рис. П.8

приведена амплитудно-фазовая характеристика

q(a) +jqi(a), построенная по кривым рис. II.7. Здесь же показано построение вектора OB — q(a)+jql(a) для а = 30°, ^ = 2, а = 20°,

М0= 10 кГм:

0д = я* (а) + /<?1 (о) tga

ОС = ?* (a) + f a (а);

CB = ?(a) = i^ L ;

я а

OB = q (а) -\-jq1 (а) = q* (а) + q~а) + fa (а).

В частных случаях, когда МС= М0 или d = 0, определение не­

линейной функции значительно упрощается.

При Мс=Мо

<7i(a)=0 и q(a) = ^ - .

па

57

При d = О

Ч(а) =

<h (а)

4М 0 |

t g a .

ла

2

tg a

( 11.21)

 

Путем несложных преобразований можно установить, что на­ личие момента сопротивления на исполнительном двигателе дей­ ствует благоприятно на устойчивость следящей системы управле­ ния рулем. Покажем это на примере электрической следящей сис­ темы управления рулем.

Уравнение движения привода при учете Мс равно

J\/[ =

d(°R

■М'

(П.22)

д ~

375

dt

 

 

где Мд — момент, развиваемый двигателем

(в случае электриче­

ской следящей системы MR = cJ)\

М 'с — приведенный момент на­

грузки М 'с = Мс&мп; GD2— маховый

момент

двигателя с учетом

приведенных масс рулевого устройства, в

кГ-м2; сод — угловая

скорость двигателя.

 

 

 

 

Для электрической следящей системы управления рулем из вы­

ражения Мд нетрудно получить

 

 

 

 

 

ОР2Р2Мд

[*<«>+•*#■

 

(II.23)

 

L

 

- * 2м .Л .

 

 

375йдС[

 

 

 

 

где йч — угол

поворота

вала

двигателя; kMп = —^ — коэффициент

механической передачи;

.

Ел .

,. Еа

шд

 

кд= — ,

/гд= '^ '-

 

 

Уравнение равновесия э. д. с. в главной цепи двигателя имеет

вид

 

Er = IR + Lpl -f- йдрРд.

 

(II.24)

 

 

 

Подставляя

в это выраженгние значение / д из уравнения (11.23),

получим

 

 

q(a) I

41(«)

 

 

 

 

 

 

 

Е = BpVtnR+BTkp>P

 

 

kl.n(R + L) Рд+ ^дРРд’

откуда эквивалентная обратная передаточная

функция

двигателя

с учетом Мс равна

 

 

 

 

 

 

Уэкв.д(а, Р) = ^

= ^-(ч(а) + ^

р

) х

 

Рд

^2

 

 

 

 

х mkl п (\ + Тр) + р(1 + Вр + ВТр>),

где

 

L .

 

R

 

 

= — ; т -

R

т

= ---------------- .

 

 

 

— :

 

При Т= 0, что часто имеет место,

 

 

 

 

^ЭКВ. д(я> р) — .

д(а)

Qi(a)

 

mk2

-

p(l + Bp).

 

 

 

 

 

М .

П

 

(II.25)

(П.26)

58

Соседние файлы в папке книги из ГПНТБ