Способ градуировки пирометра излучения и измерения температуры объекта

Способ градуировки пирометра излучения и измерения температуры объекта

Иллюстрации

Показать все

Реферат

 

Изобретение относится к технике оптической пирометрии, осуществляющей бесконтактное измерение температуры с использованием предварительной градуировки по эталонному источнику излучения, и может быть использовано для построения методики градуировки пирометров излучения и модели расчета температур. Сущность способа заключается в том, что при градуировке зависимость сигнала пирометра от температуры U f(T) снимают только для одного диапазона, для которого затем определяют зависимость эффективной освещенности фбтоприемника от температуры ЕЭф. f(T), далее на основании зависимостей рассчитывают зависимость сигнала пирометра от освещенности его фотоприемника ЕЭф. f(U). На основании ЕЭф. f(U), хранящейся в памяти пирометра для остальных диапазонов , определяется только характеристика ЕЭф. f(T) в виде трех коэффициентов т, n, k для каждого диапазона. 2 з.п.ф-лы, 6 ил. сл с

СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК

ist>s G 01 J 5/00

ГОСУДАРСТВЕННОЕ ПАТЕНТНОЕ

ВЕДОМСТВО СССР (ГОСПАТЕНТ СССР) ОПИСАНИЕ ИЗОБРЕТЕНИЯ

К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ (21) 4834006/25 (22) 04.06.90 (46) 23.12.92, Бюл. М 47 (71) Минский радиотехнический институт (72) Ю.Ч.Гайдукевич, Н.И.Домаренок, А.П.Достанко, В,М,Марченко и И.Г.Мороз (56) 1. Кузмичев В.Н, Калибровочная кривая цифровых анализаторов тепловых полей и методика ее получения//Оптико-механическая промышленность. — 1979 — М 4, — с.46 — 48.

2, Andrometer F., Ardofot und Ardocol—

Strahl Ungspyrometer mit temperatur—

linearem Stromaus — gang, "Kramer Horst, Menge Kurt." Siemens Energletec nn.", 1983, 5,N 5, 246 — 247 (прототип).

3. Коган А,В., Мануйлов Э.А, Определение характеристик пирометров частичного излучения по их аппаратным функциям //

Приборы и системы управления. 1985. — М

2, — c,16-17. (54) СПОСОБ ГРАДУИРОВКИ ПИРОМЕТРА

ИЗЛУЧЕНИЯ И ИЗМЕРЕНИЯ ТЕМПЕРАТУPbI ОБЪЕКТА

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

Известен способ градуировки пирометра излучения и измерения температуры объекта I1j, заключающийся в получении Ж 1783322 А1 (57) Изобретение относится к технике оптической пирометрии, осуществляющей бесконтактное измерение температуры с использованием предварительной градуировки по эталонному источнику излучения, и может быть использовано для построения методики градуировки пирометров излучения и модели расчета температур. Сущность способа заключается в том, что при градуировке зависимость сигнала пирометра от температуры U = f(T) снимают только для одного диапазона, для которого затем определяют зависимость эффективной освещенности фотоприемника от температуры Езф. =

=fQ, далее на основании зависимостей рассчитывают зависимость сигнала пирометра от освещенности его фотоприемника Ез . =

=f(U). Ha основании Еэф. = f(U), хранящейся в памяти пирометра для остальных диапазонов, определяется только характеристика Е ф, = f(T) в виде трех коэффициентов m, и, k для каждого диапазона. 2 з.п.ф-лы, 6 ил. зависимости сигнала пирометра от температуры эталонного излучателя, хранении грэдуировочной характеристики и измерении по ней действительной температуры объекта.

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

