WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

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

МИХАЙЛОВ Владимир Сергеевич

МЕТОДИКА ПАРАМЕТРИЧЕСКОЙ

ИДЕНТИФИКАЦИИ МОДЕЛИ ПРОЦЕССОВ

ГАЗООБМЕНА ДВУХТАКТНЫХ ДВС

Специальность: 05.04.02 - Тепловые двигатели

АВТОРЕФЕРАТ

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

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

Уфа - 2010

Работа выполнена на кафедре «Двигатели внутреннего сгорания»

ГОУ ВПО государственный авиационный технический «Уфимский университет»

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

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

Еникеев Рустэм Далилович

Официальные оппоненты: доктор технических наук, профессор Неговора Андрей Владимирович, каф. тракторы и автомобили, Башкирский государственный аграрный университет доктор технических наук, профессор Гимранов Эрнст Гайсович, каф. прикладной гидромеханики, Уфимский государственный авиационный технический университет ГСКБ «Трансдизель (г. Челябинск)

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

Защита диссертации состоится «25» февраля 2011 г. в 10-00 часов на заседании диссертационного совета Д-212.288.05 при Уфимском государственном авиационном техническом университете по адресу: 450000, г. Уфа, ул. К. Маркса, 12, УГАТУ, зал заседаний ученого совета (1 корпус) Телефон (347) 273-77-92, факс (347) 272-29-18, e-mail: admin@ugatu.ac.ru

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО «Уфимский государственный авиационный технический университет»

Автореферат разослан « 11 » января 2011 г.

Ученый секретарь диссертационного совета д-р техн. наук, профессор Ф.Г. Бакиров

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Поэтому развитие методик применения одномерных моделей для описания протекающих в ГВТ процессов представляется важным и актуальным.

В настоящем исследовании предложена и реализована методика параметрической идентификации одномерной рациональной модели течения рабочего тела в ГВТ 2-тактных ПДВС. Определен количественный уровень и качественный характер отклонений моделирования показателей этого процесса.

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

Цель и задачи исследования. Целью исследования является разработка методики параметрической идентификации модели газообмена в ГВТ 2-тактных ПДВС, включающей этапы идентификации и верификации идентифицированной модели.

В соответствии с целью формулировались задачи исследования:

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

2. Исследовать и доказать возможность параметрической идентификации модели процессов в ГВТ 2-тактных ПДВС по предложенной методике, используя для этого экспериментально определенные величины интегральных показателей и индикаторные диаграммы базовых образцов исследуемых двигателей.

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

Методы исследования. При выполнении работы использованы следующие методы исследования:

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

- экспериментальные исследования на полноразмерных ПДВС;

- методы решения задач многопараметрической оптимизации.

Исследование носило расчетно-экспериментальный характер. Методика параметрической идентификации одномерной модели процессов газообмена в ГВТ 2-тактных ПДВС проверялась на достоверность сравнением прогнозируемых расчетных показателей с экспериментальными данными.

Экспериментальные данные получены измерениями на полноразмерных ПДВС.

Научная новизна 1. Впервые разработана методика параметрической идентификации рациональной одномерной модели процессов в ГВТ ПДВС, в которой в качестве варьируемых параметров выступают параметры замыкающих соотношений вида уравнений в обобщенных переменных.

2. Теоретически обоснована правомерность использования интегральной характеристики продувки-наполнения рабочей камеры для параметрической идентификации модели газообмена в ГВТ 2-тактного ПДВС.

3. На классе 2-тактных ПДВС с волновыми эффектами в ГВТ доказано, что методика позволяет в несколько раз повысить точность расчетных интегральных показателей двигателей как для базового варианта ГВТ, так и для варианта с существенно отличными конструктивными размерами.

Практическая ценность Разработанная методика, а также результаты расчетных исследований внедрены в промышленности – ОАО УМПО (г. Уфа) и в учебный процесс УГАТУ (г. Уфа).

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

1. Повысить точность прогнозирования показателей двигателя при сквозном моделировании по одномерной модели процессов в ГВТ.

2. Использовать предложенную методику в рамках расчетных работ, направленных на параметрический синтез ГВТ.

3. Улучшать характеристики 2-тактных ПДВС при проектировании и доводке для удовлетворения требованиям к их форсированности и экономичности.

На защиту выносится:

1. Методика параметрической идентификации модели процессов в ГВТ ПДВС, в которой в качестве варьируемых выступают параметры зависимостей в обобщенных переменных, привлекаемых для замыкания уравнений рациональной одномерной модели процессов в ГВТ.

2. Использование параметров продувочной характеристики, замыкающей двухзонную модель течения в рабочей камере в период продувки, в качестве варьируемых параметров для параметрической идентификации модели процессов в ГВТ 2-тактных ПДВС.

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

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

- корректном использовании фундаментальных уравнений механики жидкости и газа, теории рабочих процессов ДВС и рациональных одномерных моделей газообмена;

