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

книги из ГПНТБ / Основы автоматизации для металлургов

..pdf
Скачиваний:
13
Добавлен:
23.10.2023
Размер:
33.25 Mб
Скачать

где

X,

Z — нелинейные

функции;

 

 

 

 

 

 

X — отклонение

регулируемой

величины;

 

 

 

z — отклонение

возмущения.

 

 

 

Это уравнение п-ного порядка можно заменить системой я диф­

ференциальных

уравнений

первого

порядка:

 

d.X\

 

V /

 

 

 

 

 

 

 

 

\

 

d t

=

-Л 1

, х 2

,

• • •>

Х п ,

Z \ ,

Z 2 ,

• •>

Zm)>

 

—£f- =

X2(xl,

x2,

 

xn,

Z i , z2,

...,

zm),

(151)

dt " ~~~Xn(xx,

x2,

...,

xn,

zu

z2,

 

...,

zm),

(152)

где

 

 

 

 

 

 

 

 

 

 

 

 

Xx—x,

x2

d t

 

, ...,xn

 

d(ll_l

 

.

 

 

 

Переменным

X\, x2,

xn

придаем

геометрическое

представление,

т. е. считаем их координатами точки в n-мерном пространстве. Тем самым мы получим фазовое пространство, так как каждой точке соответствует определенная фаза или состояние системы. Напри­ мер, для контура, выраженного уравнением второго порядка, по­ лучим фазовую плоскость. При изменении t точка перемещается по кривой, называемой фазовой траекторией, выражающей пове­

дение системы в динамике. Направление и скорость

перемещения

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

'

dxn

 

 

,

 

 

 

 

—^- , называемого фазовой скоростью:

 

 

v =

ix' + jx".

 

 

 

 

 

(153)

В качестве примера рассмотрим нелинейное уравнение второго

порядка:

 

 

 

 

 

 

*" +

£ ( * ' ) + / ( • * ) =

0.

 

 

 

(154)

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

образом:

а. Получим параметрическое выражение уравнений, описываю­

щих

решение,

и

его

первую

производную

[x = x\(t), x' = x2(t)].

После исключения

параметра t

получим уравнение фазовой траек­

тории в переменных х\,

х2.

 

 

 

б. Подстановкой

х = хи х' = х2 получим

вместо

первоначаль­

ного уравнения

систему

уравнений:

 

 

-^-=-g(x2)-f(x1),

 

 

 

 

 

(155)

 

 

 

 

 

 

 

dt

- Л-о t

 

 

 

 

 

 

 

 

 

 

 

 

 

При делении этих выражений получим

 

 

-1хТ~

~2

 

 

 

(156)

166

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

менных Х\, х2. Графическое

построение фазовой

траектории про­

изводится следующим образом: на фазовую плоскость хи

х2 нанесем

кривые хі = g(x2 ) и x2 = f(x]),

представляющие

собой

нелинейные

функции (рис. 105). Заданный угловой коэффициент фазовой тра­

ектории в любой

точке (х°\,

х°2) фазовой

плоскости получим,

про­

водя в этой точке параллельные осям хь

х2 прямые. Затем на па­

раллельной оси

Х\ прямой

отложим от

оси х2 величину

a=J(Xj).

В результате получим точку А. Из точки пересечения прямой, па­

раллельной оси

Xj, с кривой

Хі = g(x2 )

проведем

перпендикуляр

к

оси х\.

Из

полученной

точки

пересечения

этого

перпендикуляра

с

осью

Х\(С)

проведем

соединительную

прямую

с точкой А.

Пер­

пендикуляр

к этой

соединительной

пря­

 

 

 

 

 

 

 

мой

в точке

А

будет

угловым

коэффи­

 

 

 

 

 

 

циентом фазовой траектории в данной

 

 

 

 

 

 

точке.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Метод

фазовой

плоскости

применя­

 

 

 

 

 

 

 

ется и при изучении систем с помощью

 

 

 

 

 

 

 

аналоговых ЭВМ. При этом регистрация

 

 

 

 

 

 

 

процесса

осуществляется

при

помощи

 

 

 

 

 

 

 

двухкоординатного

самописца таким

об­

 

 

 

 

 

 

разом,

что на горизонтальной

оси

запи-

 

 

 

 

I

 

 

сывается

значение

Xi

= x(t),

 

а Н а

В е р -

 

рис,

105.

Графическое

построе-

Т И К а Л Ь Н О Й

З н а ч е н и е

X2

= x'(t).

 

 

 

 

 

н и е

Фазовой

траектории

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так, например, при исследовании дви­

 

 

 

 

 

 

гателя,

используемого

для

 

дистанционной

установки

положения

клапанов, исходим из нелинейного уравнения

равенства

мо­

ментов:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mmot

= Mz+MB+MT,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(157)

где

Mmot

Mmf(x)—момент

 

 

двигателя,

соответствующий макси­

 

 

 

 

 

 

 

мальному

моменту

Мт

 

при

пуске

двигателя

 

 

 

 

 

 

 

в

 

направлении,

заданном

 

функцией

f(x)

 

 

 

 

 

 

 

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

пусковым

реле);

 

 

 

 

 

 

M z = /a — момент

 

нагрузки,

определяемый

произведе­

 

 

 

 

 

 

 

нием величины момента инерции / и ускоре­

 

 

 

 

 

 

 

ния

а =

х"\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мв

= Вѵ — демпфирующий момент, равный произведению

 

 

МТ

= M0signv

 

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

трения

В

и скорости

ѵ = х'\

 

 

— момент

 

трения,

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

нормаль­

 

 

 

 

 

 

 

ному моменту трения М0

и направлению

дви­

 

 

 

 

 

 

 

жения

 

(скорости

ѵ = х')\

 

 

 

 

 

 

 

 

 

 

 

X — угол вращения

(регулируемая

величина).

 

Подстановкой

выражений

для

моментов

M m o t ,

Mz,

Мв

и Мт

в уравнение

(157), получим

уравнение

в операторной форме

(LW),

пригодное для решения на аналоговой ЭВМ:

 

 

 

 

 

р°-х = арх

— b s'ïgnpx — cf(x),

 

 

 

 

 

 

 

 

 

 

(158)

167

где

_ß_

Ma

С

/

 

 

 

Соответствующая блок-схема и схема подключения двухкоординатного самописца показаны на рис. 106.

Самописец 6АК

Рис. 106. Исследование замкнутого контура регулирования в фазовой плоскости при помощи аналогового вычислительного устройства

Устойчивость нелинейных систем

Устойчивость линейных систем зависит только от знаков корней характеристического уравнения, т. е. от свойств самой системы. У нелинейных систем устойчивость зависит не только от свойств системы, но и от начальных условий.

О поведении контура можно судить по форме фазовой траекто­ рии:

1. Неубывающие колебания характеризуются фазовой траекто­

рией

в виде замкнутой кривой (рис. 107, а).

 

 

2.

Убывающие колебания характеризуются траекториями в виде

спирали,

асимптотически

приближающейся

к началу

координат

(рис. 107,

б).

 

 

 

3.

При нарастающих колебаниях траектория по мере увеличе­

ния

времени удаляется

от начала координат

(рис. 107,

в).

4.Апериодический затухающий переходный процесс имеет фа­ зовую траекторию, показанную на рис. 107, г.

5.Апериодический незатухающий переходный процесс показан на рис. 107, д.

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

168

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

моники.

 

 

из линейной части (FL) и не­

Пусть система (рис. 108) состоит

линейной части (Fjv). У

нелинейных

элементов

при синусоидаль­

ном входе х\ выход x2

= N(x\) является

обычно

несинусоидальным,

он зависит не только

от

частоты, но

и

от амплитуды, т. е. x2N' =

= (А, со). Запишем выражение передаточной функции нелинейной

части:

 

 

 

 

 

 

 

 

 

FN(A,

/V): * 2

! t - ' m i

 

 

 

 

 

ѵ (159);

 

хі

(Л, у»)

 

 

 

 

 

 

 

 

 

 

Общая

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

функция

 

 

 

имеет

вид

 

 

 

 

 

 

 

 

 

 

FLU*)FNX

 

 

 

 

 

 

 

Х(А,

»

=

- 1 ,

(160)

 

 

 

 

 

 

FL

 

 

 

 

 

 

 

 

 

F,

 

 

 

Рис. 107.

Ф а з о в а я

траектория

нелинейных Рис.

108.

Замкнутый

контур

регулирова-

 

контуров

 

ния

с

нелинейным

регулятором

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

 

 

 

FL(J*)FN(A,

»

+ 1 = 0 ,

 

 

 

 

 

 

(161)

при решении которого можно найти критическое значение Л и со для момента возникновения колебаний. Это решение часто выпол­

няют

графически

с применением

преобразования

- Ъ

М - Т с а

и ) .

062)

когда обе характеристики представляются в комплексной плоско­ сти (рис. 109). Их точки пересечения M и N определяют амплитуду и частоту возможных колебаний. Аналогично применяется выра­ жение

для решения в логарифмических координатах. Об устойчивости

найденного решения судят по значениям Р

À . . и |— FL(j<a)\

169

при изменении А на величину АЛ. Устойчивое периодическое реше­ ние получается при положительном значении АА, если

•>\-FLU«>)

