Найбільш характерними експлуатаційними режимами роботи автомобілів являються холостий хід, прискорення, усталений режим, уповільнення. При експлуатації автомобілів у міських умовах численними дослідженнями встановлено такий розподіл режимів роботи автомобілів (табл. 1).

З таблиці видно, що співвідношення режимів роботи в загальному балансі часу несприятливе з погляду витрати палива й токсичності відпрацьованих газів. У Харкові тривалість роботи на окремих режимах для легкових автомобілів виглядає так: холостий хід - 17 %, прискорення - 41 %, постійна швидкість - 17 %, уповільнення - 25 %. При цьому швидкості руху в основному лежать у діапазоні 25...55 км/год.

Динамічна залежність комплексного критерію якості забезпечується вектором вагових коефіцієнтів k. Визначення вектора вагових коефіцієнтів критерію якості є важливим завданням, тому що дозволяє перерозподіляти пріоритети між елементами, що входять до функціоналу якості й відповідно визначають характер поведінки системи залежно від умов експлуатації автомобіля.

Визначення значень вектора вагових коефіцієнтів основане на розв'язанні задачі прийняття рішень. Про поточну умову експлуатації автомобіля можна судити за параметрами його руху - швидкість, навантаження, частота й величина прискорень та уповільнень. Як вхідні значення цих параметрів використовуються сигнали датчиків швидкості, положення дросельної заслінки й прискорення.

Рішення задачі прийняття рішень для визначення поточної умови експлуатації автомобіля за методом головного критерію показало, що найбільш адекватну оцінку поточній умові експлуатації дає середня швидкість автомобіля.

Ідентифікація значень вагових коефіцієнтів побудована на основі системи нечіткого виводу, що відповідно до конкретної умови експлуатації автомобіля встановлює певні значення вагових коефіцієнтів. Для побудови системи нечіткого виводу була прийнята така класифікація умов руху: екстенсивний міський, міський, магістральний, з підвищеною потужністю.