- сопоставлении результатов расчетов с данными экспериментов на реальных ПДВС в стендовых условиях.

Апробация работы. Результаты работы докладывались и обсуждались на пятой всероссийской зимней школе-семинаре аспирантов и молодых ученых «Актуальные проблемы науки и техники» (г. Уфа, 2010), на всероссийской молодежной научной конференции «Мавлютовские чтения» (г. Уфа, 2009, 2010), на межрегиональной научно-технической конференции на тему «Повышение эффективности колесных и гусеничных машин многоцелевого назначения (г. Челябинск, 2010), на XXXVI международной молодежной научной конференции «Гагаринские чтения» (г. Москва, 2010), на Юбилейной научно-технической конференции «Двигатель - 2010» (г. Москва, 2010), на межотраслевой научно-технической конференции «Актуальные проблемы развития поршневых ДВС» (г. Санкт-Петербург, 2010).

Личный вклад соискателя в разработку проблемы: все основные идеи работы сформулированы лично автором. Материалы диссертации основаны на исследованиях автора в период с 2005 по 2010 годы.

Публикации. Основное содержание диссертации опубликовано в печатных работах, в том числе в 3 публикациях в центральных рецензируемых журналах, рекомендованных ВАК.

Объем и структура работы. Диссертация состоит из введения, четырёх глав, основных выводов, приложения и списка литературы. Содержит страниц машинописного текста, включающего 123 рисунка, 8 таблиц и библиографический список из 83 наименований.

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

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

В первой главе проведен анализ работ, связанных с темой научного исследования. Рассмотрены и проанализированы научные труды, посвященные проблемам моделирования рабочих процессов и процессов газообмена 2-тактных ПДВС таких ученых, как А.С. Орлин, М.Г. Круглов, Б.П. Рудой, Ю.А. Гришин, Н.В. Лобов, А.А. Черноусов, а также ряда зарубежных исследователей, таких, как G.P. Blair, T. Petrie, H.S. Ricardo, C.F. Caunter, C.F.

Taylor, H. List, M. Kadenacy и др. Проведен обзор и анализ моделей, а также способов моделирования рабочих процессов и процессов газообмена в ГВТ ПДВС.

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

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

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

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

Выше изложенное позволяет выдвинуть гипотезу, что для 2-тактного двигателя решение задачи синтеза оптимального ГВТ, обеспечивающего желаемое массовое наполнение, можно рационализировать, опираясь все же на быстросчетную одномерную модель процессов в ГВТ и применяя методику параметрической идентификации используемой модели по набору параметров, выбранному таким образом, чтобы идентифицированную модель можно было бы (с «равномерно» уменьшенной ошибкой моделирования) применять к объекту, с большим успехом предлагая удачные конструктивные мероприятия для достижения требуемых показателей изделия. Основой для этой гипотезы служат данные о том, что соответствующие вспомогательные модели нестационарного течения по длинному каналу и квазистационарного течения через местное сопротивление весьма адекватны описываемым явлениям.

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

Условием сохранения массы компонента k в произвольном объеме является уравнение (для k = 1, …, K):

С учетом массовых сил закон сохранения количества движения для смеси можно записать в виде:

где Div П ij – операция дивергенции тензора, аk = аk (r, t) – ускорение массовой силы для k-го компонента.

Уравнение сохранения энергии движущейся смеси в интегральной форме имеет вид:

теплопроводности Фурье:

где = (, T, Y1,..., YK 1 ) – коэффициент теплопроводности смеси.

Течение рабочего тела в тепловом двигателе происходит, как правило, в турбулентном режиме. Сквозной расчет течения в ГВТ теплового двигателя на ЭВМ непосредственно по законам сохранения, дополненным общепринятыми моделями кинетики и молекулярного переноса, в принципе возможен. Такой подход к численному расчету турбулентного течения называется прямым численным моделированием (Direct Numerical Simulation, DNS). Однако такой подход неприменим для практических расчетов рабочих процессов тепловых двигателей, так как требует вычислительных ресурсов, превышающих возможности современных супер-ЭВМ.

Первыми моделями, позволившими рассчитывать турбулентные течения при решении инженерных задач, стали модели на основе уравнений НавьеСтокса, осредненных по Рейнольдсу (Reynolds Averaged Navier-Stokes, RANS).

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

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

Описаны модели, выведенные из законов сохранения, справедливых в одномерном приближении, описывающие рабочие процессы и процессы газообмена ПДВС. Условия сохранения масс К компонентов смеси и внутренней энергии смеси в переменном объеме выражаются уравнениями:

В приведенной системе законов сохранения учтены: объемное тепловыделение от внешних источников энергии dQвнеш / dt на m кг смеси, в форме теплоты по механизму теплообмена со стенкой с температурой Tw (t ), и обмен с внешней средой энергией в форме работы сжатия – расширения объема V (t ).

Уравнения (1) и (2) образуют систему обыкновенных дифференциальных уравнений, которые используются для описания процессов в элементах ГВТ ДВС, таких как рабочая камера, кривошипная камера, ресивер.