Недостатком этого способа является то, что можно производить контроль температур только тех объектов, коэффициент черноты которых соответствует коэффициенту черноты эталонного излучателя. Для остальных же обьектов температура, определяемая по градуировочной характеристике, называемая условной температурой Т(, будет отличаться от истинной температуры Т на величину ЛТ, определяемую различием коэффициентов излучения я объекта и эталонного излучателя, Для расширения области применения пирометра необходимо расширять диапазон контролируемых им температур. При контроле широкого диапазона температур из-за ограниченного динамического диапазона сигнала пирометра весь температурный дйапазон необходимо разбивать на ряд диапазонов, Это может производиться путем дйскретного. неселективного ослабле- ния излучения или установкой селективных фильтров, причем второй способ в некаторйх случаях является предпочтительней, т.к. измерения в различййх "спеКтральных диапазойах позволяют обеспечить максимальное температурное разрешение для данного температурного диапазона, а также выбирать для измерений onðåäåéåííûå окна прозрачности промежуточной среды.

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

Наиболее близким к заявляемому является способ градуировки пираметра излучения и измерения температуры объекта (21. заключающийся в получейии зависимостей сигнала пирометра от температурй эталонного излучателя U=f{T ) и определении зависимости истинной температуры от условной Т = f{T ) для всех диапазонов измерений, отличающихся спектральными характеристиками, хранении данных зависимастей в памяти пирометра и использовании их для расчета истинных температур.

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

5 = 1{Т1), где Т) — температура эталонного излучателя с известной излучательнай способностью я. После этого определяется зависимость измеряемой температуры от коэффициента черноты объекта F., напри10 мер в виде Т=Т.1 (1+ д Т), где DT = f (c).

Снятие данной зависимости основано на том положении, что влияние коэффициента черноты объекта F.,на измеряемую температуру эквйвалентно влиянию коэффициента

15 пропускания т. Коэффициент пропускания

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

20 объектива Н будет эквивалентно изменению е, . Математически это будет показано ниже. Такал зависимость Т от е, используется при измерении истинной температуры в виде Т = f(Ty ), где Т вЂ” условная темпера25 тура, которая берется непосредственно из градуировачной характеристики U = f{T ) .

{фиг.2), Снятие градуировачнай характеристики

U = f(T) ) и определение зависимости Т =

=f(T ) производится в каждом температурном диапазоне, отличающемся спектральной характеристикой фильтра. При большом количестве диапазонов процесс градуировки становится достаточно трудоемким, Кроме этого, для проведения измерений во всем температурном диапазоне необходимо хранить характеристики Т = 1(Т ) для каждого спектрального диапазона, что при табличном описании данной характеристи40 ки требует наличия значительного объема памяти для их хранения; аналитическое же описание данных характеристик с достаточной точностью не всегда эффективно, т.к. зависимость U = f{T) ) определяется совокупностью таких характеристик, как зависимость эффективной освещенности фотоприемника пирометра ат температуры .объекта Езф = f(T), описываемой законом

Планка, и характеристикой свет-сигнал пираметра U=f(E q), которая включает в себя световую характеристику приемника излученйя, имеющую для некоторых датчиков форму, трудна поддающуюся описанию, и передаточную характеристику всего элект55 раннаго тракта пирометра, Целью изобретения является ускорение процесса градуировки.

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

1783322 ры, в первом диапазоне определяют зависимость эффективйой освещенности от температуры E y=f(T), и ка основании ее и зависимости О=f(T) рассчитывают зависимость сигнала пирометра от освещенности 10

Т—

20

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

U=f(E q), (1)

5 где U — выходной сигнал;

Е ф — эффективная освещенность, создаваемая на приемной площадке пирометра;

f(E-„.y) — функция преобразования свет

-сигнал, включающая в себя световую характеристику приемника излучения и передаk — 1л — Е

50 (2) температуры объекта, заключающемся в получении зависимостей сигнала пирометра от температуры эталонного излучателя

О=f(T ), запоминании полученных зависимостей и определении истйкной температуt„ro фотоприемника U=f(E+y), а для остальных диапазонов определяют только зависи- . мость E>y— = f(T), используя при этом едйную для всех диапазонов характеристйку