(164)

 

Теории нелинейных контуров посвящена специальная литера­ тура, обзор которой приведен в ра­ боте [22].

 

 

 

 

 

 

 

 

8.

АГРЕГАТНАЯ

СИСТЕМА

 

 

 

 

 

 

 

 

 

 

РЕГУЛИРОВАНИЯ

УРС

 

 

 

 

 

 

 

 

 

Унификация

элементов

систем

 

 

 

 

 

 

 

управления дает не только значи­

 

 

 

 

 

 

 

тельную

экономию капитальных

за­

 

 

 

 

 

 

 

трат

(массово

выпускаемое

обору­

 

 

 

 

 

 

 

дование

 

дешевле),

ремонтных

рас­

 

 

 

 

 

 

 

ходов

(взаимозаменяемость

узлов,

 

 

 

 

 

 

 

меньшие

требования к ассортименту,

 

 

 

 

 

 

 

хранящемуся на складах), но и об­

 

 

 

 

 

 

 

легчает

 

собственно

проектирование

 

 

 

 

 

 

 

и расчет

систем регулирования.

 

 

 

 

 

 

 

 

Этими

свойствами обладает

вы­

Рис.

109.

Частотные

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

в

пускаемая

в

ЧССР

агрегатная

 

комплексной

плоскости

 

 