Екстенсивний міський рух має на увазі рух з великою щільністю транспортного потоку (затори, година пік, об'їзд ремонтних ділянок), у такому русі переважає режим холостого ходу, діапазон швидкості руху 0...25 км/год. Міський рух відповідає руху із середньою щільністю транспортного потоку, тут є в наявності всі режими руху автомобіля з їхньою регулярною зміною, швидкісний діапазон 25...60 км/год. У магістральному русі переважає режим руху зі сталою швидкістю й повністю відсутній холостий хід, діапазон швидкостей 60...100 км/год. Рух на підвищеній потужності є окремим випадком руху й має на увазі режим прискорення при обгоні або набір швидкості перед подоланням підйому.

Для нормалізації вагових коефіцієнтів їхні значення зведені до одиниці:

; ; ;

де aNe, age, aE - значення коефіцієнтів на виході системи нечіткого виводу.

Залежність вагових коефіцієнтів комплексного критерію від умови експлуатації автомобіля показано на рис. 1.

Математична модель ДВЗ представляє деякий оператор

,

який кожному вхідному вектору

,

де в - кут положення дросельної заслінки;

n - частота обертання колінчатого вала ДВЗ;

И - кут випередження запалювання;

fi - частота іскроутворення;

tr- тривалість розряду;

; ; ; ; ,

ставить у відповідність вихідний вектор

,

тобто

.

Через складність математичного опису робочих процесів ДВЗ у вигляді системи диференціальних рівнянь запропонована математична модель описує статичні характеристики ДВЗ. У цьому випадку розмірність вхідного вектора моделі скорочується:

,

а математична модель набуває вигляду

.

Тоді функціонал якості залежить від параметрів САУЗ і не залежить від функції часу. Зведемо до вигляду:

,

де - вектор вагових коефіцієнтів.

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

Завдання оптимального управління полягає в пошуку такого вектора u2*=[И* fi* tr*], що мінімізує величину J у всьому діапазоні змін вектора u1=[в n] при заданих обмеженнях, тобто

.

Поставлену задачу було розділено на дві самостійні підзадачі. Розв'язання першої полягало у визначенні залежностей, що описують взаємозв'язок параметрів вхідного й вихідного вектора . Другої - рішення задачі оптимізації .

Для визначення необхідних залежностей був розроблений метод двоетапної ідентифікації нелінійних залежностей за допомогою систем нечіткого виводу. На першому етапі виконувалася структурна ідентифікація. Вона являє собою формування нечіткої бази знань, що грубо відображає взаємозв'язок "входи-вихід" за допомогою лінгвістичних правил. Ці правила генерує експерт або отримують у результаті екстракції нечітких знань із експериментальних даних. На другому етапі відбувається параметрична ідентифікація досліджуваної залежності шляхом знаходження таких параметрів нечіткої моделі, які мінімізують відхилення результатів моделювання від експериментальних даних.

Чисельне рішення задачі оптимізації проводилося методом градієнтного спуску, шляхом підстановки вхідного вектора u1*, сформованого у вигляді числового ряду

Для кожного набору (в, n), було визначено такий набір значень (ДИ, fi, tr), при якому функціонал якості J має найменше значення.

Результат рішення задачі оптимізації наведено на рис. 2-4. Кожен параметр представлено трьома поверхнями для різних наборів вагових коефіцієнтів комплексного критерію: з пріоритетом за паливною економічністю, з пріоритетом за екологічною безпекою і з пріоритетом за потужністю.

У третьому розділі викладено методики експериментальних досліджень, описані установки й апаратура для проведення експериментів, наведені тарування приладів та оцінка похибок вимірювання, викладено методики обробки експериментальних даних.

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

Серія експериментів проводилася за єдиною спеціально розробленою методикою. Вірогідність результатів перевірялася за показниками контрольно-вимірювальних приладів і комп'ютерної системи збору даних. У процесі експериментальних випробувань проводився вибір та обґрунтування параметрів системи, за критерієм найбільшого впливу на зміну основних паливно-екологічних показників автомобіля. Проводився порівняльний аналіз результатів, отриманих при штатній системі запалювання та САУЗ. За результатами експериментальних даних проводилося довизначення й навчання нечітких систем математичної моделі.

Для забезпечення проведення експериментальних досліджень було вирішено ряд практичних завдань:

- сконструйовано експериментальний зразок системи адаптивного управління запалюванням, ґрунтуючись на аналізі схем існуючих СЗ і розроблено спосіб формування іскрового розряду;

- розроблено комп'ютеризований моторний стенд для проведення експериментальних досліджень;

- проведено оцінку похибок застосовуваної вимірювальної апаратури;

- підготовлено математичний апарат для обробки результатів експериментів;

- розроблено програмно-апаратний комплекс для проведення дорожніх випробувань САУЗ на автомобілі.

Для проведення експериментальних досліджень багатопараметричного управління запалюванням автором була виготовлена система управління запалюванням, що реалізовує розроблений спосіб управління іскровим розрядом. Функціонально вона складається із двох блоків: блок управління робить опитування датчиків і обчислює поточні значення параметрів розряду, а силовий модуль формує серії імпульсів необхідної інтенсивності й тривалості і розподіляє їх по циліндрах. САУЗ забезпечує такі характеристики:

- інтенсивність іскроутворення 4…36 кГц;

- тривалість розряду 0,2…3 мс;

- точність встановлення КВЗ 0,05 град. п.к.в.;

- дискретність встановлення КВЗ 0,1 град. п.к.в.;

- розрахункова енергія розряду 10…250 мДж.

Дослідження проводилися моторним методом на комп'ютеризованому моторному стенді, розробленому на базі двигуна МеМЗ-307 у лабораторії кафедри ДВЗ ХНАДУ. Стенд обладнаний приладами й пристроями, необхідними для контролю режимів і показників роботи двигуна. Необхідна точність вимірювання (ДСТУ 148 46-81) забезпечується установкою комплексу сучасної вимірювальної апаратури.

Система автоматичного збору даних дозволяє в реальному масштабі часу реєструвати більше двадцяти параметрів, що характеризують хід протікання експерименту. Для реєстрування тиску в циліндрі ДВЗ розроблений спеціальний програмно-апаратний вимірювальний комплекс. Комплекс створено на основі п'єзоелектричного датчика тиску 8QP500ca фірми AVL, пристрою формування сигналів синхронізації й плати аналогового вводу-виводу L-783.

Для розрахунку індикаторних показників циклу було розроблено програму, що забезпечує пакетну обробку даних програмно-апаратного комплексу. Програма автоматично розраховує індикаторну роботу серії підряд записаних циклів і за середнім значенням визначає оцінну індикаторну діаграму, для якої проводиться подальший розрахунок.

Дорожні випробування проводилися на автомобілі ЗАЗ «Таврія», так як модель його силового агрегату (МеМЗ-307) відповідає двигуну, встановленому на моторному стенді, для якого проводилася адаптація параметрів САУЗ. Для вимірювання екологічних і паливно-економічних параметрів автомобіль був оснащений п'ятикомпонентним газоаналізатором «Інфракар 5М2Т» та системою автоматичного збору даних. Інформація із приладів, у реальному масштабі часу, записувалася на жорсткий диск комп'ютера для подальшої обробки.

У четвертому розділі наведено програму експериментальних досліджень, показані результати стендових випробувань, підготовлені й проведені дорожні випробування за європейським їздовим циклом, проаналізовано дані експериментів.

Програма експериментальних досліджень включала:

1) стендові випробування:

- вплив КВЗ на параметри процесу згоряння;

- вплив інтенсивності іскроутворення на параметри процесу згоряння;

- вплив тривалості розряду на параметри процесу згоряння;

2) дорожні випробування під час руху за європейським їздовим циклом;

3) порівняльний аналіз експериментальних даних.

Установлено, що збільшення КВЗ та інтенсивності іскроутворення приводить до зменшення тривалості фази видимого згоряння, що спричиняє підвищення коефіцієнта активного тепловиділення наприкінці фази згоряння. Це веде до зростання температури в процесі згоряння й зниження ефективності згоряння палива у фазі догоряння, внаслідок чого збільшується вміст неспалених вуглеводнів і викидів NOx.

При підвищеній інтенсивності іскроутворення САУЗ надає достатні впливи для формування стійкого первинного осередку горіння бідних робочих сумішей. З'являється можливість забезпечити термодинамічно оптимальний кут випередження запалювання навіть при несприятливих умовах протікання робочого процесу. Підвищення інтенсивності іскроутворення при зменшенні КВЗ дозволяє компенсувати негативний вплив малих КВЗ на потужнісні й паливно-економічні показники. При цьому температура не досягає критичних значень, час процесу згоряння скорочується, що дозволяє знизити емісію NOx і CH.

Перевірка вірогідності оптимальних значень адаптивної настройки параметрів САУЗ, що відповідають різним експлуатаційним умовам, полягала в порівняльному аналізі паливно-економічних та екологічних показників автомобіля, отриманих при дорожніх випробуваннях. Дорожні випробування були побудовані за європейським їздовим циклом випробувань автомобілів на паливну економічність та екологічну безпеку.

Експеримент складався із трьох заїздів за європейським їздовим циклом. Перший проводився для штатної настройки системи запалювання (рис. 5), другий - з оптимізацією параметрів САУЗ за критерієм якості із пріоритетом за паливною економічністю (рис. 6) й третій для настройки САУЗ із пріоритетом за екологічною безпекою (рис. 7).

Порівняльний аналіз інтегральних значень паливно-економічних та екологічних показників автомобіля відносно штатної настройки системи запалювання показав: для САУЗ настроєної із пріоритетом за паливною економічністю, витрата палива знижується на 5…7 %, NOx на 7…10 %, а викиди CH збільшуються на 3…4 %. При настройці САУЗ із пріоритетом за екологічною безпекою викиди CH і NOx знижуються на 10…12 % і 20…24 % відповідно при збільшенні витрати палива на 3…4 %.

Запропонований метод адаптивної настройки параметрів СЗ може бути використаний як при розробці алгоритмів управління нових електронних систем управління двигуном (ЕСУД), так і для оптимізації роботи існуючих систем. Для сучасних автомобілів з ЕСУД встановлення САУЗ не передбачає змін конструкції чи встановлення додаткових датчиків.

ВИСНОВКИ

1. На основі аналізу науково-технічної літератури встановлено, що сучасні СЗ не можуть забезпечити оптимізацію управління робочим процесом за паливно-екологічним комплексним критерієм, тому що мають тільки один керований параметр - кут випередження запалювання. При цьому значення оптимальних КВЗ для паливно-економічних показників істотно відрізняються від значень оптимальних КВЗ із погляду екологічної безпеки автомобіля. У даній роботі запропоновано концепцію багатопараметричної системи управління запалюванням, що передбачає одночасне управління декількома параметрами іскрового розряду й забезпечує можливість гнучкої настройки системи запалювання.

2. Теоретичні дослідження енергетичних параметрів іскрового розряду показали, що фаза пробою має найбільший коефіцієнт передачі енергії запалювання порівняно з іншими фазами розряду. Найбільші значення потужності розряду, що визначають запальну здатність розряду, також проявляються у фазі пробою. На підставі проведених досліджень був розроблений спосіб формування іскрового розряду, що дозволяє оперативно змінювати момент, інтенсивність і тривалість розряду. В основі способу лежить метод частотно-імпульсного управління.