U=f(E>y), на основании которых определяют истинную температуру обьекта.

В данном способе предлагается вместо зависимостей T=f(T ) в каждом спектральном диапазоне определять зависимость

Е,ф=-f(T), а зависимость U=f(Ty ) определять только для одного диапазона, по которой рассчитывают зависимость U=f(E>y), используя ее при градуировке остальных диайазонов и измерениях истинной температуры (фиг.3). На основании этого можно сделать вывод, что заявляемый способ удовлетворяет критерию "новизна", Использование в качестве градуировочной зависймости U=f(E>y) вместо U=f(T) позволяет применять ее во всех спектральных диапазоках, т.к. она содержит в себе только энергетические характеристики пирометра: световую характеристику приемника- йзлучения, передаточную характеристику тракта усиления и не зависит от спектрального диапазона, в котором производйтся измерение, поэтому отпадает необходимость снятия и хранения множества градуировочных характеристик. Зависимость же температуры от эффективной освещенности и коэффициента черноты объекта поддается достаточно точно аналогическому описанию, т,к. она основана на законе Планка.

При определении температуры йо определенному Еэф с помощью характеристикИ

Еф=Щ используют формулу где n, m, k постоянные для данного диапа.зона коэффицменты, рассчитанные в процессе градуировки пирометра по трем точкам зависимости E,y=f(T); г — коэффициент черйоты объекта.

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

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

Еэф Но (1 )

k — In

Н (1 — —.)2 те

f где Н, F, а, г, Но, Г>, ао, т>- соответственно светосила объектива, фокусное расстояние объектива, расстояние до объекта, коэффициент пропускания промежуточной среды во ввремя измерения темперагурыг"и в момент"градуировки.

На фиг.1 представлены градуировочные характеристики пирометров; ка фиг.2 — ал-, горйтм градуировки" с за пом"ин а кием

U=f(Tq) и д T=f(e ); на фиг.3 — алгоритм грэдуировкй пйрометра и измерения темпе- ратуры по предлагаемому способу; на фиг.4 — типовая относителькая спектральная характеристика пирометра; на фиг.5— три градуировочные характерйстики пираметра для одного фильтра при различных значениях светосилы-объектива Н, на фиг.б — градуировочная характеристика и три точки; используемые-фЫ определения зависимостии E>y=f(T) для тем пе ратурн ых диапазонов, кроме первого точную характеристику тракта усиления.

Эффективная освещенность на фотоприемнике пирометра описывается формулой(31 = — Н (1 — — )2тя f ã(ÀT) S {А) а2, 1 а

1783322 где Н, F, à, r е — соответственно светосила объектива, фокусное расстояние объектива, расстояние до объектива, коэффициент пропускания промежуточной среды, коэффициент излучения объекта. В общем случае t и 5

F, я вBл я юeтTtс:я ф у нHкKц и я мM1и 4 о0т Т д л1и4нHbы 1 вBоoIл1 нHы, но в пирометрических расчетах их обычно заменяют интегральными эквивалентами; г(Л, Т) — спектральная энергетическая светимость абсолютно черного тела, описы- 10 ваемая законом Планка, но для Х Т<3000

К мкм данная зависимость с достаточной точностью описывается законом Вина

С2 г(1„Т) С1 1, Ч т (9) р) 15 (10) где С1 и С2 — постоянные Планка;

А- длина волны излучения;

Т вЂ” температура тела;

3(Л) — относительная спектральная ха- 2О рактеристика системы, характеризующая спектральные избирательные свойства всего оптического тракта; объектив, фильтры, приемник излучения. Экспериментально установлено, что относительная спектральная 25 характеристика приемников излучения, в частности видиконов, достаточно точно описывается функцией вида (3)

К2 (5)

К2 К1 knax (6) — + k

45 Еэф=е T+m (12) К2

