WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


 

       На правах рукописи

СЕРГЕЕВ Сергей Викторович

МЕТОДИКА ИССЛЕДОВАНИЯ ТЕПЛОВЫХ ХАРАКТЕРИСТИК
ВЛАГОСОДЕРЖАЩИХ МАТЕРИАЛОВ НА ОСНОВЕ РЕШЕНИЯ
ОБРАТНОЙ ЗАДАЧИ ТЕПЛОПРОВОДНОСТИ

Специальность: 01.04.14 – «Теплофизика и теоретическая теплотехника»

Автореферат

Диссертации на соискание ученой степени

кандидата технических наук

Санкт-Петербург

2012

Работа выполнена на кафедре физики ГОУ ВПО «Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий»

Научный руководитель:

доктор технических наук, профессор

Платунов Евгений Степанович

Официальные оппоненты:

доктор технических наук, профессор

Федоров Александр Валентинович

кандидат технических наук, доцент

Старков Александр Сергеевич

Ведущая организация:

Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский институт метрологии им. Д. И. Менделеева», Санкт-Петербург

Защита состоится “___” ______________ 2012 г. в _____ часов на заседании диссертационного Совета Д 212.234.01 при Санкт-Петербургском государственном университете низкотемпературных и пищевых технологий по адресу: 191002, Санкт-Петербург, ул. Ломоносова, д. 9, СПбГУНиПТ.

Тел./факс: (812) 315-30-15

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

       Автореферат разослан «__» ____________ 2012 г.

Ученый секретарь диссертационного Совета,

доктор технических наук, профессор                                                Рыков В. А

Актуальность работы. Информация о теплофизических характеристиках как функциях температуры играет ключевую роль в инженерных расчетах тепловых процессов в самых разных отраслях. Особый интерес представляют дисперсные влагосодержащие материалы, обладающие важной спецификой. Они меняют свои свойства  необратимо, и особенно резко это происходит с повышением температуры. Именно поэтому для влагосодержащих материалов практически важной становится возможность комплексного измерения ТФХ как функций температуры. Для этого были разработаны комплексные динамические методы измерений, базирующиеся на теории монотонного режима, которая является обобщением квазистационарных и регулярных методов.

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

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

Использование методов некорректных обратных задач для расчетов ТФХ влагосодержащих материалов особенно актуально, поскольку они являются практически единственным инструментом, позволяющим получать информацию о ТФХ в «проблемных» диапазонах температур, где традиционные методы исчерпывают свои возможности. Указанный подход перспективен еще и с той точки зрения, что он практически не требует специальных режимов проведения эксперимента и может снять ряд ограничений как при проведении опыта, так и в самой конструкции теплоизмерительной ячейки. Это объясняется тем, что сама методика предполагает решение уравнения теплопроводности в нелинейной двумерной постановке и не накладывает ограничений на характер изменения граничных условий, в отличие от традиционной, где имеет место аналитическое решение простейших задач теплопроводности и необходимость строгого выдерживания режима близкого к стационарному за счет сохранения небольших перепадов температуры внутри образца.

Создание рассмотренной выше методики позволит в перспективе, используя «традиционную» теплоизмерительную ячейку для комплексного измерения ТФХ, полностью автоматизировать обработку экспериментальных данных и расчет коэффициентов ТФХ, не требуя промежуточных трудоемких операций с участием исследователя, а вместо этого предоставляя возможность контроля и гибкого управления алгоритмами вычислений. При этом значительно повысится экспрессность эксперимента в целом. Кроме того, совместное использование двух методик расчета коэффициентов (традиционной аналитической и предлагаемой в данной работе) позволит получать наиболее объективную картину.

Целью работы являлась разработка новой экспрессной методики для комплексного исследования температурных зависимостей теплопроводности и объемной теплоемкости влагосодержащих материалов (пищевых продуктов, грунтовых и горных пород) на основе решения нелинейной инверсной коэффициентной обратной задачи теплопроводности по данным эксперимента в условиях их замораживания в диапазоне температур (–60…20) °С.

       Для достижения поставленной цели необходимо было решить следующие задачи:

  1. Обосновать необходимость разработки нового численного алгоритма расчета коэффициентов ТФХ влагосодержащих материалов при замораживании в диапазоне температур, охватывающем неизотермический фазовый переход.
  2. Разработать математическую модель теплоизмерительной ячейки, используемой для исследования твердых, рыхлых, пастообразных и жидких материалов.
  3. Разработать численный алгоритм расчета коэффициентов эффективной теплопроводности и объемной теплоемкости на основе решения обратной задачи теплопроводности.
  4. Исследовать численным методом свойства нескольких характерных влагосодержащих материалов.

