![](/user_photo/_userpic.png)
книги из ГПНТБ / Пивоваров В.А. Проектирование и расчет систем регулирования гидротурбин
.pdfгде индексом п обозначены номинальные (расчетные) значения соот
ветствующих величин. |
Разделим выражение (1.2) на Міп и, учи- |
|||||
da> |
= CÖ„ |
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