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

книги из ГПНТБ / Капустин, К. Я. Плавучие буровые установки и буровые суда

.pdf
Скачиваний:
36
Добавлен:
21.10.2023
Размер:
9.98 Mб
Скачать

Ёсли при включенной АСС принять

скорость

выборки каната

якорными

лебедками 1 м/с, т. е. в уравнении б= kt

принять

коэффициент k=l,

то б = /.

Следовательно, можно считать изменения восстанавливающих сил

 

Р2Т) +

сц3

 

 

одинаковыми с изменениями сил натяжения якорной системы, вызванными авто­ матическим устройством

|}2х -)- сх3 = yt + dt3, т. е.

у — р2 и с = d.

Уравнение движения БС после включения АСС будет

400т)" + 6т) + 0,216т]3 ----- 100 — (6/ + 0,216/3) .

Частное решение этого уравнения по способу прямой линеаризации: 0,0236А5 + 1,85А3 + (0,925/ + 0,0334/з_ 15,4) А2 4- 36А +

4- 1,3/3 4- 34,7/ — 6С0 =

0.

 

Решение этого уравнения в виде зависимости

Л(1) приведено на рис.

93.

Как видно из рисунка, при /тах=6,56 с амплитуда

отклонения равна нулю,

при

/тах>6,56с уравнение (200) применять далее нельзя. Значение /шах соответст­ вует половине колебания, так как полное колебание 7= 2 /т а х=13с.

После включения автоматической системы и возвращения БС в зону нечув­

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

необходимое до ее затухания:

0

In

1,45.

,1 -13

7,2

 

В рассматриваемом случае

система

совершит два колебания, а время зату­

хания пТ=20 с.

 

 

 

Составление и решение уравнения

автоматической системы стабилизации судна (АСС)

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

ветствии со схемой рис. 89.

сле­

Уравнение каждого звена в операторной форме будет

дующим.

можно

Первое звено — г и р о в е р т и к а л ь . Гировертикаль

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

0вх =

1 0ВЫК-

(206)

Если же за входную величину считать горизонтальное смеще­

ние БС Хвх, в этом случае можно записать:

 

■^ВХ

^Овых ' 1 >

(2 0 7 )

где Н — практически глубина воды в месте бурения.

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

 

 

U№ — k2(0сс.д — Эссп),

(208)

где

kT— коэффициент

передачи сельсинного устройства;

0ссп —

угол

поворота ротора

сельсина-приемника; 6 c c d — угол поворота

ротора сельсина-датчика.

Так как судно должно удерживаться при бурении лишь в одном положении, то при отсчете поворота ротора сельсина-приемника (который неподвижно скрепляется с приемной осью) это поло­ жение можно принять за нулевое. В связи с этим, угол 0ссп — так как для увеличения чувствительности АСС между гироверти­

калью

и ССД может быть установлен редуктор (инерцией кото­

рого пренебрегаем), то:

 

 

 

t/сп = &1&20ССД,

(209)

где

— коэффициент механической

передачи редуктора

между

 

гировертикалью и ССД.

у с и л и т е л ь . В

первом

Третье звено — э л е к т р о н н ы й

приближении усилитель можно считать безынерционным звеном, поэтому

 

 

 

£/,.y = M£/cn-e«e • *о.с),

 

 

(21-°)

где

Нэ.у — напряжение на выходе

 

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

усиления

усилителя;

0леб — угол

поворота

барабана

лебедки;

&о.с — коэффициент

передачи ЛВТ

линии обратной

связи.

При этом существует зависимость:

 

 

 

 

 

 

 

 

U%=

0 леб

* о . с

 

 

 

( 2 П )

Четвертое звено — п о л я р и з о в а н н о е

реле. Реле

является

нелинейным звеном

в АСС

со

ступенчатой

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

(рис. 94)

и может быть записано лишь в общем виде:

 

 

 

 

 

ДР =

А(ДЭ.У).

 

 

 

(212)

Пятое

звено — и с п о л н и т е л ь н ы й э л е к т р о д в и г а т е л ь .

Если

пренебречь электромагнитной

постоянной

времени цепи

якоря, то приближенно уравнение двигателя имеет

вид:

 

 

 

 

Р(ТмР+ 1)0ДВ=

 

 

 

(213)

где

