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

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

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

Логарифмическая фазовая частотная характеристика неустой­ чивого апериодического звена первого порядка в отличие от устой­

чивого звена

изменяется

от ср — — 180° при со = 0 до ср =

—90°

при со = оо

(табл. 2, звено 12).

одного

Уместно

рассмотреть

частотные характеристики еще

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

торого определяется

как .

 

 

 

 

 

w = T ^ f i ’

 

 

(2.116)

откуда

 

к

 

кТ<а

 

W (ко) = •

• Т і 'с о и ■

ѵ =

(2.117)

Г 2со2- I - 1

Г2со3 + 1

Следовательно, амплитудно-фазовая характеристика пред­ ставляет собой полуокружность, которая расположена в первом

квадранте комплексной плоскости. Из (2.117) имеем

 

 

А = r . k

<p = —360° +

arctg7V

(2.118)

]Лг2ша + Г

I

S

\

f

Отсюда видно, что выражение для А совпадает с (2.48.) Зна­ чит, логарифмическая амплитудная частотная характеристика

звена

/г/( 1 — Тр)

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

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

апериодического звена первого

порядка,

а фазовая

изменяется

от ср = —360° при

со = 0 до

ср = —270°

при

со =

оо (табл. 2

звено

13).

 

 

 

второго порядка мо­

2.

Неустойчивые апериодические звенья

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

T\p2~ T lP-f l F

^

 

 

k

 

(2.119)

 

Ңр^ + т ^ - 1 ’

U73 =

_____ k_____

!

 

 

T\P2Txp —-1

 

причем в выражения W2 и

W3 соотношения между постоянными

времени Т г и Т 2— любые.

 

 

 

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

5 В . А . Пивоваров

65

Для передаточной

функции

W1получаем:

 

ГДгсо) =

к_______

 

/г(1 - 7 f o 2 + 7yto)

 

 

~

(1 — Т\О)“)2 + Т2ш2

—Г2со2 — Гр'со + 1

 

* 0 - Ф

8)

.

( 2. 120)

 

( і _ Г 2М2)2 +

Г2Со2

 

I/

______кТха>______

 

 

 

( 1 - Г 20)2)2 +

Г2ш2

J

 

 

 

 

 

чему соответствует кривая в табл. 2 (звено 14), которая является зеркальным отображением относительно оси абсцисс амплитудно­ фазовой характеристики звена 3, поскольку в данном случае при всех значениях со >> 0 V — положительная величина.

Чтобы построить логарифмические частотные характеристики, разложим знаменатель передаточной функции на сомножители:

^ 1== (Тяр — 1) (Г.,р — 1) ’

1Г(К° ) =

 

где Т 3 и Т4 определяются из (2.61). В результате имеем:

 

]Л г2со2+

і]Л г2ш2+ 1

(2.121)

Ф = —360° + arctg Т3со arctg Т4со.

 

Выражение для А совпадает с (2.67) и, следовательно, логариф­ мическая амплитудная характеристика L = 20 lg А будет соот­ ветствовать (2.69), а фазовая частотная определяется как сумма фазовых характеристик двух последовательно соединенных не­ устойчивых апериодических звеньев, первого порядка (2.113).

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

и W3. Так, для звена 2

W 2

(гео) =

k

 

k ( 1 +

Т 2 с о 2 +

Т і < о )

— т\CD2 +

Г р ш — 1

(1 +

Г2со2)2 +

г У ’

откуда

 

 

 

 

 

 

 

 

 

U = —

k ( l + 7 > 2)

 

 

 

 

(1 -р TjCcr)" -j- Tjo

 

 

 

 

 

 

(2. 122)

 

 

 

кТхы

 

 

 

 

Ѵ = —

 

 

 

 

 

(l + Г;со2)2

г2ш

 

 

66

а для звена 3

 

К

 