Представленные модели пониженной пространственной размерности, которые были разработаны на кафедре ДВС УГАТУ, позволяют описывать протекающие в рабочей камере процессы с учетом смены заряда и переработки исходных компонентов рабочего тела в продукты сгорания. В работе используется представление рабочего тела в виде смеси двух компонентов – свежего заряда и продуктов сгорания.

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

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

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

На втором этапе методики производится параметрическая идентификация модели газообмена по условию соответствия выбранному набору данных об объекте. Осуществляется обоснованный выбор варьируемых параметров идентификации. В силу рациональности применяемой модели варьируемыми параметрами могут быть только те, что определяют протекание эмпирических зависимостей, служащих замыкающими модельными соотношениями для элементов ГВТ ПДВС. Использование в качестве варьируемых параметров характеристик моделей граничных сечений недопустимо. Модель выгорания сама может быть «откалибрована» путем параметрической идентификации, и если далее процесс газообмена смоделировать в точности, итоговая модель процессов в ГВТ, предположительно, будет выдавать весьма близкую к измеренной индикаторную мощность. Показано, что в случае 2-тактного двигателя для использования в качестве варьируемых остаются параметры модели продувки, а именно – продувочной характеристики, замыкающей двухзонную модель и позволяющей с более высокой точностью рассчитывать как массовое наполнение, так и потери свежей смеси при продувке рабочей камеры. Именно данная зависимость была выбрана для решения обратной задачи подбора варьируемых параметров. В качестве таковых, были приняты семь параметров продувочной характеристики (рис. 1): x0 – точка отрыва кривой (предельный относительный объем зоны свежей смеси, при котором из выпускных окон истекает газовая смесь только из зоны продуктов сгорания), и координаты трех промежуточных точек кривой для интерполяции продувочной характеристики многочленом четвертой степени. Предполагается, что такое сочетание параметров продувочной характеристики обеспечит необходимое для параметрической идентификации разнообразие вариантов ее протекания.

Рисунок 1 - Выбор варьируемых параметров продувочной характеристики Оптимальное протекание продувочной характеристики определяется путем идентификации по результатам измерений на двигателе в условиях стенда.

В качестве целевой функции берется сумма квадратов относительных разностей расчетных и экспериментальных значений N i и Gв в точках ВСХ:

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

Рисунок 2 - Схема расчета по одномерной модели с использованием методики Таким образом, в работе решается задача идентификации модели рабочего процесса в ГВТ двигателя, выполняемая для повышения точности и надежности предсказываемых моделированием параметров ДВС.

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

Предлагаемая методика параметрической идентификации проверена на двух 2-тактных двигателях с выраженными волновыми явлениями на выпуске:

АПД-800 – двигатель с принудительным воспламенением и внешним смесеобразованием разработки Уфимского моторостроительного производственного объединения, и ЯАЗ-М204А – двигатель с воспламенением от сжатия и внутренним смесеобразованием, производства Ярославского моторного завода.

В третьей главе представлены результаты стендовых испытаний двигателя АПД-800, оборудованного выпускными системами, настроенными на различные диапазоны оборотов. Следовало подтвердить, что использование предложенной методики параметрической идентификации модели газообмена, повышает точность расчета показателей моделируемого 2-тактного ПДВС. С этой целью измерены показатели расхода воздуха Gв (n) и индикаторной мощности N i (n) при работе по ВСХ, выполнено индицирование двигателя.

Для проведения экспериментов использована исследовательская установка (рис. 3, а) в составе двигателя АПД-800 с настроенной выпускной системой, электромеханической динамометрической тормозной машины АРА 100 фирмы AVL и специализированного оборудования фирмы AVL для выполнения индицирования рабочей камеры. Все измерения соответствуют требованиям ГОСТ 14846–81. Измерения проведены на двигателе с настроенной выпускной системой базового варианта конструкции (рис. 3, б), используемой для идентификации модели, и с выпускной системой, используемой для верификационных расчетов (рис. 3, в).

Рисунок 3 - а) двигатель АПД-800 в составе исследовательской установки;

б) выпускная система для идентификации; в) выпускная система для верификации В результате проведенных испытаний получен набор данных, по которым построены кривые изменения давления в рабочей камере двигателя при работе по ВСХ для вариантов с выпускной системой базовой конструкции и в комплектации с выпускной системой, используемой для верификации.

Для построения зависимости индикаторной мощности от оборотов N i (n) выполнялась обработка результатов измерений с помощью специально разработанной программы indic. Эта программа позволяет обрабатывать результаты измерения индикаторного давления в рабочей камере в каждой точке ВСХ и получать характеристику индикаторной мощности двигателя.