3. Проведено експериментальні дослідження впливу параметрів іскрового розряду на показники процесу згоряння двигуна. Дослідження показали істотний вплив інтенсивності розряду на швидкість згоряння, це дозволяє знижувати токсичність викидів ВГ та підтримувати паливно-економічні показники на досить високому рівні при зменшенні КВЗ. Так при зменшенні КВЗ на 10…15 град. п.к.в. і збільшенні інтенсивності й тривалості розряду до 24…32 кГц і 1,5 мс відповідно, викиди NOx і CH знижуються до 42% і 26% при цьому втрати паливної економічності не перевищують 3%.

4. Розроблено метод синтезу математичної моделі об'єкта управління, особливість якого полягає в об'єднанні методик моделювання класичної оптимальної САУ й САУ на основі нечіткої логіки й експертних систем, а також у локалізації зони пошуку оптимуму. Це дозволяє знизити обсяг математичних обчислень і використовувати градієнтні методи оптимізації.

5. Визначено комплексний критерій якості для параметричної оптимізації САУЗ. Як параметри комплексного критерію використовуються потужнісні, паливно-економічні й екологічні показники автомобіля. Це дозволяє налаштовувати САУЗ на конкретні умови експлуатації автомобіля за допомогою зміни вектора вагових коефіцієнтів k. Визначено залежність вагових коефіцієнтів від умов експлуатації автомобіля. Для проведення експериментальних досліджень розраховано субоптимальні значення параметрів САУЗ для трьох різних наборів вагових коефіцієнтів функціонала якості: з пріоритетом за паливною економічністю, з пріоритетом за екологічною безпекою та з пріоритетом за потужністю.

6. Проведено дорожні випробування автомобіля з адаптивною настройкою системи запалювання під час руху за європейським їздовим циклом. Порівняльний аналіз результатів показав, що при настройці САУЗ із пріоритетом за паливною економічністю, витрата палива знижується на 5…7 %, NOx на 7…10 %, а викиди CH збільшуються на 3…4 %. При настройці САУЗ із пріоритетом за екологічною безпекою викиди CH і NOx знижуються на 10…12 % і 20…24 % відповідно при збільшенні витрати палива на 3…4 %.

Список опублікованих праць за темою дисертації

1. Автоматизированный стенд для исследования и доводки газового малолитражного быстроходного двигателя внутреннего сгорания / А.Н. Туренко, Ф.И. Абрамчук, А.Н. Пойда, А.А. Дзюбенко [и др.] // Автомобильный транспорт: сборник научных трудов. - Х.: ХНАДУ, 2008. - Вып. 23. - С. 89-94. - здобувачем розроблено та реалізовано систему автоматичного збору даних, запропоновано спосіб варіації кута випередження запалювання мікропроцесорної системи управління двигуном.

2. Дзюбенко А.А. Управление параметрами искрового разряда системы зажигания ДВС / А.А. Дзюбенко // Вісті Автомобільно-дорожнього інституту: науково-виробничий збірник. - Горлівка: АДІ ДонНТУ. - 2009. - №1(8). - С. 69-73.

3. Сериков С.А. Идентификация математической модели двигателя внутреннего сгорания с использованием системы нечеткого вывода / С.А. Сериков, А.А. Дзюбенко // ДВС: Всеукраинский научно-технический журнал. - Харьков: НТУ «ХПИ». - 2009. - №1. - С. 14-18. - здобувачем запропоновано метод ідентифікації математичної моделі двигуна внутрішнього згорання з використанням системи нечіткого виводу.

4. Дзюбенко А.А. Параметрическая оптимизация системы автоматического управления зажиганием / А.А. Дзюбенко // Автомобильный транспорт: сборник научных трудов. - Харьков: ХНАДУ. - 2009. - Вып. 25. - С. 13-17.

5. Микропроцессорная система управления двигателем внутреннего сгорания: каталог аннотаций на разработки по материалам первого Украинско-китайского форума «Наука - производство», 23-25 мая 2007 г. / А.В. Бажинов, Ю.Н. Бороденко, А.А. Дзюбенко, С.А. Сериков. - Х.: ХНУРЕ, 2007. - С. 158. - здобувачем запропоновано схему мікропроцесорної системи для управління процесами запалювання.

