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

книги из ГПНТБ / Пивоваров В.А. Проектирование и расчет систем регулирования гидротурбин

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

где индексом п обозначены номинальные (расчетные) значения соот­

ветствующих величин.

Разделим выражение (1.2) на Міп и, учи-

da>

=

dx

получим

 

 

тывая, что -jß-

 

 

 

 

 

 

 

/ш„

dx

 

 

 

 

Миг

~di — mt tng,

(4.2)

Величина ■

■= Ta,

измеряется

в с и называется постоянной

Min

 

Физически

Та— это время, в

течение

времени гидроагрегата.

которого агрегат достигает номинальной скорости вращения при

номинальном моменте.

В практических

расчетах величину Т а

удобнее определять по

формуле

 

 

 

Т1 а=

о п Ч

(4.3)

 

365Рп

где GD2 — маховой момент ротора агрегата в тс-м2; Рп — номи­ нальная мощность агрегата в кВт; пп — в об/мин.

В результате в относительных величинах уравнение машины

имеет вид

 

Taxp — mt nig.

(4.4)

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

 

 

 

Я

_Q

h =

н - и п

 

(4.5)

 

 

 

Qn

Н п

 

 

 

 

 

 

 

 

 

Из

(4.1)

и

(4.5)

имеем:

 

со =

сол (1

+ х ); Q =

qQn, Н =

= Нп (1 + h);

гр = %ціп. Подставив

эти значения в

выражение

(1.3), получим

 

 

 

 

 

 

 

 

 

 

 

М / = ( У 0 ^ ) £ ( 1 + Д І Х_

(4.6)

В (4.6) член в скобках представляет собой номинальное зна­

чение момента турбины Mtn. Значит,

 

 

 

 

 

 

 

mt =

 

Я

X.

'

(4.7)

откуда

следует,

что

момент

іщ является

функцией

нескольких

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

 

Щ = тю+

 

А/п,;

)

(4.8)'

 

 

I л

 

J

 

тв = тео• +

A /rtg ,

 

где mt0

и mg0 — значения моментов

в

рассматриваемой точке

(в точке

линеаризации), причем

 

mt0

— tng0.

 

185

Изменение момента нагрузки можно представить в виде

Amg = ^ A x + Az(t),

(4.9)

где Az (t) — изменение момента нагрузки из-за включения и от­ ключения потребителей электроэнергии.

Величина dmjdx, согласно выражению (1.13), представляет собой коэффициент саморегулирования es нагрузки. Подставив (4.8) в (4.4), получим

Т аАхр = Amt — Amg.

(4.10)

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

функцией

четырех

независимых переменных

Апг( = f (Ay, Дер,

Ah, Ах).

Здесь Ау

и Дер — относительные

изменения соответ­

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

Ч' -- о

о )

^max

Фтіп

оо

где Фшіп, фтах — минимальный и максимальный углы разворота лопастей.

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

dm, .

dm,

dm,

dm, .

 

Ami==~W АУ + ^

Аѵ + ~жАІІ + ~Ш~Ах-

(4Л1)

Поскольку изменение расхода Аq зависит от тех же величин, что и Amt, т. е. Aq = f у, Дф, Д/г, Ах), то по аналогии с (4.11) можно записать уравнение расхода, без которого невозможно определить влияние гидравлического удара на динамику регули­

рования

1

 

-

і* = | - Л Н - £ Д ф + І Д А + Ц- Д* .

(4.12)

В выражениях для Аш( и А q частные производные представ­ ляют собой постоянные коэффициенты, которые определяются значениями этих производных в точке линеаризации.

Из гл. 1 (см. рис. 2) и из рассмотрения принципиальных схем регуляторов скорости поворотнолопастных турбин (см. рис. 70, 72) известно, что величины Аф и Ау связаны между собой комбина­ торной зависимостью, а сервомотор рабочего колеса, с помощью

186

которого осуществляется разворот лопастей, является следящей системой (гидроусилителем) со статическим коэффициентом усиле­

ния kK= и постоянной времени Тк. Следовательно, влияние

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

[Ткр + 1) Лер = kKAy.

(4.13)

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

Рис. 77. Структурные схемы управления сервомотором рабочего колеса поворотнолопастной гидротурбины

схем (см. рис. 72), у которых управление кулачком комбинатора осуществляется сервомотором направляющего аппарата. Если же кулачок управляется от промежуточного сервомотора, что соот­ ветствует схеме регулятора, показанной на рис. 71, то в этом слу­ чае входной величиной для гидроусилителя рабочего колеса является смещение поршня Ау \ промежуточного сервомотора. Последний одновременно управляет также и гидроусилителем на­ правляющего аппарата с коэффициентом усиления, равным еди­ нице, и постоянной времени Ту.

Структурная схема такой связи приведена на рис. 77, а. Ис­ пользуя правила преобразования, ее можно представить в виде, показанном на рис. 77, б. Отсюда имеем:

(TKp + l)A<p = kK(Tsp + \)Ay. '

(4.14)

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

Ah = WyAq.

(4.15)

187

Таким образом, динамические характеристики гидроагрегата как объекта регулирования определяются следующей системой уравнений:

ТаАхр = Amt Artig-,

\

,

dm,

dm,

 

dm,

dm,

Am‘ = ~W Ay + ~W Лф+

~дІГ AA+

~W Ax’

 

Aq = ж

А^ + ж Аф+ ж

 

(4.16)

 

А/г+ і І Ах;

 

 

A h =

WyAq]

 

 

 