система

 

УРС,

которая

охватывает

 

 

 

 

 

 

 

широкий

 

ассортимент

агрегатных

аналоговых

и

цифровых

узлов в

унифицированном

конструктив­

ном

исполнении, чаще

всего

на

печатных

 

схемах.

Температура

окружающей среды при использовании этой

системы

равна

5—

50° С,

а относительная

влажность

воздуха

45—80%-

Ниже

пере­

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

узлы.

 

 

 

 

 

 

 

 

 

 

Т р а н з и с т о р н ы е

у с и л и т е л и

 

AZ1, AZ2

системы служат

для

усиления малых напряжений

постоянного тока, поступающих

от датчиков, например, термопар. Усилитель AZ2 является

самым

важным узлом, так как его используют

в качестве операционного

усилителя. С помощью обратной связи его передаточную

функцию

можно

преобразовать

в пропорциональную,

 

дифференцирующую,

интегрирующую

или нелинейную. Для

 

этого

используются

комби­

нации

сопротивлений

и конденсаторов

 

в

цепи

обратной

связи.

По

схеме,

представленной

на

рис. ПО,

получают

элемент

П И Д

с передаточной

функцией

 

 

 

 

 

 

 

 

 

 

 

 

 

F(P)=-7-

1

Я . +

/ ? 2 +

2/?2

с.

 

 

 

 

 

 

1

 

 

 

(165)

 

 

 

 

 

 

 