Значения индикаторной работы циклов Li определялись для каждой точки зависимости N i (n), причем каждая зависимость индикаторного давления в рабочей камере от угла поворота коленчатого вала (ПКВ) представлена как средняя по пятидесяти циклам. Численное интегрирование велось с шагом в один градус ПКВ методом прямоугольников второго порядка с учетом номинальной кинематики кривошипно-шатунного механизма исследуемого двигателя. Значение давления p РК ( ) определялось линейной интерполяцией без какого-либо сглаживания сигнала. За начало цикла принимался момент нахождения поршня в нижней мертвой точке, получаемый автоматически с помощью датчика положения коленчатого вала AVL 365C-SET-MULTIPLIER.

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

Структурные схемы модели ГВТ создавались в графической среде системы имитационного моделирования (СИМ) Horsepower Lab 1D с учетом их основных особенностей, и задавались исходные данные, включая характеристики элементов ГВТ. Данная СИМ была использована для всех численных расчетов процессов в ГВТ.

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

Полученные из расчета показатели исследуемого двигателя (рис. 4, а) сравниваются с характеристиками, полученными на двигателе в условиях стенда. Из представленных результатов видно, что величина отклонения расчетных показателей в некотором диапазоне оборотов имеет неприемлемое значение, а именно, отклонение расчетной характеристики от измеренной достигает для расхода воздуха 60%, а для индикаторной мощности 40%. Такая величина отклонения при выполнении мероприятий по доводке двигателя недопустима.

На следующем этапе выполнялся расчет с применением методики идентификации модели ГВТ выбранного объекта исследования, при этом использованы результаты измерений двигателя, установленного на стенд в варианте, когда он был укомплектован выпускной системой, используемой для идентификации модели. Полученные зависимости представлены на рис. 4, б.

Рисунок 4 - Расчетные и экспериментальные значения интегральных показателей:

а) расчет без идентификации модели; б) расчет с идентификацией модели идентифицированной модели газообмена для двигателя, работающего по ВСХ.

Для оценки величины отклонения моделирования, показатели ДВС, полученные в результате этих расчетов сравниваются с показателями, измеренными на двигателе с другой выпускной системой (рис. 5).

Рисунок 5 - Расчетные и экспериментальные значения интегральных Анализ полученных результатов показывает, что применение методики параметрической идентификации повышает точность моделирования параметров ПДВС. Погрешность расчета индикаторной мощности двигателя АПД-800 снижена в четыре раза, а характеристики расхода воздуха – в шесть раз (рис. 4). Расчеты с выпускной системой, настроенной на другой диапазон оборотов (верификационный расчет), показывают, что величина отклонения моделируемых интегральных показателей от измеренных в условиях стенда не превышает 10% (рис. 5).

В четвертой главе представлены результаты стендовых испытаний двигателя ЯАЗ-М204А, также оборудованного двумя (идентификационной и верификационной) выпускными системами. Исследование проводилось с целью подтверждения эффективности методики параметрической идентификации для расчетов 2-тактных ПДВС с внутренним смесеобразованием, для чего выполнялись измерения показателей работы двигателя по ВСХ (расхода воздуха Gв (n) и индикаторной мощности N i (n) ) и индицирование рабочей камеры.

Для проведения экспериментов была собрана исследовательская установка (рис. 6, а) в составе двигателя ЯАЗ-М204А и электромеханической динамометрической тормозной установки 1DS 541 N фирмы VSETIN. В рабочем состоянии была оставлена ближняя к валу отбора мощности секция двигателя, прочие секции двигателя выключены из работы. Выпускные каналы сообщены с тупиковым коллектором, выпускные клапаны извлечены, а каналы в направляющих втулках заглушены, насос-форсунки отсоединены от подающей и сливной топливных магистралей и от топливной рейки, и зафиксированы в закрытом положении. Все измерения соответствуют требованиям ГОСТ 14846–81. Единственная рабочая секция ЯАЗ-М204А снабжалась специально спроектированными трубопроводами (рис. 6, б и в), используемыми для идентификации и верификации соответственно.

Рисунок 6 - а) двигатель ЯАЗ-М204А в составе исследовательской установки;

б) выпускная система для идентификации; в) выпускная система для верификации В результате проведенных измерений получены данные, по которым построены характеристики изменения давления в рабочей камере двигателя при работе по ВСХ для вариантов с выпускной системой, используемой для идентификации, и в комплектации с выпускной системой, используемой для верификации.

Построение зависимостей индикаторных мощностей от оборотов N i (n) для двигателя ЯАЗ-М204А выполнялось так же, как и для двигателя АПД-800, с помощью программы indic.

В ходе проведения подготовительного этапа моделирования произведен расчет характеристик граничных сечений ГВТ в пакете FIRE, а также интегральной характеристики газообмена в рабочей камере. После этого создавалась структурная схема модели ГВТ двигателя ЯАЗ-М204А в СИМ Horsepower Lab 1D, задавались исходные данные и проводился расчет.