Новизну работы составляют и выносятся на защиту следующие результаты и положения:

  1. Обоснование преимуществ новой методики расчета коэффициентов ТФХ влагосодержащих материалов в условиях замораживания по сравнению с приближенно аналитическим методом.
  2. Численный алгоритм расчета комплекса ТФХ влагосодержащих материалов в диапазоне температур, охватывающем фазовые и структурные превращения.
  3. Определение точности предложенной методики.
  4. Результаты расчета численным методом объемной теплоемкости и теплопроводности группы пищевых продуктов и грунтовых пород.

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

Личный вклад: разработка методики проведения опыта при последовательном замораживании и размораживании образца, обработки экспериментальных данных опыта; проведение комплексного измерения ТФХ пищевых продуктов и грунтовых пород в области температур (–60…20) °С; разработка численного алгоритма расчета теплоемкости и теплопроводности на основе решения обратной задачи теплопроводности.

Апробация работы и публикации. Основные материалы и результаты работы опубликованы в 5 печатных трудах, в том числе 2 работы в изданиях рекомендуемых ВАК РФ, и доложены на 34-ой, 35-ой и 36-ой НПК профессорско–преподавательского состава докторантов, аспирантов и сотрудников университета, Санкт-Петербург, (2008…2010 гг.); V Международной научно-технической конференции «Низкотемпературные и пищевые технологии в XXI веке», Санкт-Петербург, 2010 г.; Международная научно-техническая конференция «Современные методы и средства исследований теплофизических свойств веществ», Санкт-Петербург, 2010 г.

Объем работы. Диссертация состоит из введения, четырех глав, заключения, списка литературы (130 наименований) и 2 приложений. Работа изложена на 120 страницах машинописного текста, содержит 2 таблицы и 28 рисунков.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Для определения ТФХ влагосодержащих материалов используют динамические методы свободного охлаждения. Экспериментальная часть метода представляет собой ячейку, которая выполнена в виде многослойного сосуда, состоящего из двух вставленных друг в друга металлических стаканов 2 и 4, разделенных тонкостенной теплоизоляционной прослойкой 3. Наружный стакан 2 имеет плоское основание, которое в опыте контактирует с массивным металлическим блоком 1 через теплоизоляционную прослойку 11. Его наружная боковая поверхность окружена теплозащитной оболочкой 10. В полость внутреннего стакана 4 перед опытом вставляется по скользящей посадке съемная тонкостенная металлическая ампула 5 с находящимся в ней образцом 9.

Ампула с образцом после установки в ячейку закрывается металлической крышкой 7. К наружному стакану сверху крепится теплоизоляционное кольцо 6. На нем в опытах размещается тонкостенный прозрачный полимерный колпак 8. Крышка 7 в центре имеет направляющее отверстие для наконечника центральной термопары, регистрирующей в опытах температуру образца в его центральной зоне. Наконечник второй термопары, регистрирующей температуру внутреннего стакана 4, размещается в отверстии его стенки. Наконечник третьей термопары, предназначенной для регистрации температуры наружного стакана, монтируется в том же отверстии.

                       а)                                                                б)

Рис. 1

а) схема тепловой ячейки; б) характерный вид зависимости температуры от времени в опыте с дистиллированной водой

Температуры в опыте регистрируются с интервалом в 5–10 с. В результате получается два массива экспериментальных данных температуры центра и поверхности образца. После сглаживания имеем две экспериментальные кривые температура – время, рис. 1б.

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

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

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

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

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

Рис.3 Перепады температур внутри различных образцов в процессе опыта.

Для определения удельной теплоемкости образца используется формула :

, (1)

где М – масса образца; - скорость изменения среднеобъемной температуры образца; – суммарный тепловой поток, проходящий через тепломер; K(tя) – эффективная тепловая проводимость теплоизоляционной прослойки, отнесенная условно к температуре ядра; ϑт(τ) – перепад температуры между ядром и блоком.

Для определения теплопроводности образца используют соотношение:

.  (2)

Анализ кривых ТФХ как функций температуры, рассчитанных указанным методом, выявляет физически необъяснимое поведение в частности кривой коэффициента теплопроводности в диапазоне температур фазового перехода, выражающееся в резком (почти дельта-образном) возрастании с последующим большим разбросом значений, утрачивающим всякую закономерность. Наиболее явно это выражено у влагосодержащих материалов с неизотермическим фазовым переходом (рис. 4). Причем подобная ситуация наблюдалась при исследовании различных образцов.

На рис. 4 приведен характерный вид рассчитанного аналитическим методом коэффициента теплопроводности, когда «проблемные значения» удалены из рассмотрения. При этом гладким участком кривой показано наиболее вероятное реальное поведение зависимости в диапазоне -5 °С+5 °С. Подобное предположение легко объясняется тем, что значение эффективной теплопроводности двухфазного материала не может превышать большего значения одной из фаз в отдельности.