(TKP + l)Aq> = kKAy

 

или

 

 

 

 

 

 

(TyP -f- 1) Аф =

к к {Т ир + 1) At/.

 

Из этой системы, исключая соответствующие переменные, при условии отсутствия возмущений [Az (t) = 0 ] можно получить несколько передаточных функций, характеризующих динамичес­ кие свойства агрегата в различных режимах его работы:

1) передаточную функцию гидротурбины при постоянной ско­

рости

Wty = AmtIAy

при

X = const;

 

 

 

2)

передаточную функцию гидроагрегата при постоянном

открытии направляющего

аппарата Wax =

AxlAmg

при

у —

const;

функцию гидроагрегата

Wa =

Ах!Ay.

 

3)

передаточную

Wty

На

практике наиболее важное значение имеют функции

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

Для большей наглядности, используя уравнения (4.16), соста­ вим структурную схему гидроагрегата как звена системы регу­ лирования. Эта схема для случая управления комбинатором от сервомотора направляющего аппарата показана на рис. 78. При этом учтено, что dmjdx = et и dmglдх = eg.

Как видно из рисунка, гидроагрегат представляет собой довольно сложное звено системы регулирования и его целесооб­ разно разбить на более простые звенья. Это можно сделать, если представить изменение момента турбины Amt в виде суммы Amt = = Атіу 4- Апгіх, где Аmty — изменение момента в условиях постоянной скорости вращения, которое зависит только от поло-

188

жения регулирующих органов, а Аmtx— приращение момента, вызываемое отклонением скорости при у = const (Ау = 0). Ис­ пользуя правила преобразования структурных схем, изложенные в гл. II, разделим на рис. 78 каналы воздействий по открытию и по скорости. В результате получим преобразованную структурную

Рис. 78.- Структурная схема гидроагрегата как звена системы регулирования

схему гидроагрегата,

показанную на

рис. 79, а. На этой схеме

с целью сокращения записи введены обозначения:

 

dmt — т,

dmt

= т,

dmt

Щк,

 

ду

■tyi

dq>

•йрі

dh

dg

_

dq___

_

dq___

_

dq___

dg

~ ЧУ'

0ф —

9<p’’

dh ~~ 9л’

дх ~

Из схемы легко можно получить передаточные функции про­ стейших звеньев. К ним относятся: передаточная функция Wg = = Дх/А/Пяд. генератора с энергосистемой

гос =

___і___ •

(4.17)

ß

Tap + eg ’

 

передаточная функция Wix = Am^/Ax турбины при постоянном открытии регулирующих органов

WiX = et + ЧхЩк \ - qJhW- >

(4-18)

а также передаточная функция Wiy. Структурная связь“ между этими функциями показана на рис. 79, б. Из этой схемы следует, что

W

= — У‘—

(4.19)

а*

l — WtxWg’

 

1 189

а передаточная функция \Ѵа = Ах/Ау

гидроагрегата

= WtyWax.

(4.20)

Определим значение функции Wty. Из рис. 79, а имеем

