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

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

.pdf
Скачиваний:
7
Добавлен:
20.10.2023
Размер:
8.39 Mб
Скачать

10.2. Добавочные потерн, вызванные пульсацией магнитного потока

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

При расчете потерь в массивном роторе турбогенератора можно также применить следующую приближенную формулу.

Использовав формулу (8. 5) с учетом того, что в рассматривае-

мом случае

/ ѵ

=

— Д,

% =

— ,

где % — число

пазов

статора,

для

 

 

 

 

 

Р

 

zi

 

 

 

 

 

 

 

 

 

 

 

 

расчета

 

потерь

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(g0 /cr v )a

ДЗ

 

// . \1 . 5

 

 

 

 

,,,ПІ9\

 

 

 

 

Сы = Ъ

^ -

а

і

{

£ )

Ä„„C e I

D T ,

 

(10-12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7l2.55Ql.57

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I I " " ' U V / - * - ' " J

 

 

 

 

 

 

 

 

 

 

 

 

 

Г Д Ѳ

 

 

f p =

4 \ / ^

'

( 1 0 - 3 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I i i

+

iîi

10.

4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и учитывает отклонения деистви-

 

 

0.3

 

OA

0.5

 

O.B

тельной

формы

пульсации

 

поля

 

 

 

 

ôn/t^d.e

 

 

 

от

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

синусоиды

 

 

 

 

 

 

 

 

 

(рис. 10-2).

 

 

 

 

 

 

Рпс. 10-2.

Зависимость frn00 отbjt^

 

 

 

Коэффициент Са

в

(10. 2)

учи­

 

 

 

 

 

 

 

 

 

тывает увеличение индукции из-за

зубчатого

строения

ротора.

Зависимость

Ca=f

(b.2/^~)

приведена

на рис.

8-8.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Б0

В современных мощных синхронных машинах амплитуда

 

составляет

порядка

1 0 0 - Ю - 8

вб/см2 .

С

учетом

значения

 

 

BScv

имеем f i Ä d 5 0 0 - 1 0 8

гн/см.

При

р = 2 - 1 0 " 5

ом-см,

fi=1500x

Х Ю " 8

гн/см и величинах воздушного зазора S ^> 6.0 см коэффициент

реакции

Ä r v = 0 . 9 5 . При обычно

принимаемых

соотношениях

 

раз­

меров

пазов

ротора

и

статора

 

п0та 1.1,

С^таі.О.

Поэтому

 

 

 

 

 

!2z0 :

 

 

 

Б'п

Я?

/ / л и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В отечественной расчетной практике используется формула, подобная формуле (10. 5). Однако численный коэффициент в фор­ муле принимается равным 5 . 1 - 10 7 . Такая величина коэффициента будет иметь место при fx=2000• Ю - 8 гн/см.

Аналогичные потери в шихтованном полюсе гидрогенератора равны

160

В

формуле (10, 6)

принято р д = 1 6 0 0 - 1 0 - 8 гн/см,

р = 2 х

Х І О - 5

ом-см, к„т=1.2,

Са=1Л, &о 6 р =1.3. Для определения

потерь

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

значениях и форме

записи будет иметь

вид

 

02,0 =

13 • 10ЧРУ^a^ltBl(tÇj'\

вт.

(10.7)

Формула (10.7) не учитывает особенностей проникновения маг­ нитного поля в расслоенную сталь, и результаты расчета потерь по этой формуле оказываются завышенными в 1.2^ раза по срав­ нению с результатами расчета по формуле (10. 5), которая, как показали экспериментальные исследования, дает результаты, близкие к опытным данным.

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

При этом не принимается во внимание влияние ширины зубца статора, хотя и приводятся рассуждения о том, что при

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

пульсации в зубцах статора наблюдаются и при указанном выше соотношении размеров.

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

насыщения при расчете индукции в зубце.

 

Пульсацию индукции в зубце і?р более

правильно рассчиты­

вать через пульсацию

потока

 

 

 

Дрі =

До2Г77

тЛ

(10.8)

 

 

1

'эфф

'ь 0 л>р

 

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

в (10. 8) значения индукции в среднем се­

чении зубца

статора

 

 

 

 

 

5 p l

= B

0 2 - ^ ^

s i n - ^ .

(10.9)

Потери в

зубце статора, вызванные Вв1,

могут быть найдены

по формулам (9. 12) или (9. 13), а также (9. 16), если вместо і?ѵ подставить в эти формулы Вп, согласно (10. 8) или (10. 9), и учесть синусоидальное распределение основного потока в зазоре с помо­ щью коэффициента 1/2.

1/2 11 Я. Б. Данплевич

З а к л ю ч е н и е

Характерной особенностью турбо- и

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

большой мощности

является

быстрое увеличение линейной токо­

вой нагрузки при

мало

изменяющихся

магнитной

индукции

в зазоре и габаритах машин.

Одновременно с

ростом

мощности

турбо- и гидрогенераторов увеличивается

воздушный

зазор.

В результате в машинах большой мощности значительно уве­

личиваются добавочные

потери в торцовой

зоне генераторов

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

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

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

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

тельных величин. Так,

для турбогенератора мощностью 800 Мвт-

в режиме

номинальной

нагрузки потери составляют в среднем

173 вт/см,

в отдельных

точках увеличиваясь до 400—550 вт/см.

162

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

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

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

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

Исследования влияния транспозиции на потери от циркуля­ ционных токов показали, что применение транспозиции на 540° в пазовой части обмотки позволяет существенно снизить добавоч­ ные потери в стержне от циркуляционных токов, наведенных в эле­ ментарных проводниках внешними магнитными полями. Для существенного снижения и токов, вызванных собственными по­ лями рассеяния, необходимо также транспонировать элементар­ ные проводники в каждой лобовой части. Пазовая часть стержня при этом может быть транспонирована на угол 360°. Исследования показали, что наиболее рационально выполнять транспозицию лобовой части стержня таким образом, чтобы каждая лобовая часть содержала три участка: два равных по длине нетраиспонированных (один на выходе из паза), а третий — транспонирован­ ный на 90° участок и расположенный (у головки) между нетранспонированными участками.

11*

П р и л о ж е н и е 1

Связь между S* и основным магнитным полем

Еолп представить

 

*-*+*. !

, Ш Л )

и предположить, что сложение основного и аксиального магнитного полей происходит в среде с и. = (і0 , то уравнения электромагнитного ноля для Й0, Нв, Ê0 п ÊB запишутся в виде

rot Н0 =

0,

rot # „ =

-(£,

)

*

\

'*

( Ш - 2 )

Из ( Ш . 2) для Ь» имеем

 

 

 

 

В активной зоне генератора существует только составляющая Ел (токи статора и ротора имеют только осевую составляющую), поэтому

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

 

 

it

ff, =

2 t f J , V W

( ш - 5 )

 

 

 

 

 

 

 

где Аф =

^,- — коэффициент

формы

для синусоидальной

кривой поля,

1{ — активная длина машины.

 

 

 

 

Д л я

амплитуды напряженности поля

имеем

 

 

 

 

 

U,

 

 

 

 

 

E!0m

=

—rf-

=

2-.fBym,

( Ш . 6 )

откуда

из

( Ш . 4 ) получим

 

 

 

 

 

164

а ^ = - і * В у т .

( Ш . 7 )

 

 

П р и л о ж е н и е

2

 

 

 

 

 

 

 

 

 

 

 

Вывод уравнений (2. 24)

 

 

 

 

 

 

 

 

П р о ек т и руя

уравнения

поля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B =

votÄ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е = —jwA -

grad у, 1

 

 

 

 

 

где — скалярный потенциал,

получим

 

 

 

 

 

 

 

по

оси а:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—J

Ах -

у д£

=

тоЬх {(р . ^ 1 ) rot À)=Yy

 

{ . A _ 1 rot Â)„

-

 

 

 

à

A Л

 

 

b [дАу

 

дАх\

 

j

д

/дАЛ

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A n

0

0

 

 

 

 

 

 

 

 

 

 

 

 

л J =

 

0

k2i

0

— обратный

тензор,

k u

= —

k

k

L

 

1

0

0

k33

'

 

 

 

 

 

^

 

 

^

 

^

 

 

 

 

 

ö(/2 л 22 Ö Z 2

da;

( *

3 3 ^

+

^

=

^ г ;

 

(2П.2)

 

 

 

 

 

V

ду

1

ç J

p

 

 

по

оси I/

 

 

 

 

 

 

 

 

 

 

 

 

 

- f ^ - I S = r o t * < ^ r o t ^ = 5- ^

r o t i } * - Й ^

r o t À

)• =

 

 

 

 

 

1 1 öz V ду

dz J

d J

da: \

дх

ду ) '

 

 

 

 

 

кп

+

к33

d^J>-±

33дАі+Г)=!?

p J

p

Al.

 

(2П. 3)

 

 

 

 

1 1

Öz2 ^

3 3

öa;2

ду \

3 3 дх ^

 

y

 

 

Т а к

ка к div A = 0,

то

 

 

 

 

 

 

 

 

 

 

и в результате получаем

следующие уравнения для 4 Х и . < 4 у ( ^ = 0)

(см. выражение 2. 24)

 

 

 

 

 

 

£ 4 * 4 -

dJA^,

1

â*Ax

іщ>.я

 

дх* "»" <ty2 t -

^ / ^

Ö Z 2

р

X,

, ^ £

,

 

дг2

_±_д=М^А

дх*

<ty2

"""{^/р.,

р

у

165

П р и л о ж е н и е

3

 

 

Формулы для расчета плотностей вихревых

 

 

 

 

 

токов

в

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

области

 

 

 

 

 

 

 

 

 

Д л я

сокращения

времени

счета

на

ЭЦВМ

при

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

числен­

ного метода определения функции ае(х,

у)

(см. раздел

2.6) за

начальное

ноле

значений

этой

функции

принималось аналитическое решение урав ­

нения (2. 34)

для

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

области,

ограниченной

стенками

паза

и

шлица . В этом

случае

решение

(2. 34)

 

представляется

в

виде ряда

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ыВп /ch

la

— 1

 

 

 

 

 

\

 

 

 

 

 

 

 

 

 

 

 

 

 

2Т^2 " Ч

shla

 

 

s h

l

x ~

c h

l x + 1

J s

m

І

У г

 

 

( З П > *)

где

Ь

 

,

В„

=

I Bj{y)

sin lydy

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

разложения а к с и а л ь -

 

 

Ъ

 

 

 

Ъ оJ

 

 

 

 

 

 

 

 

 

по sin ly,

 

hj

 

 

 

 

ной

составляющей

индукции

в

ряд

Фурье

 

— расстояние

от

коронки

зубца

до

нижнего

к р а я нажимной

плиты,

b =

h-\-ha

— высота

сегмента

статора,

а— ^

 

половина

от

ширины зуоца.

 

 

 

 

 

 

Из выражения (ЗП. 1) можно

получить

расчетные

 

формулы

для токов

и потерь

в

области

коронки

зубца,

 

причем

потери

 

на

левом

и

правом

к р а я х коронки

зубца при

і =

0

в

а; =

 

а равны между

собой:

 

 

 

 

 

 

 

 

 

Ѵ *-«, -

 

Z l

? Л

 

sh la

J\ : sin

ly

 

 

 

 

(ЗП.

2)

 

 

 

 

 

 

 

.1=0

 

 

L „=i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a плотность

тока

имеет

следующие

составляющие:

 

 

 

 

 

 

 

 

 

 

 

 

••—

p

^

> В„Гс\\1а

 

1sh Ix — ch Ix -(- 1

cos

ly,

 

 

 

 

 

 

 

я=1

I

L

shla

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ЗП.

3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ p I

 

ch la

— 1• ch Ix

— sh

Ix

 

sin

ly.

 

 

 

 

 

 

 

> у т _

 

 

sh la

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П р и л о ж е н и е

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формулы

для расчета

составляющих

поля

 

 

 

 

 

 

в случае немагнитного внешнего щита

 

 

 

 

 

 

 

Рассмотрим к-іі элемент

с током (рис.4П-1). Тогда граничные

условия

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

составляющих плотности тока

будут такими, как это показано на рис. 4П-1.

 

 

 

 

 

 

 

 

 

Найдем составляющую векторного потенциала Ахк,

 

 

другие

составляю­

щие

ук

 

и

Агк)

 

могут

быть

найдены

 

аналогичным

способом.

 

 

 

 

Представляя Ахк

в в и д е *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4П.

1)

 

* В

дальнейшем

индекс

к

оп}'скается.

 

 

 

 

 

 

 

 

 

 

166

д ля Uх будем иметь

diß "т" dz'

Н-05а.-

(4П. 2)

Представим Ux в виде

пку

(4П. 3)

ь п х cos ß ,

 

71=0

ay

 

 

 

 

 

 

 

 

 

 

z

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3y = 0

 

Рис . 4П-1. К расчету поля в торцовой зоне генератора для слу­

 

 

 

чая, когда внешний щит выполнен немагнитным.

 

а — гранпчпые условия для Sx;

б — то же для Sy, о — то же для S z .

 

Коэффициенты ряда

Фурье

в (4П . З),

согласно

[12, 16], равны

 

 

 

 

 

 

 

22

в

 

ппу

 

 

 

 

 

 

Ѵ-дСпх

г

 

dy.

 

 

(4П. 4)

 

 

ТJ ^xcos -jr

 

 

 

 

 

H

' H

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Продифференцировав

 

(4П.4) no z как по параметру

 

 

 

(x0

d*Cnx_

2

[dW?dWx x

nnyпъу

 

 

(4П. 5)

 

 

H

 

dz°- ~II

J dz

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

и умножив

обе части уравнения

 

 

2

пщ)

 

(4П. 2) на ß-cos—jj-

аУ> проинтегрируем

полученные

соотношения

 

в пределах

от

0 до Н. Тогда с учетом

(4П. 5)

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

lx

d2C„• я » . 22гд*иtd*Uх

x

___mtyп%у

/fnyt4>„и \ И о -

 

2ц-j-o

 

0

 

x

 

 

 

 

 

а

 

 

0

(-6)

H

 

dz 2~ + F J ~dyTcos~Wd'J-[j)

 

Wc^~-TrJ"x'

167-

где

 

 

Jn*=

\ S*

cos -g-dy

 

 

= — 2 5 '

» Я

пк

 

nah.

[z—

 

cos » —

cos г - A J - sin

0

1

-

I

 

,

w

it>it s

i n

о и sm

 

 

птс

H

K

Iii

\

ck

 

 

 

•z— Ак

 

,

• ^ 0 x = — S'mK COS » S i n

 

 

 

 

 

-

A,. Ѳтс

 

 

-

f0

d - ß ) T

 

2 тс-

значения 5,'„ и других величин приведешь

в гл. 1.

 

 

 

 

 

 

 

Представляя

интеграл

в (4П. 6)

в

виде

 

 

 

 

 

 

 

 

 

 

 

 

2

г

d^Ux

тѵку

 

2 dUv

 

пку'

+

 

 

 

 

 

 

я

J

-dyTC0Sirdy

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, Г 2 /ITC

гатсіЛН

 

/гетс\ 2 и п

 

 

 

 

(^7>

 

 

 

 

 

+[я" Я ^s i n irJo " ( я )

я с «

 

 

 

 

и учитывая, что первые два члена

в

правой

части

(4П. 7)

 

равны

нулю,

получим

окончательно

d2Cnx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-JÏT—VCnx

 

= -2J„x,

 

 

 

 

 

 

 

(4П . 8)

Решением (411с „ . =8) оявляется^ * + D

2выражениеB - ^ - Y

J ; „ х sh

 

\{z~t)dt.

 

 

(4П. 9)

Так

как

при

z =

-f-m

величина

Л х

ограничена,

а

при z =

0

^

= 0 ,

то D1 =

ö 2

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D1 =

T

J Jnxz-Udt.

 

 

 

 

 

 

(4П. 10)

Подставляя в

 

(4П. 9) и (4П. 1 0 )

значение

/ п х ,

в

зоне

I

будем

иметь

 

 

4 Я

 

 

 

/п-ккЛ. .

/птсй0 \

 

chXz

 

Г

 

 

тсЗ

 

 

 

 

 

+

Ç s i n | E

- ^ - b - x

^ ) '

 

 

 

 

 

(4П. 1 1 )

 

 

 

 

 

 

 

 

 

TCZ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c h T

 

тс

х " ' с

 

Зтс .

 

 

 

 

С 0 ж =

 

2S'mha

COS & -Г-ГТ2

7 | т 7 \ 2

Е

 

COS-5-

+

 

 

 

 

 

 

, STC

ВТС -:лк

 

тс

-^fc+ofc)

 

 

 

 

 

(4П. 1 2 )

 

 

 

 

+ т З І П 2

 

 

- т Е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

168

П р и л о ж е н и е

6

 

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

 

 

 

по головкам зубцов

 

 

 

 

 

При расчете статор машины принимается

 

гладким, а зубчатость ста­

тора

учитывается

с помощью коэффициента

зазора к^. Д л я

машин боль­

шой

мощности отношение

зубцового шага £2

к

диаметру ротора D2

мало,

и кривизной поверхностей

можно пренебречь.

В результате

для

отдель­

ного

паза ротора

получим

расчетную схему

(рис. 511-1, а).

Рассмотрим

лишь

ее половину, например левую. Расположим ось х плоскости z вдоль

поверхности

ротора, а

ось

j'y— по

осп симметрии паза. Точки

Ал, А2,

А3 и

Л 4 па

плоскости

z

отобразим

точками

Ç2>

1 1 и а

верхней

®

/ / / / / / / / ,

y77777^

V

Рис. 5П - 1 . К расчету проводимости рассеяния по головкам зубцов.

полуплоскости Ç = Ç - j -

/ij

(рис. 5П-1, б). Тогда для преобразующей

функции

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hi)

c

 

 

 

 

 

2ft

arc

t:

1 /fc2Ç

 

 

 

 

 

 

 

 

 

 

 

 

g У

 

 

 

 

 

Определив

постоянные

С

и

/с из соответствия точек А2

и А3

на

плос­

кости г и ?2 и

?з на плоскости

Ç,

д л я

z окончательно будем

иметь

 

 

 

ІЪ»

ІП

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5П.

1)

Поверхности А2А3АІ

и

 

^

4

^

являются

эквипотенциальными для

потоков рассеяния по головкам

зубцов,

причем

поверхность

А2А3А±

имеет

12 Я. Б. Дашшевич

169

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