Тм— электромеханическая

постоянная

времени

двигателя;

&4— коэффициент передачи; 0ДВ— угол поворота

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

Шестое звено — л е б е д к а .

Пренебрегая

инерционностью ре­

дуктора и барабана лебедки:

 

 

 

 

 

 

 

 

 

 

0деб = ^50дв-

 

 

 

(214)

14 Капустин К. Я-

209

Учитывая диаметр барабана лебедки для вытравленного или выбранного троса и исходя из приведенных ранее соображений,

получим: хтрг= D0-ie6,

обозначая: kb = k$D,

(215)

получится хтр = /гБ0дп.

 

Седьмое звено — я к о р н а я

с и с т е м а . Линеаризованное

уравнение движения судна при действии АСС на якорную систему запишется в виде:

Мр2+ 2NP + Кхвых = Afcip.

(216)

Рассматриваемая автоматическая система стабилизации судна

включает в себя несколько

нелинейных звеньев (четвертое и ше­

 

стое) .

При

этом реле

(четвертое зве­

 

но) — существенно нелинейно.

 

В таких системах, как известно, мо­

 

гут появляться

автоколебания, ампли­

 

туда и частота которых зависят от па­

 

раметров системы. В связи с этим при

 

рассмотрении приведенной выше АСС

 

могут

возникнуть следующие задачи:

 

определение

границы

возникновения

 

автоколебаний, определение зависимо­

 

сти амплитуды и частоты автомехани­

 

ческих колебаний от параметров АСС,

 

оценка

качества

переходного процесса

 

в зоне отсутствия колебаний.

 

Следует отметить, что возникнове­

 

ние в

системе

АСС

автоколебаний

Рис. 94. Характеристика реле

(тем более устойчивых симметричных)

как звена АСС.

нежелательное

явление, так как это

 

приводит к длительной

знакоперемен­

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

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

210

лагаются близкими к гармоническим колебаниям, состоит в том,

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

которого ограничиваются членами первой кратности. В рассматри­

ваемой системе АСС выходная величина 0 Р реле

при любом за­

коне изменения входной величины будет иметь

прямоугольную

форму. В случае периодических

колебаний с

основной частотой

со функцию б/р можно разложить в ряд Фурье:

 

Up =

00

 

 

Апsiri (nat + ct„),

 

 

П~\

 

 

где А п — амплитуда гармоники

п\ ап — фаза

гармоники.

Так как F (Uay) нечетна и однозначна, то гармоники будут нечетные п= 1, 3, 5, 7...

На выходе линейного звена пять получится:

оо

0 дв = 5j Ь"А ‘п (m t + “ я + Р»)'■

(217 )

л = 1

 

Из уравнения (217) можно получить характеристики для зве­ на пять:

а) амплитудная характеристика

Ь„ =

(218)

Р ( Т „ р + 1) р - j n O )

/Ш |

7'-/220)“ -|- 1

б) частотная характеристика

 

 

ф = Arg

h

arctg Гм(о.

(219)

 

Р(ГыР-И)

Как видно из формулы, для Ьп амплитуды высших гармоник (3-й, 5-й и т. д.) на выходе линейного звена пять проходят с го­ раздо меньшим усилением, чем первая гармоника (га=1). Высшие гармоники практически не пропускаются. При прямоугольной входной величине сама амплитуда Ап для высших гармоник тем меньше, чем больше номер гармоники. В связи с этим можно огра­ ничиться лишь первой гармоникой разложения, величину Одр считать близкой к синусоиде:

0дв = M i sin (ait -faj-t- Pi).

(220)

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

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

14* 211

колебания, весьма близкие к простым гармоническим колебаниям с той же частотой, т. е. квазигармонические.

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

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

(221)

пА

 

Тогда уравнение четвертого звена запишется в виде

 

Up = qUs.n

(222)

После линеаризации всех звеньев системы можно составить характеристическое уравнение системы, поскольку решение мож­ но найти в форме x=Asmu>t с постоянной амплитудой А. Из этого следует, что коэффициенты q, К и Xi (216), (221) так же величины постоянные. В связи с этим уравнения (207) — (218) представляют собой системы линейных дифференциальных урав­ нений с постоянными коэффициентами. В отличие от обычных систем в этой системе постоянные коэффициенты q, к и 7,1 неиз­ вестны, и они могут быть определены после определения ампли­ туды А.

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

(MlP* + 2Np + Г) [р (Тмр + 1) + k3kMo,q] + -lfeafe3y ^

g. -

0. (223)

Если обозначить kl kl ka ki къ —Ki и &3й4&5&о.с —K2,

то

оконча-

н

 

 

теЛьно получится

 

 

212

MJnP* + (2NTM+ Mx) , f + (MxK2q + 2N + TJ?) p2 +

+ (2NK2q + Я2) p + K2q№ + KjZq = 0.

(224)

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

 

 

х (q, Я, Яъ со) -f jy (q, Я, Яь со) = 0;

 

 

(225)

 

х = МхТисо4- (MxK2q + 2N + TJ*) о)2+