6. Дзюбенко А.А. Разработка многоискровой системы зажигания с регулированием энергии разряда / А.А. Дзюбенко // ЕКОВАРНА 2008: сборник докладов / XIV научно-техническая конференция с международным участием «Транспорт, экология - устойчивое развитие» 8-10 мая 2008. - Варна: ТУ-Варна, 2008. - С. 258-263.

паливний економічність автомобіль транспортний

Анотація

Дзюбенко О.А. Поліпшення паливно-економічних та екологічних показників автомобіля на основі оптимізації параметрів системи запалювання. - Рукопис.

Дисертація на здобуття наукового ступеня кандидата технічних наук за спеціальністю 05.22.20 - експлуатація та ремонт засобів транспорту. - Харківський національний автомобільно-дорожній університет, Харків, 2010 р.

Дисертацію присвячено обґрунтуванню та розробці методу адаптивної настройки системи управління запалюванням залежно від експлуатаційних умов автомобіля.

В дисертації наведено теоретичне узагальнення та нове рішення задачі поліпшення паливно-економічних та екологічних показників автомобіля. Задача рішається шляхом адаптивної настройки системи запалювання залежно від експлуатаційних умов автомобіля на основі оптимізації параметрів системи управління запалюванням за комплексним критерієм.

Проведено теоретичні та практичні дослідження впливу параметрів іскрового розряду системи запалювання на робочий процес двигуна, за результатами яких запропоновано концепцію багатопараметричного управління запалюванням. Розроблено математичну модель об'єкта управління, що дозволяє отримувати оптимальні параметри СЗ залежно від комплексного критерію, який визначається умовами експлуатації автомобіля. Розроблено метод адаптивної настройки СЗ під конкретні умови експлуатації автомобіля та експериментально доведено доцільність його використання.

Ключові слова: система запалювання, експлуатація, комплексний критерій, адаптивна настройка, паливно-екологічні показники.

АННОТАЦИЯ

Дзюбенко А.А. Улучшение топливно-экономичных и экологических показателей автомобиля на основе оптимизации параметров системы управления зажиганием. - Рукопись.

Диссертация на соискание научной степени кандидата технических наук по специальности 05.22.20 - эксплуатация и ремонт средств транспорта. - Харьковский национальный автомобильно-дорожный университет, Харьков, 2010 г.

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

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

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

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

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

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

Ключевые слова: система зажигания, эксплуатация, комплексный критерий, адаптивная настройка, топливно-экологические показатели.

Abstract

A.A. Dziubenko. Improvement of Vehicle's Fuel-Saving and Ecological Indexes on Basis of Ignition System Parameters Optimization. - Manuscript.

Thesis for competition of scientific degree of candidate of technical science on speciality 05.22.20 -operation and maintenance of vehicle. - Kharkiv National Automobile and Highway University, Kharkiv, 2010.

The thesis deals with the grounding and development of adaptive tuning method of ignition control system depending on vehicle's operation conditions.

The theoretical generalization and the new task decision concerning the improvement of vehicle's fuel-saving and ecological factors are presented in the given work. The task is settled by means of adaptive tuning of ignition control system for a particular mode of vehicle's operation on the basis of parameters optimization of ignition operation system according to the complex criterion. A new method for ignition system spark discharge parameters operation is developed. The theoretical and practical researches of ignition system spark discharge influence on engine's working process, according to the results of which the concept of multivariate ignition operation is offered, are developed.

The mathematical model of management object that allows to obtain IS optimal parameters depending on the complex criterion which is determined by the vehicle's operation mode is developed. The IS adaptive tuning method for particular vehicle's traction mode is developed, and the expediency of its application is experimentally proved.

Key words: ignition system, vehicle's operation mode, complex criterion, adaptive tuning, fuel-ecological indexes.

Підписано до друку 05.05.2010 р. Формат 6090 1/16. Папір офсетний.

Гарнітура Times New Roman Cyr . Віддруковано на ризографі

Ум. друк. арк. 0,9.

Зам. № 344/10. Тираж 100 прим. Ціна договірна

Видавництво

Харківського національного автомобільно-дорожнього університету

Видавництво ХНАДУ, 61200, Харків-МСП, вул. Петровського, 25.

Тел. /факс: (057)700-38-64; 707-37-03

Свідоцтво Державного комітету інформаційної політики, телебачення та радіомовлення України про внесення суб'єкта видавничої справи до Державного реєстру видавців, виговидавничої продукції, серія виготівників і розповсюджувачів видавничої продукції, серія ДК №897 від 17.04 2002 р.

Размещено на Allbest.ru