Первоначально проводилось выполнение прямого расчета по одномерной модели, где для замыкания процессов в рабочей камере выполнялась идентификация закона выгорания. Полученные расчетом показатели (рис. 7, а) сравниваются с характеристиками, измеренными на двигателе в стендовых условиях. Из полученных результатов видно, что величина отклонения расчетных показателей в рабочем диапазоне оборотов имеет неприемлемое значение, а именно, отклонение расчетной характеристики расхода воздуха от измеренной достигает 20%, а характеристики индикаторной мощности - 30%.

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

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

Рисунок 7 - Расчетные и экспериментальные значения интегральных показателей:

а) расчет без идентификации модели; б) расчет с идентификацией модели.

идентифицированной модели газообмена для двигателя, работающего по ВСХ.

Для оценки величины отклонения моделирования показатели ДВС, полученные в результате этих расчетов, сравниваются с измеренными на двигателе с другой выпускной системой (рис. 8).

Рисунок 8 - Расчетные и экспериментальные значения интегральных В результате проведенного расчетно-экспериментального исследования показано, что применение методики параметрической идентификации для 2-тактных двигателей с внутренним смесеобразованием позволяет повысить точность расчета параметров ПДВС. Для двигателя ЯАЗ-М204А величина погрешности расчета индикаторной мощности снижена в шесть раз, а расхода воздуха (в рабочем диапазоне оборотов) – в два раза (рис. 7). При этом величина отклонения моделируемых интегральных показателей в рабочем диапазоне оборотов (при изменении условий однозначности, в частности, с другой выпускной системой) относительно измеренных в условиях стенда, не превышает 10% (рис. 8).

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

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

Результаты выполненного исследования позволяют сделать следующие выводы:

1. Обоснована, сформулирована и реализована методика моделирования процессов в ГВТ двухтактных ДВС с волновыми явлениями в выпускной системе, направленная на повышение точности результатов расчетов. В основу методики положена параметрическая идентификация модели течения при продувке рабочей камеры ДВС по известным (измеренным) значениям интегральных показателей и индикаторных диаграмм базового образца ДВС.

2. Проведена параметрическая идентификация моделей процессов в ГВТ двухтактных ПДВС по предложенной методике. Выполнены измерения соответствующих показателей двух полноразмерных ПДВС указанного класса, различающихся способом смесеобразования и продувки. Установлено, что отклонения интегральных параметров расхода воздуха и индикаторной мощности от измеренных не превысили 10% на участке ВСХ, подвергнутом идентификации, тогда как базовая модель давала в несколько раз большие отклонения. Таким образом, в обоих случаях применение методики идентификации модели позволило получить достаточно достоверные (для двухтактных ДВС с настроенным ГВТ) результаты.

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

Таким образом, в работе сформулирована, реализована и проверена методика параметрической идентификации одномерной модели процессов в ГВТ двухтактных ПДВС с волновыми явлениями в выпускной системе.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ

В СЛЕДУЮЩИХ РАБОТАХ:

Статьи в рецензируемых журналах, рекомендованных ВАК Рособрнадзора:

1. Михайлов В.С. Газообмен и эффективные показатели ДВС с двухконтурной системой впуска / Р.Д. Еникеев, В.С. Михайлов // Вестник УГАТУ. 2007. Т 9, № 6 (24). С. 82–97.

2. Михайлов В.С. Методика параметрического синтеза ГВТ ДВС по рациональным быстросчетным моделям процессов газообмена/ В.С. Михайлов // Вестник Академии Военных Наук, № 1 (30) (спецвыпуск), Москва, 2010.

С. 126–132.

3. Михайлов В.С. Расчетно-экспериментальное совершенствование процессов газообмена двигателя АПД-800/ Р.Д. Еникеев, В.Ю. Иванов, В.С. Михайлов, В.Ф. Нурмухаметов, С.П. Павлинич, А.А. Черноусов // Вестник УГАТУ. 2010. Т.14, № 2 (37). С. 13–20.

В других изданиях, включая труды Всероссийских и международных НТК:

4. Михайлов В.С. Опыт применения расчетных продувочных характеристик при моделировании газообмена двухтактных двигателей/ А.А. Черноусов, В.С. Михайлов // Вопросы теории и расчета рабочих процессов тепловых двигателей: Межвузовский научный сборник. Уфа. Выпуск 22, УГАТУ, 2008. С. 237-253.

5. Михайлов В.С. Расчетное исследование путей совершенствования процессов газообмена двухтактного двигателя АПД-800/ В.С. Михайлов // Мавлютовские чтения: Материалы конференции. Уфа, УГАТУ. 2009. Т.1.

С. 123-124.

6. Михайлов В.С. Методика параметрического синтеза газовоздушных трактов ДВС по одномерной модели процессов газообмена/ В.С. Михайлов // Актуальные проблемы науки и техники: Сборник трудов пятой всероссийской зимней школы-семинара аспирантов и молодых учёных. Уфа, УГАТУ. 2010.

Том 2. С. 262–266.

7. Михайлов В.С. Методика параметрического синтеза газовоздушных трактов двухтактных ДВС по одномерной модели процессов газообмена/ В.С. Михайлов // XXXVI Гагаринские чтения: Международная молодежная научная конференция. Москва, МАТИ. 2010. Т.2. С. 178–180.