(i+ryf+n^

 

Г 3(ко)

-Т'Ісо2 — Tjico — 1

 

 

 

 

 

 

 

 

k (1 + Ңш2)

v =

kT^

(2.123)

 

(1 + т у у + т***;

(і + Г=со2)2+П Ѵ

 

 

Кривые, построенные по (2.122) и (2.123), являются зеркаль­ ным отображением друг друга относительно отрицательной дей­

ствительной оси комплекс­

 

 

 

 

 

 

ной

плоскости.

Однако

в

 

 

 

 

 

 

данном случае амплитудно-

 

 

 

 

 

 

фазовые

 

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

 

 

 

 

 

 

(2.122)

и (2.123)

при всех

 

 

 

 

 

 

значениях со не пересекают

 

 

 

 

 

 

мнимую

 

ось. Более того,

 

 

 

 

 

 

при

со =

оо значение ср =

 

 

 

 

 

 

= — 180°, так как при уве­

 

 

 

 

 

 

личении со величина

V бы­

 

 

 

 

 

 

стрее

стремится

к

 

нулю,

 

 

 

 

 

 

чем

U.

Особенно

хорошо

 

 

 

 

 

 

это

видно

 

из

рассмотре­

 

 

 

 

 

 

ния

логарифмических ча­

 

 

 

 

 

 

стотных

 

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

 

 

 

 

 

 

Согласно

(2.119),

зна­

 

 

 

 

 

 

менатели

 

передаточных

 

 

 

 

 

 

функций можно разложить

 

 

 

 

 

 

на одни

и

те

же

сомно­

 

 

 

 

 

 

жители

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wa~ Г 3 -

 

 

 

 

 

 

 

 

 

~ (Т’зР + 1) іР— 1)’ (2-124)

 

 

 

 

 

 

в котором

для

W2 время

 

 

 

 

 

 

Т4 >

Т 3, а для W3, наобо­

 

 

 

 

 

 

рот, Т 3 >

Т4. Выражение

Рис. 27.

Логарифмические

частотные харак­

(2.124)

 

показывает,

что

теристики

неустойчивых

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

звенья

 

с

передаточными

 

звеньев

второго

порядка

функциями W2и W3можно

 

соединением

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

представить

последовательным

звеньев

первого

порядка,

причем

одно из

них — устойчивое, а

второе — неустойчивое.

Логарифмические

частотные характери­

стики

этих звеньев

уже

известны. На рис. 27

сплошными ли­

ниями

показан

результат

сложения

характеристик, входящих

в (2.124), при замене р на

ісо (табл. 2,

звенья 15,

 

16).

5*

67

В частном случае, когда Т 3 = Т4, звенья W2 и W3 превраща­ ются в неустойчивое консервативное звено, для которого

W = ftp*.

№ (до) = —

к

Г-со2 — 1

и =

k

 

(2.125)

ѵ =

о.

ТЧо- - f 1 ’

Из (2.125) видно, что амплитудно-фазовая частотная характе­ ристика представляет собой прямую, расположенную на отрица­ тельной действительной оси комплексной плоскости. Здесь при всех значениях со сдвиг фазы ср = — 180° (табл. 2, звено 17).

3. Неустойчивое колебательное звено описывается дифферен­ циальным уравнением

T, £

^ _

2 t T d x ^ + xBm = kxBx при

£ < і.

(2.126)

Переходная

функция этого звена для хвх =

[1]

будет

иметь

вид

 

 

 

 

 

 

 

Дых

^

-L /

cos <о0г! •

■slncö0£ \j.

(2.127)

 

 

 

 

(

ТсОп

 

 

 

Здесь вещественная часть сопряженных комплексных

корней

г; ± і V і—£2

— положительная

величина,

что

указы­

г = ------- J

 

вает на неустойчивый колебательный характер переходной функ­

ции. Из

(2.126) получаем

 

 

W = -

к

 

£ (1_Г 2м2 + 2£Г;ш) .

ТгР~- 2£7ф + 1 ; W ( i © ) =

(1 — Г2со2)2 ф-4£2Г2о/ ’

 

6(1 — Г-üfi)

 

(2.128)

U =

V

. 2/г£Гсо

( 1 — Т 2ш2) 2 - ( - 4 £ 2Г 2со2 ’

(1 — 7 2cu2)2 + 4£2T22 ‘

Сравнивая (2.128) c (2.64) видим, что эти кривые отличаются только знаком перед V. Значит, на комплексной плоскости кри­ вая, построенная по (2.128), является зеркальным отображением кривой, соответствующей (2.64), относительно действительной оси.

Модуль выражения- (2.128) совпадает с (2.71), а выражение для фазового сдвига ср отличается от (2.72) только знаком, т. е. логарифмические амплитудные частотные характеристики устой­ чивого и неустойчивого колебательных звеньев имеют одинако­ вый вид, а логарифмическая фазовая частотная характеристика неустойчивого колебательного звена изменяется от ср = -—360° при со = 0 до ср = —180° при со = оо. Для определения точных логарифмических характеристик при различных значениях коэф­ фициента следует использовать графики на рис. 23 и 24. Дина­ мические характеристики неустойчивого колебательного звена приведены в табл. 2 (звено 18).

68

Инверторные звенья. К инверторным относятся звенья, кото­ рые изменяют лишь знак входного сигнала или знак одного из его составляющих.

1. Простое инверторное звено определяется алгебраическим

выражением

 

хвых = — кхвх,

(2.129)

которое отличается от (2.36) только знаком, поэтому его динами­ ческие характеристики определяются достаточно просто. Они показаны в табл. 2 (звено 19). В данном случае при всех частотах

от со = 0 до со =

оо сдвиг фазы

—— 180°.

звено.

2. Инверторное

пропорционально-дифференцирующее

Динамика этого звена по аналогии

с (2.99) описывается

одним

из дифференциальных уравнений

вида:

 

 

 

■Дых

^ (kХвх

Т

^вх

\ .

(2.130)

 

dt

)

 

Дых =

( rp dxBx

Д X^ •

(2.131)

 

^ V

dt

Переходные функции для уравнений (2.130) и (2.131) можно получить простой заменой знака соответствующей составляющей на графике в табл. 2 (звено 9).

Рассмотрим частотные характеристики. Согласно (2.130), имеем:

W = k{\ —Тру,

W(/со) — k(l — Тгсо);

(2.132)

'£/ = *;

V = кТа,

 

а из (2.131)

 

 

W = k(T p — iy, ,

W (гео) = k(Tiw— 1);

(2.133)

U = к,

V — kTm.

 

На комплексной плоскости амплитудно-фазовые характери­ стики (2.132) и (2.133) являются прямыми, параллельными мни­

мой оси,

причем

одна из них (2.132) расположена

в

четвертом,

а другая

(2.133) — во втором квадрантах. Модули

комплексных

величин

(2.132)

и (2.133) имеют одинаковое выражение А =

=

k ] /Т 2ю2 + 1 ,

которое совпадает с (2.101),

что

указывает

на

идентичность

логарифмических амплитудных

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

обычного и инверторных пропорционально-дифференцирующих

звеньев.

Сдвиг

фазы для

звена W — k (1 — Тр) определяется

как ср =

—arctg Гео,

т. е.

изменяется

от ф = 0 при

со = 0 до

Ф = .—90° при

со =

оо, а

для звена

W — k ( T p -— 1)

значение

ср = —180° — arctg Т со, т. е. изменяется от ф = — 180° при со = 0 до ф = —270° при со — оо.

69

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

10. Критерии устойчивости

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

Имеется несколько критериев устойчивости. Ниже рассма­ триваются без доказательства основные из них: алгебраический критерий Рауса и Гурвица, критерий Михайлова и частотные критерии, основанные на амплитудно-фазовых и логарифмиче­ ских частотных характеристиках [27].

Алгебраический критерий устойчивости Рауса и Гурвица.

Прежде всего заметим, что для устойчивости системы необходимо (но еще недостаточно), чтобы все коэффициенты характеристи­ ческого уравнения (2.7) были положительны и не равны нулю. Если это условие не выполнено, то система заведомо неустойчива. Указанное условие достаточно для систем не выше второй сте­ пени. Для уравнений выше второго порядка требуются дополни­ тельные соотношения Рауса и Гурвица, которые можно получить следующим образом.

Пусть имеется характеристическое уравнение системы

a0z" + а1гп~1 + ----- 1- an_xz + аѣ= 0.

(2.134)

70

Из коэффициентов этого уравнения составляется специальная матрица

«1

Йо

0

0

а со

й -2

й і

й 0

й 5

. .

0

0

0

Й 4

. .

0

0

0

 

 

 

 

а з

. .

0

0

0

Й о . . .

0

0

0

0

0

0

.

■ • й „ _ з

« п - 1

0

0

0

0

. • • й ; , _ 4

й / j—2

й „

по правилу: в первую строку выписываются все коэффициенты

характеристического уравнения с нечетными номерами и в конце

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

ментов было равно степени п характеристического уравнения,

во вторую строку выписываются коэффициенты с четными но­

мерами, начиная с й0. Затем первая и вторая строки сдвигаются

вправо на один

номер и т. д., причем

на

место коэффициентов

с порядковыми

номерами, большими п

и

меньшими нуля, ста­

вятся нули. Всего должно быть написано п строк. По этой ма­ трице составляется п определителей Рауса и Гурвица

йі а3

Ді — ах\ Л2=

ао а2 ’

откуда

II со <

ЙІ й3 йв й0 й2 а 1

0 ЙІ й3

Ді = йь Ао = Й1Й2 — Оойз! Аз = йз(йіЙ2 — йойз) —

ах(айаъ—oyaj и т. д.

(2.136)

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

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

71

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

Критерий устойчивости Михайлова. Кр'итерий Михайлова основан на свойстве годографа характеристического уравнения, известного под названием кривой Михайлова. Чтобы получить кривую Михайлова, нужно в характеристическом уравнении си­ стемы (2.134) заменить г на чистое мнимое значение ісо, как это делалось для определения частотных характеристик. Тогда полу­ чим многочлен

L (ісо) = а0(«о)" + % (ко)'1" 1 - f ----- j- ап_г(гео) + ап,

(2.137)

в котором можно выделить вещественную и мнимую части

L (too) — U {со) -f іѴ (со),

(2.138)

где

а„_4со4 -

 

U (со) = ап— а„_2со2 +

(2.139)

К (со) = ап_хсо — а„_3со3 +

а„_5со5

 

Для каждого значения параметра со величину L (ісо) можно изобразить на комплексной плоскости в виде вектора. При изме­ нении параметра со в интер­

вале

от

со = 0

до

со

= ОО

конец

этого

вектора

опи­

шет кривую, которая и назы­

вается

 

кривой

Михайлова

(рис.

28, й).

Если

условия

устойчивости

выполняются,

т. е. если вещественные части

всех корней

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

ского

уравнения отрицатель­

ны, эта кривая обладает сле­

дующими

свойствами.

 

1) Начало кривой (со = 0) всегда расположено на веществен­

ной оси на расстоянии, равном аѣ от начала координат,

а ее ко­

нец (со = оо) всегда уходит в бесконечность. Это следует из вы­ ражения (2.139).

2) Если 11 — степень характеристического уравнения, то кри­ вая Михайлова пересекает вещественную ось пі2 раз при четной п и (п + 1)/2—при нечетной /г. Мнимая ось пересекается кривой п/2 раз при четной п и (/г—1)/2 раз при нечетной п. Значения со, при которых кривая пересекается с вещественной и мнимой осями,

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

V (со) =

0 и

U (со) = 0.

 

регу­

Критерий Михайлова гласит следующее: чтобы система

лирования п-й степени была устойчивой, необходимо и доста­ точно, чтобы кривая Михайлова, построенная для характеристи­ ческого уравнения данной системы, поочередно проходила п квадрантов против часовой стрелки, все время окружая начало

72

координат (рис. 28, б, кривая /). Когда система неустойчива, кривая Михайлова минует какие-либо квадранты комплексной плоскости (рис. 28, б, кривая 2).

Границей устойчивости является наличие нулевого

корня

2 = 0 или чисто мнимых корней zli2= ± m . Если один

из кор­

ней характеристического уравнения z = 0, то это означает, что свободный член ап = 0. В этом случае кривая Михайлова будет

исходить

из начала

координат, так как

по

(2.139)

при ш = 0

U (ш) =

0 и V (со) =

0 (кривая 3 на рис. 28, б). Когда же имеется

пара чисто мнимых

корней zli2= ±/co,

то

кривая

проходит

через начало координат при некотором

значении со

= со 0, по­

скольку

наличие корня

z = /со эквивалентно

приравниванию

нулю выражения для кривой Михайлова,

т. е.

L (ісо) = 0 или

U (со) =

0 и V (со) =

0, причем оба эти выражения должны удов­

летворяться одновременно

при

одном и

том же значении со =

= со0 (кривая 4 на

рис. 28, б).

Значение

со

= со0, которое соот­

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

находящейся на границе

устойчивости.^

П р и м е р ы п р и м е н е н и я к р и т е р и я М и х а й л о в а .

1. Допустим, что получено

характеристическое уравнение системы в виде

 

2г3 +

Юг2 +

5 г + 8 = 0.

Заменим г на

г'ш

 

 

 

—2/со3 — Юсо2 +

5і'со - { - 8 = 0 .

В результате

получим:

 

 

 

U (со) = — Юсо2 + 8;

V (со) = —2со3 5со.

 

Определим точки пересечения кривой с осями1U и V. Ось U пересекается

кривой

при

V (со) =

0, т. е. когда —2со3 +

5со = 0,

откуда сл1 = 0

и

со2 =

= ч-

 

=

±1,58 . Поскольку кривая Михайлова

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

только

при

положительных

значениях

со, то

со2 =

1,58, при

со =

0

U (0) ,=

8; при

со =

1,58 U (1)58) =

— 17.

 

 

 

 

 

 

V при

 

Найдем

теперь

значения со, при которых кривая

пересекает ось

U (со) =

0,

т. е.

когда — Юсо2 +

8 =

0, откуда со =

0,9.

Для

со =

0,9

V (0,9) = —2- 0,93 + 5- 0,9 = 3,05.

На рис. 29, а показана кривая, построенная для различных значений со. Данная система устойчива, так как кривая на рис. 29, а удовлетворяет крите­ рию Михайлова.

2. Построим кривую Михайлова для системы пятой степени с характери­ стическим уравнением

5гв + 4 zi + Юг3 + 7га + 5г + 9 = 0.

После замены г на гсо имеем:

 

 

U (со) =

4со4 — 7со2 +

9;

V (со) =

5со8 — 15со3 +

5со.

Определяем точки пересечения с осями U и V. Кривая пересекает ось U, когда У (со) = 0, т. е. при

5со6 — 15со3 + 5со = 0,

откуда для положительных частот сох = 0, со2= 0,635, со3= 1,62.

73

Для этих Wзначения U (со) будут: U (0) = 9; U (0,635) = 6,8; U (1,62) = 18.

Далее, кривая пересекает ось V, когда U (со) = 0, т. е. при 4со4— 7со2+

+ 9= 0.

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

построенной по точкам. Она сразу проходит из первого в четвертый квадрант, минуя второй и третий.

Частотные критерии устойчивсти. Рассмотренные выше кри-

терии устойчивости применимы лишь в тех

случаях,

когда из­

 

 

 

 

 

 

вестно характеристическое

 

 

 

 

 

 

уравнение

системы

или,

 

 

 

 

 

 

что то же самое,

 

когда по-

 

 

 

 

 

 

лучены дифференциалы-іые

 

 

 

 

 

 

уравнения

всех

 

звеньев

 

 

 

 

 

 

входящих

 

в данную систе­

 

 

 

 

 

 

му. На практике эти

урав­

 

 

 

 

 

 

нения

не

всегда

удается

 

 

 

 

 

 

составить,

однако

имеется

 

 

 

 

 

 

возможность (особенно

в

 

 

 

 

 

 

электрических

 

' схемах)

 

 

 

 

 

 

собрать

макет

 

 

звена

и

 

 

 

 

 

 

опытным

путем

 

получить

 

 

 

 

 

 

его

частотные

характери­

8

10

П

іб

18 U

стики.

В

таких

 

случаях

для

определения

 

устойчи­

Рис. 29. Кривые Михайлова

 

для устойчи­

вости системы использует­

 

ся

частотный

критерий,

вой (о)

и неустойчивой

(б) систем

 

 

 

 

 

 

позволяющий по частотной

 

 

 

 

 

 

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

разомкну-

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

 

 

Аналогично

(2.12)

передаточная

функция

разомкнутой си-

стемы

 

 

 

 

Д (Р)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 ( Р )

,

 

 

 

 

 

 

(2.140)

 

 

 

 

 

Q(р)

 

 

 

 

 

 

 

 

 

 

где R (р) и Q (р) — правая и

левая

части

дифференциального

уравнения

разомкнутой

системы.

в правую часть уравнения

Условимся не вводить знак

минус

(2.140), как

это делалось для

замкнутой системы.

 

 

 

 

 

 

Формулируемые ниже частотные критерии справедливы для - систем, у которых степень числителя R (р) передаточной функции G (р) меньше или равна степени 'знаменателя Q (р). Это условие выполняется у большинства систем, в том числе и у систем регу­ лирования гидротурбин.

Чтобы определить устойчивость системы, необходимо построить на комплексной плоскости амплитудно-фазовую частотную ха­ рактеристику разомкнутой системы регулирования G (ш). На этой плоскости отметим точку С с координатами (—1, і-0), ко­ торая называется критической точкой.

74

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