
книги из ГПНТБ / Капустин, К. Я. Плавучие буровые установки и буровые суда
.pdfЁсли при включенной АСС принять |
скорость |
выборки каната |
якорными |
лебедками 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. Ве личина этого коэффициента равна
4с |
(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) |
|
|
Зсо У |
9Г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