книги из ГПНТБ / Березин С.Я. Системы автоматического управления движением судов по курсу
.pdfИсходя из этого, исследование устойчивости систем автомати ческого управления движением судов по курсу целесообразно раз бить на два этапа — исследование устойчивости следящей системы управления рулем и исследование устойчивости всей системы.
Исследование устойчивости следящих систем управления рулем.
Следящие системы управления рулем являются нелинейными. Эти системы могут быть релейными и непрерывного действия. В пер вых имеется нелинейность типа реле, а во вторых, как правило, нелинейность типа насыщения, так как мощность исполнительных двигателей в таких системах ограничена.
Кроме этих нелинейностей следящие системы управления рулем имеют нелинейность, обусловленную моментом сопротивле ния на баллере руля, и нелинейность типа люфта в механических передачах.
Рис. 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.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 имеем |
|
|
|
|
Uл —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