С{гр

 

 

 

где первый член правой части уравнения является пропорциональ­ ной составляющей, второй — дифференцирующей, а третий — инте­ грирующей составляющей передаточной функции.

170

Закорачиванием конденсатора Ci получают элемент PD, a ne- редаточная функция преобразуется в форму

 

 

 

 

 

 

 

 

 

 

 

 

 

(166)

При условии

 

R\=R2

= 0

получим

интегрирующий элемент

с

пере­

даточной функцией

 

 

 

 

 

 

 

 

 

F{p)

1

 

 

 

 

 

 

 

 

 

 

 

(167)

 

С\гр

 

 

 

 

 

 

 

 

 

 

 

 

Усилители AZ3, AZ4 служат в качестве усилителей мощности

для управления

магнитными

усилителями

и т. д. Если

необходима

большая

мощность,

к

ним

под­

 

 

 

 

 

 

соединяется

усилитель

AZ05.

 

 

 

 

 

 

Усилитель

 

AI02

имеет

вы­

 

 

 

 

 

 

сокое входное

сопротивление. Он

 

 

 

 

 

 

служит

для

создания

элементов

 

 

 

 

 

 

с большими

 

постоянными

вре­

 

 

 

 

 

 

мени. Модулятор AD04

преобра­

 

 

 

 

 

 

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

в мо­

 

 

 

 

 

 

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

 

переменный

с ча­

 

 

 

 

 

 

стотой 50 гц.

 

 

 

функциональ­

РиСі и о

Операционный

усилитель AZ2

Для

получения

 

 

системы УРС

 

 

ных зависимостей используют ло­

 

 

 

 

 

 

гарифмический

элемент

AL01,

электромеханический

узел

A M I

(при монтаже

 

интегратора

или программных

устройств)

и

квад­

ратичный элемент ANaOl.

 

 

 

 

 

 

 

Для

гальванической

развязки

сигналов

используют

раздели­

тельный трансформатор AV1; для разъединения обратных

связей—

изоляционный элемент АѴ02а.

 

 

 

 

 

 

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

Г л а в а V

АВТОМАТИЧЕСКИЕ ЭЛЕКТРОННЫЕ ВЫЧИСЛИТЕЛЬНЫЕ МАШИНЫ (ЭВМ)

1. ОБЩИЕ ПОЛОЖЕНИЯ

Современная автоматическая ЭВМ представляет собой кибер­ нетическое устройство, способное осуществлять с большой скоро­ стью математические и логические операции с большим числом

171

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

Со времени установки первой коммерческой

цифровой

ЭВМ

в 1951 г. скорости вычислений возросли у ЭВМ

в 10 000

раз, а

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

Цель

настоящей

главы — дать читателю основную

качествен­

ную информацию о

принципах работы

автоматических

ЭВМ, об

основах

их программирования и использования. Более подробные

данные

приводятся

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

литературе, на

которую

даны ссылки.

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

Способы программирования, использования, расчета и т. д. цифровых и аналоговых ЭВМ различны, поэтому эти вопросы бу­ дут рассмотрены отдельно. О гибридных ЭВМ, применение кото­ рых в настоящее время ограничено, будет сообщено на с. 0000.

Цифровой вычислительный процесс .

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

172

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

Формализованная запись

программы вычисления

у—(а + Ь)2\

а = 7,281, Ь = 16,31 показана ниже:

 

 

Номер строки

 

 

 

 

 

в тетради

 

 

 

Команда

 

 

 

. Возьми число из 7-й строки и запиши его

в регистр

арифмо­

2.

 

метра

 

 

 

 

. Прибавь к показанию

арифмометра число

из 8-й строки

3.

. Показание

арифмометра запиши во вспомогательную строку 6

4.

. Показание

арифмометра умножь на данные строки 6

 

5.

. Показания

арифмометра перепиши в 9-ю

строку

 

6.

'.

7.281

 

 

 

 

 

 

 

 

 

 

.

16,31

 

 

 

 

Работа оператора, проводящего эти вычисления, показана на рис. 111. Участие человека в вычислительном процессе создает воз-

 

 

 

Входное

 

 

 

 

 

 

устройство

 

 

 

 

 

РасчетчиЯ

 

 

 

 

Про-

О

Регистр\

 

 

 

 

\гртяа\

 

 

 

 

уінные]