Ко Апах е виэх, одним из достоинств которой является то, что ее коэффициенты m, n u k могут быть получены практическим путем из градуировочных характеристик по трем точкам, Достоверность выражения (12) можно доказать путем сравнения выражений (11) и (12). Приравнивая правые части уравнений (11) и (12), а также выражений их первых и вторых производных и решив систему из

55 трехуравнений, находим зависимости коэффициентов п1 и Й ОТ К1 и К2

1 (1 3)

2 — +—

С2 Т

Я(А)=КО А е (. (4)

Значения К>, К1 и К2 функционально связаны со значением S1 относительной спектральной характеристики (фиг.4) по следующим формулам (математический вывод данных соотношений дан в приложении

1): где ilmax — длина волны максимальной спектральной чувствительности;

Л1 — длина волны на длинноволновой ветви спектральной характеристики, на которой спектральная чувствительность равна

S1;

С учетом этого формула (2) запишется в виде

1 F э — (К1+В)

Š٠— H (1 — — )2 Т е С1 Ko f À х

4 а о хе Т бА д(— + K2) (8)

После раскрытия интеграла, используя преобразования Лапласа, выражение принимает вид

Ьа= — Н (1 - — )2 t e C1 Ko (— +

1 Р С2

4 а Т

{ К1 +4)

К2) Г(К1 +4). где Г(К1+4) — гамма функция, которая для целых К1 может быть записана формулой (К1+3)!. Обозначив

Кз- — Н(1- — )2гес1 К., F

4 а эффективную освещенность можно выразить формулой — (K1+4) эф Кз (— Ò +К2 ) Г (К1 + 4 ). (11)

Данная формула точно описывает зависимость эффективной освещенности от температуры, но значения К1 и К2 должны быть заранее определены по спектральным характеристикам приемника излучения, оптики и промежуточных фильтров. В большинстве случаев не имеется точных данных обо всех спектральных характеристиках, а решение системы из трех уравнений (11) для трех температур, при»x практическом определении из градуировочных характеристик с целью нахождения К1, К2 и Кз, приводит к трансцедентному уравнению.

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

1783322

n = 2m(1+ — )(К1+1), Т (14)

Коэффициент k является масштабным и определяется отдельно, Коэффициенты и! и и зависят от температуры — это означает. что выражение (12) аппроксимирует выражение (11) с достаточной точностью в некотором температурном диапазоне, По практическим расчетам точ. ность соответствия формулы (12). формуле

{11) достаточно высокая.

Коэффициенты m, n uk определяются из градуировочных характеристик по tp6M точ-кам.

Значение эффективной освещенности

Еэф, как видно из формулы (9), пропорционально светосиле объектива Н, йоэтому зависимость эффективной освещенности от температуры можно определить из зависимости температуры от светосилы объектива на одном уровне выходного сигнала пирометра.

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

1, Путем полного открытия диафрагмы светосила объектива Н устанавливается:на максимальное значение, равное Н1. Затем постепенно увеличивая температуру эталонного излучателя, точку за точкой, снимается зависимость выходного сигнала U от температуры Т во всем динамическом диапазоне пирометра, Из данной характеристики берется точка на уровне сигнала Up, для которой определяется Т1. После этого диафрагма прикрывается до значения, соогветствующего Нг, и температура эталонного излучателя увеличивается до тех пор; пока выходной сигнал не станет равным Up. Данное значение температуры обозначается Тг.

Температура Тз получается аналогичным образом для светосилы обьектива Нз. На фиг.5 представлены три зависимости U=f(T), полученные при разных значениях светосилы объектива H. Снятие двух полных кривых для Нг и Н3 необязательно. Достаточно определить температуры Тг и Тз при одинаковых уровнях выходного сигнала пирометраUp, соответствующего температурной точке

Т! на градуировке Н1. Необходимость проведения измерений на одном уровне сигнала Up, для трех градуировок, а не трех 0 на одной градуировке, вызвана тем, что для пирометров функция преобразования светсигнал Ре(Еэф) в формуле (1) является в большинстве случаев нелинейной и, вообще не известной. Для исключения влияния Этой нелинейности на результат градуировкй измерения Т1, Тг, Т3 необходимо проиэводить на поСтоянном уровне 4.

2. Имея три температуры Т1. Тг, Тз и три значения светосилы объектива Н1, Нг. Н3. с учетом выражений (12) вычисляют значения

m, n u k no следующим формулам:

С} Тз 12