K2qW + K&q = 0;

(226)

 

у = — (2NTM+ Мх) со3 + (2NK2q + Я2)со =

0.

 

(227)

 

После преобразований и подстановки получится

 

одно

урав­

нение:

 

 

 

 

 

 

-

8N3K2TMq -

2NKlM\q* - 4N*K2Miq -

4NMxTwK2ql2+ M\K2q}? +

+

W T l K t f q + W T ^M J C ^q +

-

4 А 2Г МЯ2 -

2NT4M%*-

 

 

M\K2q>?2NMX№ = 0.

 

 

(228)

Приведенные

выражения (226) — (228)

можно

использовать

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

 

Ь ,= - ----— =

=

 

(229)

 

Зсо У

2со2+ 1

 

 

Если потребовать, чтобы коэффициент усиления первой гар­

моники был bi^>b3, то

 

 

 

 

 

 

------

К

» -------

k4

(230)

со V

Г2со2+1

Зсо У

 

+1

 

После подстановки в (230) значения со получится:

1

80 (2NK2q + Я2) Т \ +

8 (2NTM+ MJ » 0.

(231)

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

Определение областей устойчивого

равновесного состояния системы

Чтобы найти зависимость между амплитудой и частотой авто­ матических колебаний от различных параметров системы АСС, нужно рассмотреть изменение трех основных характеристик: Ki

213

<

коэффициента усиления ACC; К2 коэффициента обратной связи; с — постоянную реле. Остальные параметры можно считать не­ измененными. Так как исследование в общем виде (без числен­ ных значений параметров) затруднительно, в дальнейшем в ка­ честве примера принимается якорная система бурового судна, рассмотренная ранее (на стр. 225), а значения других параметров системы АСС принимаются: Гм=1 с; 2А = 40 тм/с2.

Рис. 95. Зависимость

амплитуды

и

Рис. 96. Зависимость

амплитуды и

частоты

автоколебаний

от коэффи­

частоты

автоколебаний

от коэффи­

 

циента Ki.

 

 

циента К^.

 

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

2 — частота

ко­

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

2 — частота

 

лебаний.

 

 

 

колебаний.

 

При

подстановке указанных

численных

значений

параметров

в уравнение (228) получается следующее уравнение с неизвест­

ными коэффициентами Ki и К2'

 

 

 

 

 

Ав +

2774А4 + 4900Л:2А3 — 297 • 108/СИ8+

106 •

ЮМ2 +

+

892 • 103/С2Л - 1 1 ,6 - 1 0

6/СИ + 6,4 •

10®/С| = 0.

(232)

 

 

 

40/г2. ——+ 6 + 0.154/42

 

 

 

 

 

 

Сй2 = --------А------------------ .

 

 

 

 

 

 

 

440

j

 

 

 

Если задаться

различными численными

значениями

Ki и К2,

то можно

решить

эти

уравнения

относительно

А

и со,

причем,

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

внимание

лишь действительные

положи­

214

тельные корни уравнений (так как только они соответствуют физическому смыслу задачи). На рис. 95, 96 показана зависимость амплитуды и частоты автоколебаний от коэффициентов Ki и Кг, найденные из уравнений (232).

Аналогичные

уравнения

получаются для параметра

с, если

принять Лл=1

и /(2= 0,046.

 

 

 

 

 

А6+ 2780А4 — 351 •

ЮМ3 +

107 • 103А2—'96 •

104Ас + 22 •

103с2 = 0.

 

 

 

 

 

 

 

 

 

 

 

(233)

 

 

 

со2 =

0,0053 — +

0,0136 + 35 • 10-5Л2.

 

 

 

 

 

 

 

А

 

 

 

 

На рис. 97 показана зависимость

A=f(c)

и со= / (с).

 

На рис. 95—97 выделена

 

 

 

 

область отсутствия

автоко­

 

 

 

 

лебаний. Как отмечалось,

 

 

 

 

эти колебания весьма не­

 

 

 

 

желательны для АСС буро­

 

 

 

 

вого судна, и значения па­

 

 

 

 

раметров

системы

должны

 

 

 

 

находиться в зоне, не имею­

 

 

 

 

щей автоколебаний. Поэто­

 

 

 

 

му в

дальнейшем

система

 

 

 

 

рассматривается

 

лишь

в

 

 

 

 

этой зоне.

 

 

 

 

 

 

 

 

 

 

При этом следует отме­

 

 

 

 

тить следующее:

автоматиче­

 

 

 

 

Из

теории

 

 

 

 

ского

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

нели­

 

 

 

 

нейных систем известно, что

 

 

 

 

для линеаризованной

систе­

 

 

 

 

мы АСС

можно

применить

 

 

 

 

критерии устойчивости, ана­

 

 

 

 

логичные

с

применяемыми

 

 

 

 

для линейных систем, на­

 

 

 

 

пример, критерий

Гурвица.

 

 

 

 

Вне

зоны

автоколебаний

 

 

 

 

при любых начальных уело- *

Рис. 97.

 

 

 

ВИЯХ колебательный процесс

Зависимость амплитуды и час-

будет

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

зату-

 

тоты от параметра с.

 

хающим или расходящимся.

В первом случае характеристическое уравнение должно удов­ летворять этому критерию при всех значениях амплитуды; во вто­

ром случае оно не будет удовлетворяться.

Уравнение (224) про­

веряется по критерию Гурвица.

 

 

 

Запишем четыре условия устойчивости:

 

1)

М}Тм > 0~2NTM+