8. Михайлов В.С. Методика параметрической идентификации газовоздушных трактов двухтактных ДВС по одномерной модели процессов газообмена/ В.С. Михайлов // Мавлютовские чтения: Материалы конференции.

Уфа, УГАТУ. 2010. Т.1 – с. 81 – 82.

9. Михайлов В.С. Параметрическая идентификация одномерной модели газообмена в газовоздушном тракте двухтактных двигателей/ В.С. Михайлов // Двигатель - 2010: Сборник научных трудов международной конференции, посвященной 180-летию МГТУ им. Н.Э. Баумана. Москва. МГТУ. 2010.

С. 239-243.

10. Михайлов В.С. Методика параметрической идентификации модели процессов газообмена в ГВТ двухтактных ДВС/ В.С. Михайлов // Актуальные проблемы развития поршневых ДВС: Материалы межотраслевой научнотехнической конференции. Санкт-Петербург. СПбГМТУ. 2010. С. 68-69.

МЕТОДИКА ПАРАМЕТРИЧЕСКОЙ

ИДЕНТИФИКАЦИИ МОДЕЛИ ПРОЦЕССОВ

ГАЗООБМЕНА ДВУХТАКТНЫХ ДВС

Специальность: 05.04.02 - Тепловые двигатели диссертации на соискание ученой степени Подписано в печать 27.12.2010. Формат 60 х 84 1/ Бумага офсетная. Печать плоская. Гарнитура Times New Roman.

Усл. печ. л. 1,0. Усл. кр. - отт. 1,0. Уч. - изд. л. 0,9.

ГОУ ВПО Уфимский государственный авиационный Центр оперативной полиграфии 450000, Уфа – центр, ул. К.Маркса,

 


Похожие работы:

«Демьянов Владимир Александрович РАСЧЕТНО-ЭКСПЕРИМЕНТАЛЬНОЕ ОБОСНОВАНИЕ КОНСТРУКЦИИ ЭКОЛОГИЧНЫХ ПОВОРОТНО - ЛОПАСТНЫХ ГИДРОТУРБИН Специальность 05.04.13 - Гидравлические машины и гидропневмоагрегаты АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург - 2013 Работа выполнена в ОАО Силовые машины. Научный руководитель : доктор физико-математических наук, профессор, член - корреспондент РАН, Петреня Юрий Кириллович. Официальные оппоненты...»

«Лыков Алексей Викторович ВЫБОР И РАСЧЕТНОЕ ОБОСНОВАНИЕ ХАРАКТЕРИСТИК УТИЛИЗАЦИОННОЙ ПАРОТУРБИННОЙ УСТАНОВКИ ДЛЯ ВЫРАБОТКИ ЭЛЕКТРОЭНЕРГИИ НА СОБСТВЕННЫЕ НУЖДЫ ГАЗОПЕРЕКАЧИВАЮЩИХ КОМПРЕССОРНЫХ СТАНЦИЙ Специальность 05.04.12 – Турбомашины и комбинированные турбоустановки АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург – 2014 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального...»

«МУКТАРОВ Орынгали Джулдгалиевич ВЛИЯНИЕ РЕЖИМОВ ОБРАБОТКИ УСКОРЕННЫМИ ИОНАМИ АЗОТА НА СТРУКТУРУ ЭЛЕКТРОПЛАЗМЕННОГО ПОКРЫТИЯ И ФИЗИКО-МЕХАНИЧЕСКИЕ СВОЙСТВА ТИТАНА Специальности 05.09.10 – Электротехнология 05.02.07 – Технология и оборудование механической и физико-технической обработки АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Саратов – 2013 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего...»

«МАЖИТОВ БАХРИДДИН ЖАМИЛОВИЧ МЕТОДЫ СНИЖЕНИЯ ТОКСИЧНОСТИ ОТРАБОТАВШИХ ГАЗОВ ДИЗЕЛЯ И ТЕПЛОНАГРУЖЕННОСТИ ТОРМОЗНОЙ СИСТЕМЫ АВТОМОБИЛЕЙ ПРИ ЭКСПЛУАТАЦИИ В ГОРНЫХ УСЛОВИЯХ Специальности: 05.04.02 - Тепловые двигатели 05.22.10 - Эксплуатация автомобильного транспорта АВТОРЕФЕРАТ Диссертация на соискание ученой степени кандидата технических наук Санкт-Петербург 2011 2 Работа выполнена в ФГОУ ВПО Санкт – Петербургский государственный аграрный университет Научные руководители:...»