m = (15) (16) (Т! — Тг )

n— и с .. (—

T1+m где

5 . 1.12 T1 — Т3

Ц

L13(T1 — T2) (17) (18) L12 = In

Еэ Т1 Н1

= In

Е,ф (T2)

1-13 = In

Еэ Т1 H1

=In

Еэф (T3) НЗ (19) (20) Эффективная освещенность точки

Еэф{Т1) принимается равной 1 для расчета k из условия нормировки.

Таким образом, в результате этих вы= числений определяется градуировочная зависимость Еэф=f(T) для первого спект30 рального диапазона описывается выражением (12) с учетом данных коэффициентов

m, n и К которые заносятся в память пирометрической системы.

3. Для определения зависимости выходного сигнала от температуры необходимо определить функцию преобразования светсигнал в формуле (1). Данная функция определяется из снятой по п.1 характеристики

0=1(Т) путем пересчета значений температуры Т на этой кривой в эффективность освещенность по формуле (12) с учетом рассчитанных по п.2 коэффициентов m, n u

k. Полученная таким образом градуировоч45 ная зависимость U=f(Esp) записывается в память пирометра.

4. Зависимость U=f(E>y) является постоянной для данного пирометра, Поэтому при градуировке пирометра в остальных спектральных диапазонах нет необходимости повторного снятия денной характеристики.

Расчет значений n, m u k для других спектральных диапазонов выполняется также по формулам (15-17), но при этом расчет производится по точкам Т1, Т2 и Т3, взятым на одной температурной кривой U=f(T) (фиг.6), а не на одном уровне трех характеристик Н1, Нг и Нз, как в первом диапазоне.

Это позволяет ускорить и упростить про1783322

15 и п Еи (24) 20

40 объекта измерения

Еэф с (21) 45

50 цесс градуировки во всех последующих спектральных диапазонах.

Для этого в оптическую систему пирометра устанавливается соответствующий данному диапазону светофильтр и повторяется процесс снятия зависимости U=f(T), аналогично описанному в п.1, Для ускорения процесса градуировки достаточно определить только три точки этой характеристики

U(T>), U(T2) и 0(Тз). По градуировочной зависимости U=f(E>y) для этих точек определяются соответствующие Е ф(Т ); Еэф(Т ) и

1 ф(Тз), которые и используются для перерасчета спектральных коэффициентов m, и и k по формулам (15) — (20). Рассчитанные коэффициенты m, n u k записываются в память и фактически описывают соответствующую градуировочную зависимость E>y-=-f(T) с помощью выражения (12).

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

0=1(Еэф) и трех числовых коэффициентов m, n и k для всех спектральных диапазонов, необходимых для описания соответствующих градуировочных зависимостей Езф=f(tj.

Реализация описываемого способа в процессе йзмеренйя температуры-при onределении истинной температуры по полученным градуировкам производится в следующей последовательности: а) по полученному значению выходного сигнала U, через хранящуюся в памяти градуировочную зависимость 0=1(Еэф), определяют значение эффективной освещенности Еэф, б) полученное значение осввщенности корректируется по коэффициенту черноты или, если характеристики обьектива, расстояние до объекта, коэффициент пропускания промежуточной среды в процессе измерений отличается от соответствующих параметров в процессе градуировки, то Е необходимо корректировать по формуле где г„, т, — определяется по коэффициентам в. момент градуировки и в момент измерения температуры по формуле г= Н(1 — — Яrя (23)

F а в) по скорректированному значению Еи определяется температура с учетом хранящихся в памяти расчитанных в процессе градуировки в конкретном спектральном диапазоне коэффициентов m, n u k