(,»+‘-,»>(ѵйл5'«р+1)"'Л

W;іу

(TKp + \ ) ( \ - q h\Vy)

 

+

 

(mty+ kЛ ?) - -> ■

Г р + Л

 

+

---------------

Акт лр

К

/

(4.21)

ГкР+1

 

 

 

 

 

 

 

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

Рис. 79. Преобразование структурной схемы гидроагрегата

постоянной скорости вращения. Однако для практических расче­ тов оно сравнительно сложно. Чтобы упростить выражение для

Wty,

положим

mty

^

Яи

что весьма близко

тіу-\~ К.ті<ѵ

 

Яу “Ь

 

 

 

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

 

W

іу

— Р

ТкР + 1

П +

(Xmm — qi,)Wy

(4.22)

 

 

 

у

Ткр + 1

[

1

- q hW y

 

 

где еу =

miy -f

kKmi{f — коэффициент

усиления

турбины по

моменту

при постоянной скорости,

 

 

 

 

 

 

%=

^ + V Ф

 

г ' — ПЦу Т

к •

 

 

 

 

 

 

1

к --- ~~р і

 

 

 

 

 

 

 

 

 

еу

 

 

190

Если подставить выражения (4.17) и (4.18) в (4.19), то с учетом,

что es et — еп,

будем иметь

 

W

=

___________ 1 -qnWy__________

(4.23)

wад: —

(ТаР + еп) — (qiiTap + qsmlh -j- enqh) Wy

 

Зная Wty и Wax и подставляя их значения в (4.20), легко можно получить передаточную функцию гидроагрегата с поворот­ нолопастной турбиной в общем виде, которая учитывает все ос­ новные параметры, оказывающие влияние на его динамические характеристики,

 

 

 

Ц7а =

еу

т'кр+ 1

 

 

 

l + (Knuh - q h)Wy

(4.24)

 

 

 

Т к Р + 1

ТаР +

вп— (дпТар + qxmih + enqh) Wu

Прежде чем перейти к рассмотрению передаточной функции

гидроудара

WtJ, выразим mth и qh через

величины, легко

опреде­

ляемые по данным универсальной ха­

 

 

 

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

 

модели

рабочего

ко­

 

 

 

леса,

 

которая

обычно имеется в рас­

 

 

 

поряжении

конструктора. Вид такой

 

 

 

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

для

модели

ради­

 

 

 

ально-осевой

турбины

уже

 

рассма­

 

 

 

тривался в гл.

1 (см. рис. 3). Для мо­

 

 

 

дели

 

рабочего

 

колеса

поворотноло­

 

 

 

пастной

турбйны

подобные

харак­

 

 

 

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

 

 

 

углов

разворота

лопастей

(обычно

 

 

 

углы задаются через 5°). Они назы­

 

 

 

ваются пропеллерными, т. е. в этом

 

 

 

случае имеется серия характеристик,

 

 

 

связывающих между собой Q(, п\, ам,

 

 

 

г|м при Я

=

const и ф =

const.

 

 

 

 

 

Чтобы определить mth и qh, а

 

 

 

также

другие

коэффициенты

(А,, еу,

 

 

 

et), входящие в выражение (4.24),-

 

 

 

необходимо по пропеллерным

харак­

 

 

 

теристикам

построить

моментные

и

 

 

 

расходные

характеристики модели и

 

 

 

турбины.

Их

 

вид

для двух I углов

 

 

 

(ср =

0 и ф =

10°) показан на рис. 80

 

 

 

и на рис. 81.

Они

получены следую­

 

 

 

щим образом. На обоих графиках по

 

 

 

горизонтальной оси отложены отно­

0. Л,

 

сительные

приведенные

скорости

 

/

 

' f f

 

тт

 

 

 

 

 

 

 

Рис. 80. Моментные хнрнктери-

Xi

П\ІП\п-

По

ОСИ ординат на МО-

стики модели поворотнолопаст-

ментной

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

отложены

 

ной гидротурбины

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

величины

приведен­

 

расходной — относитель­

ных

 

моментов

 

іп\ =

М\ІМ\п,

а

на

 

ные

 

приведенные

расходы

 

q{ — Ql/Qi,,-

 

191

I

Каждая кривая на графиках соответствует определенному от­ носительному открытию а = a ja wl = а/ап, взятому с пропеллер­ ной характеристики. За номинальное значение аып (ап) принято открытие направляющего аппарата, соответствующее максималь­ ному ходу поршня сервомотора. Такой выбор ап обусловлен тем,

что основные параметры регу­ лятора скорости (Ьр, bt) отне­ сены к полному ходу сервомо­ тора. Величины пі,и М [,„ Q[n соответствуют номинальным- (расчетным) приведенным зна­ чениям скорости, момента и рас­ хода, причем значения М\ и М[п определяются по формулам:

 

 

 

 

М[

ЗОу

 

Qi

■Пи;

(4.25)

 

 

 

 

_

ЗОу Qi»

 

 

 

 

 

ЛмпГ

 

 

 

 

М1п =

 

Пп

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

пг[

 

 

'Хм-

(4.26)

 

 

 

Построение

_

зависимостей

 

 

т[ = f

(х[),

ql =

f

(x'l)

ведется

 

 

таким же образом, как

это де­

 

 

лалось

в гл.

1 при

построении

 

 

статических характеристик тур­

 

 

бины.

 

 

 

 

 

 

 

 

 

Однако в выражение (4.24)

 

 

входят

не Приведенные

величи­

 

 

ны

модели,

а

значения пара­

 

 

метров

натурной

турбины, по­

Рис. 81. Расходные характеристики

этому

необходимо выразить т\,

х{,

q{

через

соответствующие

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

им

значения

mh

х,

q. Это

бины

 

 

 

можно сделать,

если использо­

вать формулы пересчета (1.4), (1.5)

и (1.7). Действительно,

разделив Q на Q», п на

пп и М ( на М (п, будем иметь:

 

 

'

л Г Ж

п

'

 

н<і\

Ч =

/

ігп > ^

= х'Ѵ т г п->

=

Но учитывая,

что Н!Нп — 1 + /г,

п!пп = 1 + х,

= т{/%м и полагая It =

Х„,

получим:

 

 

Ші =

nit

: _

\+ x .

?i =

 

1+/J ’

X i

V l + h ’

у П + л '

(4.27)

а qi/xi =

(4.28)

192

Таким образом, зависимости т{ = / (xi), qi — / (х{)> представ­ ленные на рис. 80 и 81, одновременно являются статическими ха­ рактеристиками натурной турбины:

т

1+ х ' и ------..= f ( 1 'М при а — const,

1+ h

Ѵ Т + Т і,

по которым и определяем /га,Л и qh. Пусть, например, турбина работает в режиме, соответствующем а = 0,647 и ср = 10°. Для номинального напора в установившемся положении х = 0, h = 0 mt = tnt0 и q = q0. В точках с указанными координатами прове­ дем касательные к кривой а = 0,647, тем самым заменяя нелиней­ ные зависимости линейными. Уравнения этих касательных будут:

пц

__

дпц

1+ х

1 - 1- л

17110 —

дх

/Г + 1

я

 

дс} f

(4.29)

 

1 - j - X

КГ+л

9о — дх (j/Т + л

откуда

р

 

 

-gjj- = mth = tnt0— 0,5e<;

 

(4.30)

Jjf ==qh = 0,5(qo — qx).

Теперь рассмотрим передаточную функцию гидроудара Wy. Вначале решим задачу без учета упругости стенок трубопровода и воды, т. е. будем считать удар жестким. Воспользуемся основ­ ным законом механики F = mdv/dt, которому подчиняется дви­ жение воды в водоводах. Здесь в данном случае пг — масса, а d vld t•— ускорение воды. Предположим, что по всей длине L трубопровода площадь поперечного сечения одинакова и равна А. Тогда:

_AL_. dv_ __

1dQ '

F =■— AHA.

m ~

g

dt ~

А dt '

 

 

Учитывая, что

Q =

qQn и

АН = 1іНп, получим

 

 

h — — ®nL

dq

(4.31)

 

 

 

gHnA

dt

 

В этом выражении коэффициент перед

имеет размерность

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

Т

Qn

U

(4.32)

* w

gHn

2 Ai

 

 

 

13, В. А. Пивовароі

1 9 3

где L(- и A t — длина и площадь поперечного сечения отдельных участков трубопровода.

Физически Тш— время, в течение которого

расход

воды

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

напора

изме­

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

Итак, в форме жесткого удара передаточная функция имеет

вид

 

Wu = - T aP.

(4.33)

Если же учитывать упругие деформации стенок трубопровода и воды, то любое изменение давления или расхода распростра­ няется по длине трубопровода не мгновенно, а с некоторой конеч­ ной скоростью а, величина которой зависит от диаметра трубо­ провода, а также от толщины и модуля упругости его стенок. В трубопроводах гидроэлектростанций обычно а = 700-4-1100 м/с. Для этих условий напорный трубопровод'турбины описывается известной системой дифференциальных уравнений в частных про: изводных:

ди __

дИ

дН___ а- ди

(4.34)

' д Г ~ ё

~дТ'

~ m ~ ~ g " d l '

 

где у и Н — скорость течения воды и напор в _произвольной точке, являющиеся функциями времени t и координаты I вдоль трубопровода.

Решение системы (4.34) приводит к следующему выражению [6, 10]:

А = -

А

Л

1

~Тгр

(4.35)

-—

.

 

 

 

1 +

е гР

 

где hw = T J T r — ударная

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

трубопровода;

Тг = 2Ыа — фаза гидравлического

удара, или

время пробега

ударной волной двойной

длины

трубопровода.

 

В (4.35) показательную функцию можно выразить через гипер­

болическую

-V

в форме упру­

th 0,5Тгр. В результате

гого удара

1 + е -тгР

 

передаточная функция Wy будет

 

 

Wy = - 2 h wthO,5Trp.

(4.36)

Таким образом, согласно выражениям (4.33) и (4.36), вид пере­ даточной функции гидроагрегата (4.24) зависит от того, в какой форме будет принята передаточная функция гидроудара. Поэтому прежде всего необходимо выяснить, можно ли вообще и, если можно, то при каких условиях считать гидроудар жестким.

194

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