«Антоненков Максим Александрович ОБОСНОВАНИЕ ПРОЕКТНЫХ РЕШЕНИЙ ГИДРОСТАТИЧЕСКИХ ПОДШИПНИКОВ ГЛАВНЫХ ЦИРКУЛЯЦИОННЫХ НАСОСОВ РЕАКТОРОВ НА БЫСТРЫХ НЕЙТРОНАХ, ОХЛАЖДАЕМЫХ СВИНЦОВЫМ И СВИНЕЦ-ВИСМУТОВЫМ ТЕПЛОНОСИТЕЛЯМИ Специальность 05.04.11 – Атомное реакторостроение, машины, агрегаты и технология материалов атомной промышленности АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Нижний Новгород 2013 Работа выполнена на кафедре Атомные, тепловые станции...»

«КАСАТКИНА Елена Геннадьевна ПОВЫШЕНИЕ КАЧЕСТВА ПЛАТИНИТА СОВЕРШЕНСТВОВАНИЕМ ТЕХНОЛОГИИ ЕГО ПРОИЗВОДСТВА Специальность 05.02.23 – Стандартизация и управление качеством продукции (металлургия) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Магнитогорск – 2006 2 Работа выполнена в ГОУ ВПО Магнитогорский государственный технический университет им. Г.И. Носова Научный руководитель доктор технических наук, профессор Гун Геннадий Семенович Официальные...»

«БУЯНКИН ПАВЕЛ ВЛАДИМИРОВИЧ ОЦЕНКА УСТОЙЧИВОСТИ ПЛАТФОРМ И НАГРУЗОК В ОПОРНО-ПОВОРОТНЫХ УСТРОЙСТВАХ ЭКСКАВАТОРОВ-МЕХЛОПАТ Специальность 05.05.06 – Горные машины Автореферат диссертации на соискание ученой степени кандидата технических наук Кемерово - 2014 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Кузбасский государственный технический университет имени Т. Ф. Горбачева. Научный руководитель - доктор...»

«Галатов Кирилл Станиславович СОВЕРШЕНСТВОВАНИЕ ОБОРУДОВАНИЯ ДЛЯ РАЗДЕЛЕНИЯ ПЕРОПУХОВОГО СЫРЬЯ НА ФРАКЦИИ С РАЗРАБОТКОЙ УЗЛА МЕХАНИЧЕСКОЙ ДЕФОРМАЦИИ СТЕРЖНЯ КУРИНОГО ПЕРА Специальность 05.02.13 – Машины, агрегаты и процессы (коммунальное хозяйство и сфера услуг) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Шахты – 2013 2 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования...»

«Гайнов Алексей Александрович ПОВЫШЕНИЕ НАДЕЖНОСТИ И РЕСУРСА РАБОЧИХ ЭЛЕМЕНТОВ СУДОВЫХ ГАЗОТРУБНЫХ КОТЛОВ Специальность: 05.08.05 – Судовые энергетические установки и их элементы (главные и вспомогательные) Автореферат диссертации на соискание ученой степени кандидата технических наук Нижний Новгород – 2011 Работа выполнена в Федеральном бюджетном образовательном учреждении высшего профессионального образования Волжская государственная академия водного транспорта (ВГАВТ)...»

«Савченко Андрей Владимирович СОВЕРШЕНСТВОВАНИЕ СКВАЖИННОЙ ТЕХНОЛОГИИ ГИДРОИМПУЛЬСНОГО ВОЗДЕЙСТВИЯ НА ГОРНЫЕ ПОРОДЫ ПРИ ДОБЫЧЕ ПОЛЕЗНЫХ ИСКОПАЕМЫХ Специальность: 25.00.22 – Геотехнология (подземная, открытая и строительная) 05.05.06 – Горные машины АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2009 Работа выполнена в Учреждении Российской академии наук Институте горного дела Сибирского отделения РАН академик РАН, профессор Научный...»

«УДК 621.981.1 Гудков Иван Николаевич РАЗРАБОТКА ТЕХНОЛОГИИ ИЗГОТОВЛЕНИЯ ТОНКОСТЕННЫХ ПЕРФОРИРОВАННЫХ ПРОФИЛЕЙ МЕТОДОМ ИНТЕНСИВНОГО ДЕФОРМИРОВАНИЯ Специальность: 05.03.05 – Технологии и машины обработки давлением АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Нижний Новгород – 2009 Работа выполнена в Ульяновском государственном техническом университете Научный руководитель : доктор технических наук, профессор Филимонов Вячеслав Иванович...»

«ЯБЛОНЕВ АЛЕКСАНДР ЛЬВОВИЧ ОСНОВЫ ПРОГНОЗИРОВАНИЯ ВЗАИМОДЕЙСТВИЯ ПНЕВМАТИЧЕСКОГО КОЛЕСНОГО ХОДА С ТОРФЯНОЙ ЗАЛЕЖЬЮ Специальность 05.05.06 Горные машины Автореферат диссертации на соискание ученой степени доктора технических наук Тверь 2011 2 Работа выполнена на кафедре Торфяные машины и оборудование ФГБОУ ВПО Тверской государственный технический университет. Научный консультант : Доктор технических наук, профессор Зюзин Борис Федорович Официальные оппоненты : Доктор...»