Л42 > 0 ;

 

 

2)

4№ТЫ-|- 2ЛТ^Я2+

М\Кгд + 2NM1

>

0;

215

3) (2NK2q : ^) Ла - (2NTM-f Mf) (K2qV + Ktfq) = A3 > 0; (234)

4) (ktfW + Kitiq) A3 > 0.

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

На рис. 98, 99 приведено численное решение уравнения (234) третьего условия в зависимости от различных значений коэффи-

Рис.

98. Численное решение уравнения

Рис. 99. Численное решение

уравне-

(234)

в зависимости от Ki (при /С2==

ния (243) в зависимости от

л 2 (при

 

=0,046).

(4 = 0,05).

 

циентов К\ и Кг- В зоне, прилегающей к оси абсцисс, Дз<0, что соответствует неустойчивому состоянию системы; в зоне, приле­ гающей к оси ординат, Дз>0, что соответствует устойчивости си­ стемы. Необходимо помнить, что устойчивая область системы обладает этим свойством при любом положительном значении амплитуды колебания — это соответствует затухающему переход­ ному процессу.

216

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

Оценка качества переходного процесса

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

х — а0е^л' sin (сол/ -f- ф0).

(235)

Считается, что в нелинейной системе функцию x(t) ввиду плав­ ного ее характера можно приближенно также описать синусоидой

спеременной частотой со и переменным показателем затухания

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

комплексных корней в левую часть характеристического уравне­ ния (224) подставляется значение Я = £ + /со. После подстановки получается выражение, аналогичное (225), состоящее из двух ча­ стей: вещественной и мнимой

х (q, X, Xv со, I) + }у (q, X,

со, |) =

0.

(236)

В результате получается два уравнения:

 

 

М { Г ^ - 6М1Гмсо^ 2 + М ,7>4 4- (2NTм +

• (g3-

3£о>г) +

+ (M1K2q + 2N !

T J 2) (I2- со2) -|- |

(2NK2q +

Я2) +

K2qX2+

 

+ KxX\q = 0;

 

 

(237)

M J M(4£3co-

4£co3) + (2NTM+ ЛД) . (3£2<o - со3) +

+ + 2N + TJ?) •2gco+ (2NK2q + Я2) со= 0.

217

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