Заявляемый способ может быть реализован в любых пирометрах, осуществляющих градуировку и измерение температуры под управлением ЭВМ, например, в тепловиэионном пирометре ИИ-42Т (5).

Таким образом, из вышерассмотренного можно заключить, что использование в качестве градуировочной зависимости

0=1(Еэф) вместо U=f(T) позволяет использовать ее во всех спектральных диапазонах, тем самым сокращая количество хранимых полных градуировочных характеристик до одной. Кроме этого, упрощается процесс определения зависимости измеряемой температуры от коэффициента черноты для остальных диапазонов, отличающихся спектральными характеристиками фильтров, вместо снятия трех градуировок снимаются только три точки на одной градуировке.

Формула, используемая для расчета температуры по определенному значению

Еэф, имея довольно простое написание, обладает высокой точностью, Так, по практическим расчетам при использовании в качестве датчика в тепловизионной лирометрической системе видикона ЛИ-421 с максимумом спектральной чувствительности на Лбах=580 нм и с полушириной

Ail, =180 нм (КО=70, К2=-40,5), максимальная погрешность описания эффективной осве. щенности формулой (12) в диапазоне температур 800...2000 К составляет 0,43 %.

Данная погрешность уменьшается при увеличении Т, а также коэффициента К2, который увеличивается при уменьшении ширины спектральной характеристики. Для широко распространенных узкополосных гирометров (Я х = 650 HM и ЛЛ =10 нм) максимальная погрешность для того же температурного диапазона равна 0,000014%.

Расширенная формула позволяет при измерении температуры также учесть характеристики объектива, расстояние до объекта," коэффициент пропускания проме>куточной среды.

Формула изобретения

1, Способ градуировки пирометра излучения и измерения температуры объекта, заключающийся в получении зависимостей сигнала пирометра от температуры эталон1783322 ного излучателя U=f(I), запоминании полученных зависимостей и определении истинной температуры, отличающийся тем, что, с целью ускорения процесса градуировки, в первом диапазоне определяют зависимость эффективной освещенности от температуры Е ф=1(Т) и на основании ее и зависимости U-f(T) рассчитывают зависимость сигнала пирометра от освещенности его фотоприемника U=f(E,ô), а для остальных диапазонов определяют только зависимость E>y=f(T), используя при этом единую для всех диапазонов характеристику

U=f(E>y), на основании которых определяют истин ную температуру обьекта.

2. Способ по п.1, отличающийся тем, что расчет истинной температуры Т объекта при использовании пирометра ча-. стичного излучения производят по формуле.

Т =п

" fTI

k — In — Е3

Е где п, m, k — постоянные для данного диапазона коэффициенты, рассчитанные в процессе градуировки пирометра по трем точкам зависимости Е y=f(T); я — коэффициент черноты объекта.

5 З.Способ попп.1,отличающийся тем, что, с целью учета характеристик объектива, расстояния до объекта, коэффициента пропускания промежуточной среды, отличных от значений соответствующих парамет10 ров, в момент градуировки расчет истинной температуры производят по формуле

Т

Еэф Но (1 — — ) то

Io ао

k — ln

Н (1 — — ) те

f а где Н, F, а, r, Н,, Fp, ао, г-соответственно

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

1783322

1783322 т т юг 8

Составитель Ю. Гайдукевич

Техред М.Моргентал Корректор Н. Ревская

Редактор

Производственно-издательский комбинат "Патент", г. Ужгород, ул.Гагарина, 101

Заказ 4507 Тираж . Подписное

ВНИИПИ Государственного комитета по изобретениям.и открытиям при ГКНТ СССР

113035, Москва, Ж-35, Раушская наб., 4/5