Авторефераты по всем темам >>
Авторефераты по техническим специальностям
На правах рукописи
Гринюк
Александр Владимирович
термоэлектрический метод контроля теплофизических свойств твЁрдых материалов
05.11.13 Приборы и методы контроля природной
среды, веществ, материалов и изделий
автореферат
диссертации
на соискание учёной степени
кандидата технических наук
Москва 2012
Работа выполнена на кафедре УАвтоматизация производственных процессовФ Новомосковского института Российского химико-технологического университета им. Д.И.Менделеева
Научный руководитель | доктор технических наук, профессор Беляев Юрий Иванович |
Официальные оппоненты: | доктор технических наук, профессор Рылов Владимир Аркадьевич главный конструктор проекта НПО Химавтоматика кандидат технических наук, доцент опатин Александр Геннадьевич Университет Российской академии образования Новомосковский филиал (НФ УРАО) |
| Ведущая организация: | ООО Научно Ц исследовательский проектно Ц изыскательский межрегиональный институт азотной промышленности, г. Москва |
Защита состоится 26 апреля 2012 года в 16-00 часов на заседании Диссертационного Совета Д 212.145.02 в Московском государственном университете инженерной экологии по адресу: 105066, г. Москва, ул. Старая Басманная, 21/4.
Отзывы в двух экземплярах, скреплённых гербовой печатью, просим направлять по адресу: 105066, г.Москва, ул.Старая Басманная, 21/4.
С диссертацией можно ознакомиться в библиотеке университета.
Автореферат разослан 24 марта 2012 г. 28 августа 2006
Ученый секретарь
диссертационного Совета
Д 212.145.02 ___________ д.т.н., доцент Мокрова Н.В.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы.
Одним из основных направлений современной науки и техники является разработка прогрессивных технологий производства новых материалов с заданными теплофизическими свойствами. Причём, потребность в таких материалах с каждым годом увеличивается.
Существующие установки измерения теплофизических свойств имеют различные недостатки, например такие как:
- Для приборов основанных на стационарном методе измерения характерными недостатками являются длительное время измерения и громоздкость теплоизмерительной системы;
- Существующие нестационарные методы хотя и обладают высоким быстродействием, но также имеют ряд недостатков:
- в методе азерной вспышки предъявляются повышенные требования к условиям эксперимента, и он является косвенным методом;
- метод монотонного нагрева нуждается в длительной градуировке;
- практическое осуществление метода мгновенного источника тепла, сопряжено с существенными погрешностями, возникающими из-за сложности точного определения значения экстремума температуры Tmax и времени max.
Поэтому данная работа направлена на устранение недостатков существующих методов путём применения в теплофизических измерениях системы элемент Пельтье - адиабатическая камера.
Цель работы. Целью диссертационной работы является решение актуальной научно-технической задачи создания термоэлектрического метода контроля теплофизических свойств твёрдых материалов, и разработка на его основе приборов контроля теплофизических свойств, а также повышение эффективности теплофизических измерений.
Для достижения поставленной цели необходимо решить следующие задачи:
- получение математической модели термоэлектрических процессов в системе элемент Пельтье - образец теплоприёмник;
- получение математической модели термоэлектрических процессов в системе элемент Пельтье - образец адиабатическая камера;
- разработка метода измерения теплофизических свойств материалов на основе математической модели системы элемент Пельтье образец - теплоприёмник и элемент Пельтье образец - адиабатическая камера;
доказать преимущества применения элемента Пельтье в теплофизических измерениях;
- экспериментальное обоснование эффективности разработанного метода контроля за счёт применения системы элемент Пельтье - образец адиабатическая камера в измерениях теплопроводности и теплоёмкости.
Методы исследования.
В работе использован метод математического моделирования, а также экспериментальное исследование особенностей протекания теплофизического эксперимента в твёрдых материалах.
Научная новизна:
- получена математическая модель стационарных термоэлектрических процессов в системе элемент Пельтье образец;
предложено применение адиабатической камеры в системе элемент Пельтье - образец для измерения теплофизических свойств материалов;
получена математическая модель стационарных термоэлектрических процессов в системе элемент Пельтье образец - адиабатическая камера;
- разработан метод измерения теплофизических свойств материалов в системах элемент Пельтье - образец - теплоприёмник и элемент Пельтье образец - адиабатическая камера;
- впервые предложено применение элемента Пельтье для измерения теплового сопротивления и теплопроводности материалов;
- экспериментально исследовано применение адиабатической камеры в для измерения теплофизических свойств материалов.
Практическая ценность:
полученные результаты использованы для разработки систем измерения теплофизических параметров;
применение адиабатической камеры в системе элемент Пельтье образец позволяет повысить точность и сократить время измерений;
разработанный метод можно использовать для разработки средств оперативного контроля качества промышленных изделий;
исследованы преимущества применения элемента Пельтье в теплофизических измерениях;
исследованы преимущества применения адиабатической камеры в теплофизических измерениях.
Реализация научно-технических результатов. Основные идеи и результаты теоретических и экспериментальных исследований реализованы в измерителе теплопроводности КИТ-02 - и в измерителе теплоёмкости теплоизоляционных и композиционных материалов БИТ-01С, внедрённом в ОАО Корпорация Московский институт теплотехники и ОАО Центральный научноЦисследовательский институт специального машиностроения.
Материалы диссертационной работы используют в учебных курсах и научно-исследовательской практике Новомосковского института РХТУ.
Достоверность работы экспериментально подтверждается совпадением результатами математического моделирования с законами Пельтье и Зеебека.
Разработанный метод контроля теплофизических свойств твердых материалов апробирован на эталонных образцах. Результатом явилось совпадение показаний прибора с сертифицированным значением эталонов для измерителя теплопроводности КИТ-02 - с погрешностью, равной 5 %, а для измерителя теплоёмкости БИТ - 01С 3 %.
Автор выносит на защиту:
термоэлектрический метод контроля теплофизических свойств твёрдых материалов;
математическую модель термоэлектрических процессов в системах элемент Пельтье - образец и элемент Пельтье - образец - адиабатическая камера;
применение системы элемента Пельтье - образец адиабатическая камера в измерениях теплопроводности и теплоёмкости;
- экспериментальное исследование системы элемент Пельтье - образец - адиабатическая камера;
- приборы, созданные на основе термоэлектрического метода контроля теплофизических свойств твёрдых материалов.
Апробация и публикации. По теме диссертационной работы опубликовано десять печатных работ, в том числе две в журналах, рекомендованных ВАК.
Результаты работы докладывались и обсуждались на конференциях:
Математические методы в технике и технологиях (Ярославль 2007, Пенза 2011), XXI Международной конференции молодых ученых по химии и химической технологии (Москва, 2007 г.), РХТУ им. Д.И.Менделеева, Новомосковский институт, 2007 - 2011 гг.
Структура и объем работы. Диссертация состоит из введения, четырёх глав, заключения и списка итературы с 156 наименованиями.
Изложена на 125 страницах машинописного текста. Включает 63 рисунка, 3 таблицы.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность и значимость работы, особенности проблем современных методов измерения теплофизических параметров твёрдых тел, вытекающие из этих особенностей цели исследования, новизна и задачи, решённые в данной работе.
В первой главе рассматриваются основные способы применения термоэлектрических элементов в качестве термостабилизаторов и измерителей теплового потока на примере серийно выпускаемых теплофизических приборов, также рассмотрены основные методы и средства измерения теплопроводности твёрдых материалов.
Во второй главе предложен термоэлектрический метод контроля теплофизических свойств твёрдых материалов, основанный на термоэлектрических эффектах Пельтье и Зеебека, возникающих в термоэлектрических модулях - элементах Пельтье.
При протекании электрического тока через элемент Пельтье в нём возникает два тепловых потока: Пельтье и Джоуля, также возникает ЭДС Зеебека, величина, которой пропорциональна перепаду температур на термоэлектрическом элементе. Эквивалентная схема элемента Пельтье показана на рис. 1 а. Взаимосвязь электрических и тепловых процессов отображают уравнения:
, (1)
, (2)
где UП - падение напряжения на элементе Пельтье, I - электрический ток, проходящий через элемент Пельтье, RЭП - электрическое сопротивление элемента Пельтье, EЗб - ЭДС Зеебека, КЗб - коэффициент Зеебека, TП - перепад температур на элементе Пельтье
|
| |
а | б |
Рис. 1. Эквивалентная схема элемента Пельтье (а) и тепловые потоки
в системе элемент Пельтье - образец - теплоприёмник (б):
UR - падение напряжения на активном сопротивлении элемента Пельтье; RtП - тепловое сопротивление элемента Пельтье; 1 - образец; 2 - элемент Пельтье; 3 - теплоприёмник; t1 - температура нижней обкладки (поверхности) элемента Пельтье; t2 - температура верхней обкладки (поверхности) элемента Пельтье и нижней стороны образца; t3 - температура верхней стороны образца; qRx, qвх , qвых - тепловые потоки в системе элемент Пельтье - образец - теплоприёмник
Эти особенности элемента Пельтье позволяют применять его для измерения теплового сопротивления.
Тепловое сопротивление на участке тепловой цепи вычисляется по формуле:
, (3)
где TRx - перепад температур на образце; qRx - тепловой поток, проходящий через образец.
На рис. 1 б показана система элемент Пельтье - образец - теплоприёмник. Тепловой поток qвых этой системы определяется выражением:
, (4)
где п - коэффициент Пельтье.
Тепловой поток qвх равен выражению:
, (5)
Эквивалентная схема системы элемент Пельтье образец приведена на рис. 2.
Рис. 2. Эквивалентная схема системы элемент Пельтье - образец Ц теплоприёмник: UП - падение напряжения на элементе Пельтье; Rб - балластное сопротивление; RtХ - тепловое сопротивление образца; RtП - тепловое сопротивление элемента Пельтье; qвх , qвых - тепловые потоки, генерируемый элементом Пельтье; TП - перепад температур на элементе Пельтье; TRx - перепад температур на образце |
Из рис. 2 следует, что qRx = qвх, тогда тепловой поток, проходящий через образец RtХ определяется выражением (5). Так как температуры t1 и t3 (рис. 1 б) зафиксированы теплоприёмниками, то t1 и t3 можно считать неизменными и равными температуре окружающей среды, т. е. t1 = t3 = const. |
Примем, что распределение температур в элементе Пельтье инейно и Джоулево тепло распространяется одинаково в обе стороны, тогда распределение температур в системе элемент Пельтье - образец теплоприёмник будет соответствовать показанному на рис. 3. |
| Рис. 3. Распределение температур в системе Пельтье - образец - теплоприёмник |
При достижении стационарного режима и при условии равенства температур t1 = t3 (рис. 3) тепловое сопротивление RtХ определяется соотношением:
. (6)
Так как TП = - TRx, а qRx = qвх то получим следующее выражение:
. (7)
Подставив соотношения (1) и (2) в выражение (7), получим:
. (8)
В связи с тем, что теплоприёмники 3 на рис 1. б имеют конечную теплоемкость, то при возникновении тепловых потоков в системе неизбежно возникает разность температур t1 и t3, что является источником методической погрешности. Для обеспечения равенства t1 = t3 необходимо применять специальные устройства. Эти недостатки можно устранить, применив при измерении адиабатическую камеру (рис. 4)
Рис. 4. Тепловые потоки в системе элемент Пельтье - образец -
адиабатическая камера:
1 - образец; 2 - элемент Пельтье; 3 - адиабатическая камера; 4 - тепловой
поток, протекающий в камере
Эквивалентная схема системы элемент Пельтье - образец - адиабатическая камера приведена рис. 5.
Рис. 5. Эквивалентная схема системы элемент Пельтье - образец - адиабатическая камера:
Rtа - тепловое сопротивление адиабатической камеры; qП - тепловой поток Пельтье, протекающий в системе; qД - тепловой поток Джоуля, протекающий в системе; TRа - перепад температур в адиабатической камере
Примем, что тепловой поток в адиабатической камере представляет собой суперпозицию теплового потока Пельтье и теплового потока Джоуля. Джоулево тепло, выделяющееся в элементе Пельтье, распространяется через две поверхности (обкладки) элемента в противоположные стороны, а выделяющийся тепловой поток Пельтье формируется в одном направлении, в зависимости от полярности приложенного напряжения. Выделившиеся тепло Джоуля идёт на разогрев адиабатической камеры, а тепловой поток qвыделенное, создаваемый эффектом Пельтье, протекающий по тепловой цепи Rtа - RtХ - RtП, определяется выражением:
. (9)
Распределение температур в стационарном режиме в системе тепловых сопротивлений Rtа - RtХ - RtП и тепловом потоке qвыделенное = const будет соответствовать показанному на рис. 6.
| Рис. 6. Распределение температур в системе Rtа - RtХ - RtП |
При тепловом равновесии системы элемент Пельтье - образец - адиабатическая камера должны выполнятся условия неразрывности теплового потока:
, (10)
где qП - тепловой поток в элементе Пельтье, qa - тепловой поток в адиабатической камере, qX - тепловой поток в образце.
Преобразовав условия (10) получим:
, (11)
. (12)
Из уравнения (12) и выразим TП получим:
. (13)
При условии << получим:
. (14)
Из выражения (4) следует:
. (15)
Подставим (15) в выражение (14) и получим:
. (16)
Подставим (1) и (2) выражения в (16) и получим расчетную формулу теплового сопротивления RtX:
. (17)
Из выражения (17) выразим UП и получим статическую характеристику системы элемент Пельтье - образец - адиабатическая камера в следующем виде:
(18)
Применение адиабатической камеры позволяет значительно упростить теплоизмерительную систему, так как отпадает необходимость в измерении температуры t3.
Таким образом следует, что термоэлектрический метод контроля теплофизических свойств твёрдых материалов заключается в том, что элемент Пельтье приводится в соприкосновение с исследуемым материалом, затем, на него подается ток и по изменению напряжения и тока элемента Пельтье определяют тепловое сопротивление материала:
Далее в главе рассматривается измерительная система с расположением элементов Пельтье с двух сторон образца.
Для определения удельной теплопроводности можно использовать следующую формулу:
, (19)
где H - толщина образца; RtX - тепловое сопротивление; S - площадь поперечного сечения образца.
В третьей главе проводится экспериментальное исследование термоэлектрического метода контроля теплофизических параметров.
Для проверки возможности косвенного измерения теплового потока, создаваемого элементом Пельтье, путём измерения протекающего через него электрического тока, был проведён опыт, суть которого заключалась в измерении теплового потока, генерируемого элементом Пельтье при различных значениях электрического тока (рис. 7).
Полученные зависимости сигнала тепломера от электрического тока для термоэлемента типа ТВ-31-1,0-2,5 при различных направлениях электрического тока показаны на рис. 8.
Рис. 7. Экспериментальная установка для исследования
термоэлектрического метода:
1 - исследуемый элемент Пельтье; 2Ц датчик теплового потока; 3 - теплоприёмник; А1 - амперметр; V1, V2 - вольтметр
Кривая 2, аппроксимированная уравнением q = 0,28I + 0,39I2, соответствует протеканию тока от минусового вывода к плюсовому, а кривая 4, - обратному направлению тока.
Средняя относительная квадратичная ошибка аппроксимации на интервале значений электрического тока от 0 А до 1,4 А составила = 2,9 % для кривой 1 и = 1,4 % на интервале от 0 А до 1,4 А для кривой 3.
Рис. 8. График зависимости сигнала тепломера от тока через элемент Пельтье: 1 Ц выделение тепла на горячей стороне элемента Пельтье; 2 - аппроксимация уравнением q = 0,28I + 0,39I2; 3 - поглощение тепла на холодной стороне элемента Пельтье; 4 - аппроксимация уравнением q = 0,28I - 0,39I2 |
Расхождения экспериментальной и аппроксимированной кривой, можно объяснить принятым допущением о инейном распределении температуры внутри элемента Пельтье.
Для определения режима измерения, при котором обеспечивается максимальная чувствительность к тепловому сопротивлению образцов, был проведён эксперимент, суть которого заключается в измерении теплового потока в системе элементе Пельтье образец при разных значениях тока, протекающего через элемент и на разных образцах. Схема включения термоэлектрического элемента представлена на рис. 7, а результаты измерения - на рис. 9 а. Измерения производились на элементе Пельтье типа ТВ-31-1,0-2,5 и на образцах диаметром 15 и высотой 5 мм.
Из рис. 9 а следует, что максимальная чувствительность находится в диапазоне значений электрического тока 0,5 - 0,6 А. На рис. 10 приведены графики переходных процессов на элементе при ступенчатом изменении тока от 0 до 0,5 А на различных образцах.
Из полученных кривых (рис. 9 б) можно определить постоянную времени элемента Пельтье типа ТВ-31-1,0-2,5, которая составляет примерно 4 с, что позволяет создавать быстродействующие приборы измерения теплофизических свойств.
а | б |
Рис. 9. Сигнал тепломера на различных образцах при различных значения электрического тока (а) и кривые переходного процесса элемента Пельтье при ступенчатом изменении тока от 0 до 0,5 А на различных образцах (б):
1 медь ( = 400 Вт/(мК); 2 - алюминий ( = 220 Вт/(мК)
3 - оргстекло ( = 0.2 Вт/(мК); 4 Цпенополистирол (=0.033 Вт/(мК), 5 - воздух (=0,02 Вт/(мК)
Высокое быстродействие системы элемент Пельтье - образец - адиабатическая камера достигается за счет того, что элемент Пельтье позволяет создавать направленный тепловой поток в адиабатической камере, который в определенных режимах (qП > qД) замыкается в адиабатической камере. В этом случае время установления стационарного режима определяется теплоемкостью адиабатической камеры, элемента Пельтье и образца. Так как теплоемкость адиабатической камеры больше теплоемкости образца, то время установления стационарного режима мало зависит от теплоемкости образца и составляет 30 - 40 секунд (рис. 10). В существующих системах время установления стационарного режима определяется установлением тепловых потоков между теплоприёмниками и окружающей средой.
Таким образом, на основании проведённых исследований выявлены следующие преимущества применения элементов Пельтье в теплофизических измерениях:
- Элемент Пельтье позволяет одновременно генерировать одномерный тепловой поток заданного направления, определять величину этого теплового потока и регистрировать температуру поверхности
- Высокая чувствительность и быстродействие позволяют создавать на основе элементов Пельтье точные быстродействующие (экспрессные) методы измерения теплофизических свойств твёрдых материалов различных геометрических размеров и тонких плёнок.
- Даёт возможность проводить теплофизические эксперименты в широком диапазоне температур от - 120 до +200 С.
- Применение адиабатической камеры совместно с элементом Пельтье привело к существенному уменьшению времени измерения теплофизических свойств в стационарном режиме с десятков минут до десятков секунд.
- Элемент Пельтье позволяет измерять термоэлектрическим методом теплопроводность образцов размерами от нескольких миллиметров.
- Элемент Пельтье обладает высокой чувствительностью к перепаду температур.
- Элемент Пельтье обладает высоким быстродействием (малое время установления теплового потока и температуры поверхности).
- егко определяется величина теплового потока, создаваемого элементом Пельтье, с помощью измерения электрического тока, проходящего через элемент и поперечного сечения элемента.
В четвёртой главе приводится описание приборов КИТ 02 - и БИТ 01С, основанных на термоэлектрическом методе измерения теплофизических свойств твёрдых материалов.
Измеритель теплопроводности КИТ 02Ц, основанный на термоэлектрическом методе контроля теплофизических параметров, является малогабаритным измерительным модулем, подсоединяемым к ПК через USB интерфейс. Прибор предназначен для измерения теплопроводности композиционных образцов из металлов, металлокерамики и других материалов в диапазоне от 0,02 до 400 Вт/ (мК) в различных модификациях.
Компьютерный измеритель теплопроводности КИТ 02 - модификации Базальт позволяет измерять образцы диаметром 15 и выстой 0,5 5 мм, с относительной погрешностью 5 % в диапазоне от 0,02 до 5 Вт/(мК). Погрешность определялась путём серии измерений на эталонных образцах и сравнения показаний прибора с сертифицированными значениями эталонных образцов. Время измерения в стандартных условиях составляет 30 40 секунд.
Сканирующий измеритель теплоёмкости БИТ 01С позволяет измерять температурозависимый коэффициент теплоёмкости образцов в форме дисков и параллелепипедов различных габаритов, в диапазоне температур от Ц50 до +50 С. Прибор БИТ - 01С показал относительную погрешность измерения на эталонных образцах не более 3 %. Время проведения измерения в диапазоне температур от Ц50 до +50 С с шагом сканирования 10 С, составляет 15 - 20 минут.
ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ РАБОТЫ.
Предложен нетрадиционный подход к измерению теплофизических свойств твёрдых материалов и новый термоэлектрический метод измерения теплового сопротивления на основе эффектов Пельтье и Зеебека, позволяющий повысить точность и сократить время измерений; доказана перспективность использования этого метода в теплофизических измерениях.
Предложена теплоизмерительная система элемент Пельтье - образец - адиабатическая камера.
Получена математическая модель термоэлектрических процессов систем элемент Пельтье - образец - теплоприёмник и элемент Пельтье - образец - адиабатическая камера.
Экспериментально показано, что конструкции приборов, основанных на термоэлектрическом методе контроля теплофизических свойств твердых материалов упрощаются, и время измерения сокращается до 30 - 40 секунд в диапазоне 0,02 - 400 Вт/(мК).
Осуществлена реализация термоэлектрического метода контроля теплофизических свойств твердых материалов в приборах КИТ - 02 - и БИТ - 01С.
ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ ОТРАЖЕНО
В ПЕЧАТНЫХ РАБОТАХ.
- Беляев Ю.И., Гринюк А.В., Внуков А.В. Компьютерный измеритель теплопроводности в диапазоне температур 20...700 С // Датчики и Системы, 2009 №10. С. 3236.
- Беляев Ю. И., Гринюк А. В., Сканирующий измеритель теплоёмкости // Приборы, №12(138) 2011. С. 13 - 16.
- Беляев Ю.И., Вепренцева О.Н., Гринюк А.В., Двикат Басам Приборы неразрушающего теплофизического контроля - Тез. докл. IX научн.-техн. конф. молодых учёных, аспирантов и студентов. Ч. II. Новомосковск, Новомосковский институт. 2007. С. 104.
- Беляев Ю.И., Вепренцева О.Н., Гринюк А.В. Применение быстродействующего измерителя теплоёмкости, Сб. трудов XX Междун. научн. конф. "Математические методы в технике и технологиях" (ММТТ-20). Т. 7. Ярославль, 2007. С. 43 - 44.
- Беляев Ю.И., Вепренцева О.Н., Гринюк А.В. Экспресс метод измерения теплоёмкости, Сборник науч. трудов. XX1. Международной конференции молодых ученых по химии и химической технологии МКХТ 2007 Москва Т. 1 2007. - С. 24-26
- Гринюк А.В., Беляева Е.Ю, Гринюк О.Н., Беляев Ю.И. Техническое обеспечение в теплофизическом эксперименте XII научно-техническая конференция молодых учёных, аспирантов, студентов. Тезисы докладов. Часть II / ГОУ ВПО РХТУ им. Д.И. Менделеева, Новомосковский институт (филиал), 2010. 228 с.; С. 7.
- Беляев Ю.И., Гринюк А.В., Иванков В.И Исследование термических свойств материалов на основе комплексного учета эффектов Пельтье и Зеебека, Вестник МАСИ. Информатика, экология, экономика. том 13 Часть 1/ Международная Академия Системных Исследований, М:, 2011. 188С., 131.
- Беляев Ю. И., Гринюк А. В., Способ определения термических свойств материалов Труды НИ РХТУ им. Д.И. Менделеева. Сер. Кибернетика, автоматизация, математика, информатизация 2011, 6(25), С., 148 - 149.
- Беляев Ю.И., Гринюк А.В. Быстродейсвующий измеритель теплоёмкости БИТ-01С, Сб. трудов XX Междун. научн. конф. "Математические методы в технике и технологиях" (ММТТ-24). Т. 7. Пенза. 2011. С. 52 - 53.
- Гринюк А.В., Беляев Ю.И., Беляева Е.Ю. Динамический метод измерения теплоёмкости XIII научно-техническая конференция молодых учёных, аспирантов, студентов. Тезисы докладов. Часть II / ГОУ ВПО РХТУ им. Д.И. Менделеева, Новомосковский институт (филиал), 2011. 228 с.; С. 45.
Форм. бум. 60х84 1/16. Объём 0,93 усл. п. . Уч.-изд. . 1,0.
Печать - ризограф. Тираж 100 экз.
МГУИЭ, издательский центр
105066, Москва, ул. Старая Басманная, 21/4
Авторефераты по всем темам >>
Авторефераты по техническим специальностям