Рис. 4. Зависимость эффективных теплофизических характеристик в опыте с яблоком

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

Можно выделить несколько основных этапов в реализации предлагаемого метода расчета коэффициентов ТФХ:

  1. Построение математической модели теплоизмерительной ячейки
  2. Оценка адекватности модели теплоизмерительной ячейки с ее тепловым режимом, который имеет место в процессе опыта.
  3. Оценка традиционной методики расчета коэффициентов ТФХ с помощью предлагаемой модели.
  4. Постановка коэффициентной обратной задачи теплопроводности.
  5. Расчет ТФХ некоторых характерных образцов в рамках решения обратной задачи.
  6. Сопоставление обоих методов расчета ТФХ

На первом этапе строится математическая модель реальной теплоизмерительной ячейки, работающей в составе стенда для исследования ТФХ материалов динамическим методом. Исследуемый образец влагосодержащего материала представляется изотропным телом в форме ограниченного цилиндра высотой и радиусом и имеющий до начала опыта температуру . C началом опыта () образец начинает охлаждаться, через его основание и боковую поверхность проходит тепловой поток плотностью , при этом температура боковой поверхности и в центре образца изменяется с течением времени. Функции и измеряются в процессе эксперимента.

Сформулированная задача сводится к решению квазилинейного двумерного уравнения теплопроводности:

.  (6)

С начальным условием

. (7)

И граничными условиями

, , , .  (8)

,  . (9)

Здесь – температура в точке с координатами в момент времени ; – скрытая теплота фазового перехода. – безразмерный параметр состояния (– образец полностью разморожен, – образец полностью заморожен).

Двумерную задачу теплопроводности (6) удобно решать численным методом переменных направлений, который требует постановки двух вспомогательных одномерных задач в нашем случае по радиальному и осевому направлениям. Таким образом, имеем:

Для осевого направления:

  (10)

Для радиального направления:

(11)

С граничными условиями (8), (9) и начальными условиями (7).

Параметр удобно использовать при моделировании неизотермических фазовых переходов, поскольку «размытость» теплового эффекта по температуре в этом случае определяется видом производной (см. рис. 5).

Частным случаем является изотермический фазовый переход, когда производная представляет собой – функцию Дирака. Параметр состояния при этом скачком изменяется от 0 до 1.  При рассмотрении однофазного материала указанная производная обращается в ноль, поскольку параметр остается постоянным и принимает значение 1 или 0 в зависимости от вида фазы. В этом случае в уравнении (6) в правой части пропадает последнее слагаемое.

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

       

                       а)                                                б)

Рис. 5.
а) качественный вид кривой параметра состояния ;
б) производной параметра состояния .

На втором этапе для оценки адекватности предложенной модели, рассматривается образец из полиметилметакрилата, который используется для калибровки теплоизмерительной ячейки. Тепловые свойства этого материала хорошо известны и являются табличными величинами. Используя табличные данные ТФХ указанного материала, рассчитываем температурное поле образца на основе предлагаемой математической модели теплоизмерительной ячейки. Из полученного поля температур выбираются две точки, пространственно соответствующие местам установки спаев термопар в образце. Тогда появляется возможность сопоставить первичные температуры, полученные из эксперимента с температурами, которые рассчитаны на основе модели. Данные сопоставления представлены на рис. 6.

       

Рис. 6. Сопоставление расчетных и экспериментальных данных
на примере полиметилметакрилата

На третьем этапе для оценки достоверности, рассчитанных по (1) и (2) коэффициентов ТФХ, предлагается следующий математический эксперимент. Моделируется температурное поле внутри образца в виде численного решения двумерной нелинейной задачи теплопроводности для цилиндрической области конечных размеров с известными из эксперимента граничными условиями (тепловой поток как функция времени). При этом ТФХ как функции температуры нам также уже известны (вычислены традиционным аналитическим методом). Решив такую задачу, мы получаем температурное поле, как функцию времени внутри цилиндрической области. Далее выделим из температурного поля две точки – в центре модели и на границе и получим две кривые температура – время. При сопоставлении их с экспериментальными температурными кривыми будет наблюдаться некоторое расхождение в области температуры затвердевания, вызванное неточностью вычисления коэффициентов ТФХ приближенно аналитическим методом. Используя расчетные данные и для яблок (рис. 4), вычисляются кривые температуры как функции времени.

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

Рис. 7. Сопоставление расчетных температур по данным и и экспериментальных.

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

  (10)

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

(11)

(12)

где- решение сопряженной задачи теплопроводности вида:

(13)

(14)

(15)

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