«ХИСАЕВА ЗЕМФИРА ФАНИЛОВНА ПОВЫШЕНИЕ СТОЙКОСТИ МЕТАЛЛА ПЕЧНЫХ ТРУБ К КОКСООТЛОЖЕНИЮ СИЛИЦИРОВАНИЕМ ПОВЕРХНОСТИ Специальность 05.02.01 – Материаловедение (Машиностроение в нефтегазовой отрасли) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Уфа - 2003 Работа выполнена на кафедре Машины и аппараты химических производств Уфимского государственного нефтяного технического университета. Научный руководитель доктор технических наук, профессор Кузеев...»

«Башаров Рашит Рамилович СОВЕРШЕНСТВОВАНИЕ ТЕХНОЛОГИИ КОНЦЕВОГО ФРЕЗЕРОВАНИЯ С УЧЁТОМ УПРУГИХ ОТЖАТИЙ ВЫСОКОЧАСТОТНОГО ЭЛЕКТРОШПИНДЕЛЯ СТАНКА 05.02.07 – Технология и оборудование механической и физико-технической обработки Автореферат диссертации на соискание ученой степени кандидата технических наук Оренбург 2012 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Уфимский государственный авиационный...»

«Костюк Инна Викторовна МЕТОДЫ ОЦЕНКИ ЭФФЕКТИВНОСТИ ТЕХНОЛОГИИ АДАПТИВНОГО РАСТРИРОВАНИЯ Специальность 05.02.13 – Машины, агрегаты и процессы (печатные средства информации) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва – 2010   Работа выполнена на кафедре Технологии допечатных процессов в ГОУ ВПО Московский государственный университет печати. Научный руководитель : доктор технических наук, профессор Кузнецов Юрий Вениаминович Официальные...»

«Гаврилов Илья Юрьевич ОПРЕДЕЛЕНИЕ ВЛИЯНИЯ НАЧАЛЬНОГО СОСТОЯНИЯ ПАРА НА ВОЛНОВУЮ СТРУКТУРУ И ПАРАМЕТРЫ ДВУХФАЗНОГО ПОТОКА В СОПЛОВОЙ ТУРБИННОЙ РЕШЕТКЕ Специальность 05.04.12 – Турбомашины и комбинированные турбоустановки АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва – 2014 Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Национальный исследовательский университет...»

«ФИГУРА КОНСТАНТИН НИКОЛАЕВИЧ ОБОСНОВАНИЕ КОНСТРУКТИВНЫХ ПАРАМЕТРОВ И РЕЖИМОВ РАБОТЫ СМЕСИТЕЛЬНЫХ АГРЕГАТОВ С ВНУТРЕННИМИ ВИБРОАКТИВАТОРАМИ Специальность: 05.02.13 – Машины, агрегаты и процессы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Братск 2013 Работа выполнена в ФГБОУ ВПО Братский государственный университет Научный руководитель : кандидат технических наук, доцент, заведующий кафедрой Подъемно-транспортные, строительные, дорожные машины...»

«• Щербаков Виталий Сергеевич НАУЧНЫЕ ОСНОВЫ ПОВЫШЕНИЯ ТОЧНОСТИ РАБОТ. ВЫПОЛНЯЕМЫХ ЗЕМЛЕРОЙНО-ТРАНСПОРТНЫМИ МАШИНАМИ 05.05.04 - Д о р о ж н ы е и с т р о и т е л ь н ы е м а ш и н ы Автореферат д и с с е р т а ц и и на с о и с к а н и е у ч е н о й с т е п е н и доктора технических наук О м с к - 2000 Г у? у 9 Работа выполнена в Сибирской государственной автомобильно-дорожной академии (СибАДИ) Официальные оппоненты : доктор технических наук, профессор Абрамснко в Э.А.; доктор...»

«МОСКОВКО Юрий Георгиевич МЕТОДИКА ПРОЕКТИРОВАНИЯ И РАЗРАБОТКА ЭНЕРГОЭФФЕКТИВНЫХ ОСЕВЫХ ВЕНТИЛЯТОРОВ С ПРОФИЛЯМИ ЛОПАТОК СПЕЦИАЛЬНОЙ ФОРМЫ Специальность: 05.04.06 - Вакуумная, компрессорная техника и пневмосистемы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург- 2011 Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Санкт-Петербургский государственный...»

«ОСИПОВ Александр Вадимович ПОВЫШЕНИЕ ЭКОНОМИЧНОСТИ ДВУХСТУПЕНЧАТОГО ОТБОРНОГО ОТСЕКА ПАРОВОЙ ТУРБИНЫ Специальность 05.04.12 Турбомашины и комбинированные турбоустановки Автореферат диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург 2002 Работа выполнена в Брянском государственном техническом университете. Засл. деятель науки и техники РФ, Научный руководитель доктор техн. наук, профессор В.Т. Буглаев. Официальные оппоненты : – Засл. деятель...»














 
© 2013 www.diss.seluk.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Методички, учебные программы»

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