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

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

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

пульсной переходной функции р(^) при />0 определяют методом трапецеидальных частотных характеристик по выражениям [28]

P(0= Pcos(*)+ Ps,n(*);

рr COS

(/) =— 2

мл,ю-

sin cot-/ sin Д^

 

 

 

 

W

 

д

 

01

I

a it

 

 

 

 

 

 

 

 

 

 

 

1 v

Voi

 

cos u>it sinAj^

 

p.,„(0 =

—-

i=l

 

(Oit A^

<s>it

n

 

 

(in .70)

(111.71)

(111.72)

где UotJ Voi — высота трапеции.

 

 

III.42).

Значения Д* и ш,

зависят от наклона трапеций (рис.

Для

вычисления

pSm

(t) и

pCos (t) по

формулам

(III.71) и

(III.72)

следует пользоваться

таблицами

s iп jc

cos х

функций ----- ,

----- ,

приведенными в [8—12, 38].

 

 

 

 

По найденной кривой р(/) опре­

 

 

 

деляется

передаточная

функция

 

 

 

В(/со). Если

В(/to) = Рв (ы) + jQB(со),

то

р в (со) = f р (/) cos сotdt;

 

о

(111.73)

оо

QB ((£>)= —J р (t) sin (otdt.

Рис. 111.42. Трапецеидальная ве­ щественная (мнимая) типовая частотная характеристика.

Очевидно, что при графоаналитическом методе вид функций

Лв(со), Qb (со) так же, как pCOs(0> Psin(0> можно найти, разложив р(/) на трапецеидальные составляющие по формулам [28]:

n / \

■V"'

о ^

sin/,o) sin A/fi>

 

Р в И

=

Роiti

J ------ т~*~ '•

 

Qb (®) — — Ро(7(

 

tid)

sin Д;со

1 ]

 

 

Д;(о

 

 

i=i

 

 

 

где через Рог обозначена

высота

трапеции;

значения А* и

опре­

деляются наклоном трапеции (см. рис. III.42).

При определении функций р(/), Рв (со) и Qb(co) для систем автоматического управления движением судов по курсу следует соблюдать высокую точность вычислений и графических построе­ ний. Это объясняется тем, что для систем автоматического регу­ лирования с большими инерционностями, которыми являются авторулевые, функция р(/) затухает через большие промежутки вре­ мени (часто больше 100 с). Поэтому высокочастотные составляю­ щие Рв (со) и QB( со) будут иметь большие амплитуды, и незначи­ тельная ошибка, допущенная при интегрировании функции р(/)

159

в области больших t, приведет к значительным погрешностям. При вычислении оптимальных передаточных функций для следящих сис­ тем, которые рассмотрены в [19, 28], функция ${t) практически затухает через несколько сотых долей секунды и приближенная an­

al

Рис. II 1.43.

Кривые

импульсной

переходной

функции 0 (0

и функции

 

 

Qb(to), вычисленные по £(<)•

 

проксимация

кривой

р(t) не

приводит

к заметным

погрешностям

в вычислении Рв( со) и Qb ( со) .

На рис. II 1.43,

а изображены две кривые импульсной переход­

ной функции р(0.

отличающиеся друг от друга только в области

больших значений t, а

на рис. III.43, б приведены вычисленные

по р(0

функции Q b ( co) .

Из сравнения кривых видно, что незначи­

тельное

отличие кривых

р(0 приводит к весьма существенному

изменению вида кривых Qb(со).

160

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

по выражениям

А

( с о )

(II 1.74)

 

Л(со) = - ^

- ;

 

^ 4 - ф

( с о )

 

 

ф ( с о ) = ф в ( с о ) — Ф ф ( с о ) .

(III.75)

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

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

5е(«) =

 

5ф (и) +

WfUm)

Sf (со), (II 1.76)

1+ г к. у(/со)

\ + WK.y (/<о)

 

 

 

где WK.y(j<a) — передаточная функция разомкнутой

системы с уче­

том введения корректирующих устройств. Например, при коррек­ ции системы последовательным устройством по первой производ­

ной она равна

(/со) • (1 —Ci/со).

WK.у (/со) = Фс. с (/со)

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

St (/со) =

1

W*. у ( / с о )

Sf (со). (III.77)

З ф ( с о ) +

1 + W K. у ( / с о )

1 +

WK. у ( / « в )

 

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

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

W'n.y (/®)= 1+Ci/co.

6 С. Я. Березин, Б. А. Тетюев

161

Значение ct необходимо найти из равенства

 

 

 

 

дс

= о,

 

 

где

 

 

 

дс±

 

 

 

 

У

Г

 

 

Wf исо)

 

 

С—

' +

Wк. у (/о))

5ф(й-

'S/ (со)| cfco

 

 

1 + W K. у (у со)

 

 

 

Пусть, например, выражение для квадрата среднеквадратичной ошибки системы автоматического управления, найденное аналити­ ческим путем, равно [19]

K{i l AK) {S9 ^ m + K ( l - A ) ] + S,(« 0 ) - ^ ] , (III.78)

где 5 ф(со)=2ро — спектральная плотность управляющего воздей­ ствия; S/(co) = jV— спектральная плотность помехи;

 

 

 

 

 

Т

=

Т 1+ Т 2\

 

 

Т\,

Г2 — постоянные времени системы;

 

 

 

 

 

 

 

А _

J j]j_ ;

 

 

к =

т ^1^3

•эквивалентный

коэффициент усиления, величину ко-

торого можно менять; klt k2,

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

отдель­

ных элементов.

 

(III.78) и обозначая через L отношение

 

Преобразуя выражение

S f ((D)

 

_

_ _ _ _

N

 

, получим

 

5 Ф (со) Г 2

которое равно

 

 

 

 

 

2|3стГ2

J

 

 

 

 

 

 

с2= рстГ

1 +

АГ(1 - - 4 ) +

L K 2

(III.79)

 

 

 

 

 

 

 

2/С (I —ЛАГ)

 

 

Дифференцируя

выражение

(III.79)

по коэффициенту

усиле­

ния К и приравнивая производную нулю, найдем выражение для коэффициента усиления, соответствующего минимуму ошибки сис­

темы:

1

 

К min

А + V a + L '

 

 

 

Подставив Кпйп в (III.79), нетрудно найти минимальное значе­

ние среднеквадратичной ошибки

 

c2min = ^ f ( A

+ 2]/^ A + L + l ) .

(III.80)

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

162

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

Синтез корректирующих устройств. После построения оптималь­ ной амплитудно-фазовой характеристики замкнутой системы необ­ ходимо перейти к выбору корректирующих устройств. Методика и порядок выбора необходимых корректирующих устройств будут такими же, как и в случае синтеза по заданному переходному про­ цессу при скачкообразном управляющем воздействии (см. § 11).

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

Рис. 111.44. К построению модуля передаточной функ­ ции ошибки системы автоматического управления дви­ жением судна по курсу.

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

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

Ф Л -дЛ =

(/со) [К о (/СО) + feo. с]

*

1 + Кр(/со)[У0(/со) + йо.с] ’

При

п(/<*>)- ^*0

^э(/“)

ф у©)

(III.81)

1 + К р (/со)

 

 

. ^0

 

* Величина k0.c в выражении для ф (/ю ) Сможет быть выбрана в

пределах

от 0,2 до 1, так как она не оказывает влияния на форму частотной

характе­

ристики.

 

 

6*

163

Для определения зависимости Ф(/а>) =/(£о) удобнее пользо­ ваться графоаналитическими вычислениями, производимыми на комплексной плоскости.

Построим

обратную

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

Ур(/о)) (рис.

III .44). Тогда

числитель в выражении (111.81) для

некоторой частоты со* получится в результате простого перемно­ жения векторов OA = Y°(j(Oi) и OB— — Jr k0 с:

 

К

6 Ъ = о А о в = (/ш()

с •

Рис. II 1.45. Кривые изменения модуля передаточной функции ошибки в зависимости от коэффициента уси­ ления.

Искомое значение модуля передаточной функции ошибки * для некоторой частоты со» и k0 определится как отношение векторов

оЪ

. Для других значений со и k0 вычисления следует повторить.

DE

рис. II 1.45 приведены кривые изменения модуля передаточ­

На

ной функции ошибок для судна с обратной передаточной функцией

/© ( 1 + 3 0 ,5 /(0 ) _

'0 ,1 (1 4 -6 ,5 /© )

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

* Метод определения ф е (/ю )= /(М при возмущающем воздействии на суд­ не изложен в § 11.

164

со = ©о = 0,07 рад/с). Из этого же рисунка также следует, что изме­ нение к0 благоприятно действует на ошибку только в некоторых пределах и излишнее его увеличение нецелесообразно, так как не приведет к желаемому уменьшению ошибки.

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

0

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

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

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

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

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

165

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

Для настройки системы в процессе наладки рекомендуется предусмотреть возможность регулировки параметров выбранных корректирующих устройств в пределах ±20—30 %• После выбора корректирующих устройств необходимо произвести контрольную проверку устойчивости системы, а на ее макете или на аналоговой электронной вычислительной машине — проверку качества проек­ тируемой системы для выбранных параметров корректирующих устройств. При этом нужно также проверить качество работы сле­ дящей системы управления рулем с энергетической точки зрения, т. е. с точки зрения удовлетворения выбранной мощности, и число перекладок руля.

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

На практике могут встретиться два случая: первый, когда за­ даны допустимая среднеквадратичная ошибка работы системы и спектральная плотность случайного воздействия, и второй, когда ни спектральная плотность, ни корреляционная функция не заданы, но известны допустимое значение среднеквадратичной или устано­ вившейся ошибки и максимальные значения скорости и ускорения внешних воздействий.

Начальным этапом синтеза является определение значений амплитуды и частоты эквивалентного гармонического воздействия.

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

Если спектральная плотность управляющего воздействия не имеет явно выраженного резонансного пика или он лежит в об­ ласти частот, не пропускаемых объектом регулирования (кривые I и IV на рис. III.46), то величина резонансной частоты амплитудночастотной характеристики объекта регулирования может быть при­ нята в качестве расчетной, так как при этой частоте ошибки ра­ боты системы будут близки к максимальным.

Если амплитудно-частотная характеристика объекта регулиро­ вания не имеет резонанса, а полоса частот, пропускаемая объектом.

166

достаточно велика (кривая II на рис. III.46), то в качестве рас­ четной частоты целесообразно принять частоту резонанса спект­ ральной плотности.

При наличии резонанса и достаточно широкой полосе пропус­ кания частот объектом регулирования (кривая III на рис. II 1.46) расчет необходимо проводить для нескольких частот между резо­ нансами кривых 5 ф(со) и Ао. per(ft>), или для наибольшей частоты с последующей проверкой величины ошибки для других частот.

Если кривые / — IV резонансных пиков не имеют, то в качестве расчетной частоты можно принять максимальную частоту, пропус­ каемую объектом регулирования.

Рис.

II 1.46.

Амплитудно-частотные

характеристики

объ­

екта

регулирования

Л 0 . р е г ( и )

(кривые IIII) и

спект­

ральная

плотность

управляющего

воздействия

S<p(co)

 

 

 

(кривая

IV).

 

 

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

/

1 °°

 

-S-J S„(<o)d<D.

(HI.82)

 

Если задана корреляционная функция внешнего воздействия, то

Ф = |Л ? / (0).

(II 1.83)

В рассматриваемом случае синтез корректирующих устройств рассчитываемой системы необходимо осуществлять в следующем

порядке:

1) установить величину эквивалентной частоты и амплитуды внешнего воздействия;

167

2) пользуясь частотными методами, выбрать корректирующие устройства, принимая за основу заданную величину среднеквадра­ тичной ошибки;

3) построить кривую | Фе (/со) |2 с учетом выбранных корректи­ рующих устройств. Перемножить ее ординаты с ординатами (для соответствующих частот) кривой спектральной плотности S,p(to). Вычислить площадь получившейся кривой (ограниченную осью абсцисс). Найти значение среднеквадратичной ошибки

4) сравнить полученную ошибку с заданной. Если значени среднеквадратичной ошибки получилось меньше заданного или равным ему, то на этом расчет может быть закончен. Если же

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

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

следует выбирать по эквивалентному гармоническому воздействию

 

^2

( а расчет нужно

с частотой созкв= Втах и амплитудой Лэкв= _

^max

16тах

 

вести по установившейся ошибке, составляющей 2,5—3 средне­ квадратичные ошибки.

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

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

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

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

168

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