Апробирование общего метода решения производилось на модельной задаче. В частности важным оказался выбор способа минимизации функционала (10). Для минимизации функционала (10) первоначально использовались три метода: градиентного метода, метода наискорейшего спуска и метода сопряженных градиентов. Лучшие результаты дал метод сопряженных градиентов. Анализ показал, что во всех трех случаях не обеспечивается равномерная сходимость решения у границ исследуемой области. Причем улучшение сходимости не наблюдается даже при интерполяции между узловыми точками функцией более высокого порядка, чем линейная, в частности квадратичная, кубическая, сплайн. Такая ситуация имеет место и при восстановлении ТФХ по данным реального эксперимента с глиной, яблоками и мясом проведенных Тамбулатовой Е. В. Такое подведение решения может быть связано с тем, что компоненты градиента функционала (10), соответствующие граничным узлам функции в несколько раз меньше компонент соответствующим внутренним узлам. Последнее обусловлено характером температурного поля в образце, системой исходных данных и величиной шага табличного представления функций и . В результате при минимизации функционала граничные точки области движутся к решению значительно медленнее остальных, что и вызывает искажение решения у границ.

       В соответствии с задачей исследования было проведено определение ТФХ ряда наиболее характерных материалов: дистиллированная вода, полиметилметакрилат, влажная глина, а также продуктов: яблоко, мясо кабана. Регистрация временных зависимостей температур и потоков в каждом опыте осуществлялась в зависимости от специфики образца с различным интервалом опроса датчиков температур и общей длительностью эксперимента. Относительные отклонения восстановленных зависимостей от справочных данных по ТФХ полиметилметакрилата не превышали 5 % как по теплоемкости, так и по теплопроводности.

                               а)                                                        б)

в)                                                г)

                       д)                                                        е)

Рис. 8 Сопоставление зависимостей и , рассчитанных двумя методами

Коэффициенты ТФХ, рассчитанные по методике, разработанной автором, проверялись на прямой задаче теплопроводности. При этом полученные кривые температур сравнивались с первичными температурами, измеренными в опытах на соответствующих образцах.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

  1. Разработан метод расчета коэффициентов ТФХ как функций температуры образцов цилиндрической формы  по известным временным функциям теплового потока через поверхность образца и температур в центре и на поверхности цилиндра при замораживании в широкой области температур в зонах резкого изменения свойств. Метод основан на численном решении инверсной задачи теплопроводности в экстремальной постановке.
  2. Адекватность метода расчета коэффициентов и подтверждена на физических экспериментах с различными типами образцов.
  3. Разработан численный алгоритм расчета функций и в среде MathCad.
  4. Исследовано влияние систематической погрешности измерения температур и потока, а так же численного метода решения ИЗТ на точность восстановления теплофизических характеристик. Исследование проводилось путем последовательного ввода некой случайной погрешности с нормальным распределением и моделируемой амплитудой в гладкие входные данные.
  5. Проанализировано влияние параметров дискретизации области численного расчета на точность результата вплоть до определения границ сходимости решения.
  6. Построена математическая модель тепловой ячейки для расчета двухмерного температурного поля внутри образца и последующей оценки корректности найденных ранее функций и по средствам решения прямой задачи теплопроводности.

ОСНОВНЫЕ  ПУБЛИКАЦИИ

  1. Сергеев С. В., Зобков П. Н., Шрамова А. Л. Влияние способов отвода теплоты от сорбента на эффективность процесса криосорбционной откачки. «Вестник МАХ», выпуск №3, 2009 г., с. 41-43.
  2. Сергеев С. В., Баранов И. В., Куслиева Е. В., Платунов Е. С. Прибор для комплексного исследования влагосодержащих материалов. Известия Санкт-Петербургского государственного университета низкотемпературных и пищевых технологий., №1, СПб, 2009 г., с. 143-146.
  3. Сергеев С. В., Платунов Е. С. Восстановление коэффициентов теплопроводности и объемной теплоемкости как функций температуры по экспериментальным данным установки для исследования рыхлых грунтовых пород. «Вестник МАХ», выпуск № 3, 2011 г., с. 27-31.
  4. Сергеев С. В. Исследование комплекса теплофизических характеристик влагосодержащих материалов как функций температуры в зоне фазового перехода на основе решения обратной задачи теплопроводности. Сб. материалов МНТК «Современные методы и средства исследований теплофизических свойств веществ» Материалы конфренции. – СПб.: СПбГУНиПТ, 2010. с.73.
  5. Сергеев С. В. Численное моделирование процесса промерзания влагосодержащих образцов различной конфигурации в области фазового перехода. Электронный научный журнал «Холодильная техника и кондиционирование»/ ГОУ ВПО «Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий» ­– Электронный журнал ­–Санкт-Петербург: СПбГУНиПТ, 2011, выпуск №1.
 






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.