. л .

 

 

 

 

 

 

 

Упоминающее

Блок

Арифметическое

Нанипулятор

Арифпометр

устройство

управления

устройство

Рис. 111. Схема

ручного цифрового вычис­

Рис. 112. Блок - схема

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

лительного процесса

 

вой

Э В М

 

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

Блок-схема, объясняющая работу цифровой ЭВМ, показана на рис. 112. Оператор в ЭВМ заменен блоком управления (устрой­ ством управления), который управляет вычислительным процессом в соответствии с программой, заложенной в память цифровой ЭВМ.

173

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

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

Аналоговый

вычислительный

процесс

Решим дифференциальное уравнение:

 

s" (0=«is' (0 + floS (0 =

К,

 

которое описывает, например, разгон регулируемого серводвига­

теля. Здесь

s(t)—-угол

поворота

выходного вала в

зависимости

от времени t,

К, ßo, ai — константы.

 

 

Это уравнение можно

решить

аналитически или

в цифровой

форме. Можно его «решить» еще и таким способом: установить оборудование и определять разгон s(t) экспериментально, изменяя положение выходного вала в зависимости от времени. Причем сле­ дящий привод можно установить иного масштаба (например, мень­ ший) и таким образом создать модель собственно оборудования, на котором проводятся измерения. Создавать модель на аналогич­

ной физической основе,

т. е. на базе микродвигателя, сервоусили­

теля и т. д., не обязательно, так как достаточно использовать

обо­

рудование, аналогичное

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

том

смысле, что

изменения

физических величин будут

качественно

(в масштабе)

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

решаемой

проблемы.

 

 

 

 

"Если приведенное выше уравнение будет описывать, например, скорость химической реакции в зависимости от времени, то изме­ нение этой скорости можно определить, составив приведенную электромеханическую модель и измерив положение выходного вала s(t). Затем на основании масштаба M модели можно «прокалиб­ ровать» результаты измерений: единица угла поворота вала в за­ висимости от времени =М- (единица скорости химической реакции).

174

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

шаемой

проблемы.

 

 

 

 

 

Устройство, используемое в качестве

аналога

качественной

стороны) различных физических

задач,

должно

быть

агрегатным

и

обеспечивать

решение целого

ряда

физических

(технических

и

др.)

проблем.

Этим требованиям удовлетворяют

аналоговые

ЭВМ, предназначенные для решения систем дифференциальных уравнений. Аналоговая ЭВМ придает независимой переменной ве­ личине при решении дифференциального уравнения вид аналоговой величины — времени (т), а зависимой переменной физической ве­ личине—вид величины электрического напряжения (U). С этими величинами можно легко проводить математические операции: из­ меняющийся во времени сигнал (напряжение) легко интегрируется, суммируется с другим сигналом (напряжением), умножается на константу и т. д.; результирующие характеристики U(т) можно легко регистрировать при помощи осциллографа, самописца, их можно легко измерять при помощи вольтметра.

Приведем наглядный пример аналогового вычисления по схеме: математическая задача->-аналоговое вычисление-крещение задачи.

входное

выходное

преобразование

преобразование

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

t =

t{x),

 

 

где

t — температура, °С,

 

 

X — длина прутков,

м.

 

Входное преобразование может быть обусловлено соотношением

U =

U(x),

 

 

где

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

электрического напряжения, в,

 

% — время, сек.

1 ° С ~ 1 в, 1 Л І ~ 1

сек.)

(Масштаб, например,

По окончании аналогового вычисления получим на осцилло­

грамме изображение зависимости U=U(x).

Выходное преобразо­

вание является обратным входному преобразованию. Проанализировав кривую (7=(7(т), получим представление об

изменении температуры

t = t(x) материала. Кривая

на

осцилло­

графе изображает решение исходного уравнения t~t(x),

а с

по­

мощью аналоговой ЭВМ

смоделирован физический

процесс

на­

грева.

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

175

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