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

Размещено на http://www.allbest.ru/

ВІННИЦЬКИЙ ДЕРЖАВНИЙ ТЕХНІЧНИЙ УНІВЕРСИТЕТ

Автореферат

дисертації на здобуття наукового ступеня доктора технічних наук

ПОБУДОВА МЕТОДОЛОГІЇ МАТЕМАТИЧНОГО МОДЕЛЮВАННЯ НЕЛІНІЙНИХ НЕСТАЦІОНАРНИХ СИСТЕМ ДЛЯ ЗАДАЧ ОЦІНКИ ЇХ СТІЙКОСТІ ТА ЧУТЛИВОСТІ В ПРОСТОРІ ПРИРОСТІВ ПАРАМЕТРІВ

Спеціальність: Математичне моделювання та обчислювальні методи

ЮХИМЧУК СЕРГІЙ ВАСИЛЬОВИЧ

ВІННИЦЯ, 1999 РІК



1. ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

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

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

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

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

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

Зв'язок роботи з науковими програмами, планами, темами.

Робота виконувалась у відповідності з проектом Державного фонду фундаментальних досліджень України ”Розробка методів аналізу нелінійних систем автоматичного керування динамічними об'єктами”, затвердженому наказом ДКНТП України від 22.03.1994, №52 (№ держ. реєстр. 0194U042191), державною науково - технічною програмою 04.00 ”Екологічно чиста енергетика та ресурсозберігаючі технології” (шифр завдання 05.51.03./419 - 92, затвердженого Постановою Комітету по науково - технічному прогресу при Кабінеті Міністрів України від 04.05.1992 р., №12 та наказом ДКНТП України від 01.03.1993 р., №15, № держ. реєстр. 0193U002771), координаційним планом науково - дослідних робіт Міністерства освіти України за фаховим напрямом ”Інформатика і нові інформаційні технології”, затвердженим наказом Міністерства освіти України від 13.02.1997 р., №37 (№ держ. реєстр. 0197U012665), а також з тематичними планами проведення НДДКР у Вінницькому політехнічному інституті (ВПІ) та Вінницькому державному технічному університеті (ВДТУ) на госп. договірних засадах з ВО ”Вінницяенерго”(№ держ. реєстр. 01880044245), з Державним концерном України ”Арксі” (м. Боярка, Київська область)(№ держ. реєстр. 0194U008791), ВО ”Хімпром”(м. Вінниця) (№ держ. реєстр. 0194U025745), Ладижинською ДРЕС (м. Ладижин, Вінницької області) (№ держ. реєстр. 0195U003915), ВАТ ”Вінницягаз” (№ держ. реєстр. 0194U044166), а також з тематичними планами проведення НДДКР у ВПІ та ВДТУ за рахунок коштів державного бюджету, узгодженими Міністерством освіти України (№ держ. реєстр. 0191000320,0193U029093, 0194U016744).

Мета і задачі дослідження.

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

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

2. Побудувати методологію розв'язку задач стійкості класу систем, що розглядаються, в запропонованому просторі опису;

3. Підтвердити ефективність запропонованого підходу шляхом розв'язку задач оцінки стійкості конкретних нелінійних систем;

4. Шляхом математичного моделювання вирішити задачу оцінки стійкості та оптимізації параметрів одного широкого класу нелінійних нестаціонарних систем в електроенергетиці;

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

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

Наукова новизна одержаних результатів.

В роботі отримані наступні наукові результати:

1. Вперше запропоновано здійснювати опис поведінки нелінійних нестаціонарних систем у новому просторі - просторі приростів параметрів. Побудована методологія розв'язку задач у запропонованому просторі опису;

2. Вперше здійснено узагальнення методу описуючих функцій для однозначних та неоднозначних типових нелінійних при різних характерах зміни їх параметрів;

3. На основі отриманих результатів побудована методологія аналізу нелінійних нестаціонарних систем та отримані нові наукові результати - розв'язана зворотне задача стійкості для деяких класів нелінійних інформаційно-вимірювальних систем із змінними параметрами;

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

5. Здійснено подальший розвиток методології Фур'є - інтегрального методу ідентифікації на задачі оцінки стійкості та діагностики стану нелінійних нестаціонарних систем, що дозволило запропонувати функціональні схеми пристроїв, технічна новизна яких захищена авторськими свідоцтвами СРСР на винаходи, патентами України та Російської Федерації.

Практичне значення одержаних результатів.

Одержані наукові результати впроваджені на ТЕЦ підприємства теплових мереж (ПТМ) ПЕО ”Вінницяенерго”, на Державному концерні України ”Арксі” (м. Боярка, Київської області), на ВО ”Хімпром” (м. Вінниця), на Ладижинській ТЕС (м. Ладижин, Вінницької області), ВАТ ”Вінницягаз” та на Вітебському телевізійному заводі ВО ”Вітязь” (м. Вітебськ, Республіка Білорусь). Впровадження підтверджуються відповідними актами.

Використання одержаних у роботі результатів дозволило розробити промисловий обчислювач об'єму газу ”Універсал”, який зареєстрований у Державному реєстрі засобів вимірювальної техніки, допущених до застосування в Україні за №У759-97, що підтверджується сертифікатом затвердження типу засобів вимірювальної техніки, виданим Державним комітетом України по стандартизації, метрології та сертифікації (Держстандартом України) за №UA-М1/1-167-97, виданим 18.11.1997 року та сертифікатом відповідності засобів вимірювальної техніки затвердженому типу, виданим Держстандартом України за №UA-М1/2-119-97 18.11.1997 року. Вищеназваний обчислювач витрат газу дозволено серійно випускати Вінницькому державному технічному університету.

Одержані наукові результати, які використані при розробці систем контролю та обліку всіх видів енергоресурсів та ТЕЦ ПТМ ПЕО ”Вінницяенерго”, Державному концерні ”Арксі” та ВО ”Хімпром”, дозволили підвищити точність та оперативність контролю за виробітком та споживанням основних енергоресурсів на промислових підприємствах (електричної енергії, перегрітої пари, холодної та гарячої води, стислого повітря, тощо), що дає можливість оперативно контролювати споживання енергоресурсів, здійснювати заходи по зменшенню їх витрат.

Впровадження систем контролю та обліку всіх видів енергоресурсів тільки на ТЕЦ ПТМ ПЕО ”Вінницяенерго” та концерні ”Арксі” дало економічний ефект більше 80 тис. рублів у цінах на 1січня 1991 року.

На Ладижинській ТЕЦ впроваджені результати роботи при створенні комп'ютерних систем підготовки операторів, що дозволило досить просто, порівняно з існуючими методами, реалізувати алгоритми оцінки стійкості та діагностики стану систем управління технологічними процесами на даній ТЕЦ та підвищити точність рішення таких задач.

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

Апробація результатів.

Викладені в дисертації результати досліджень були оприлюднені на 46 наукових конференціях, симпозіумах, серед них: ”Гагарінські наукові читання з космонавтики та авіації” (м. Москва, 1985 р.), Всесоюзна науково-технічна конференція (НТК) ”Структурне діагностування пристроїв радіоелектронної апаратури” (м. Москва, 1986 р.), IV Всесоюзна нарада ”Надійність, живучість та безпека автоматизованих комплексів” (м. Москва, 1988 р.), Всесоюзна НТК ”Технічне та програмне забезпечення комплексів пів натурного моделювання” (м. Москва, 1988 р.), ІІІ та IV радянсько-польська НТК ”Комплексна автоматизація у промисловості” (Вроцлав, 1988 р, Київ, 1990 р.), Всесоюзні НТК ”Вимірювальні інформаційні системи” (ВІС-85, 91) (м. Вінниця, 1985 р., м. Санкт - Петербург, 1991р.), ІІ міжнародна НТК ”Актуальні проблеми фундаментальних наук” (м. Москва, 1994 р.), І, ІІ та IV Всеукраїнські конференції з обробки сигналів і зображень та розпізнавання образів (УКРОБРАЗ - 92, 94, 98) (м. Київ, 1992, 1994, 1998 роки), Міжнародній НТК ”Інформаційні технології і системи” (м. Львів, 1993 р.), І-IV міжнародних НТК ”Контроль і управління в технічних системах (КУТС - 92, 93, 95, 97) (м. Вінниця, 1992, 1993, 1995, 1997 роки), Всесвітньому конгресі ITS-93 ”Інформаційні комунікації, мережі, системи і технології” (м. Москва, 1993р.), І-IV Українських конференціях по автоматичному управлінню” (Автоматика - 94, 95, 96, 97) (м. Київ, 1994 р., м. Львів, 1995 р., м. Севастополь, 1996 р., м. Черкаси, 1997 р.), Міжнародний НТК ”Математичне моделювання в електротехніці та електроенергетиці” (м. Львів, 1995 р.), Всеросійській НТК ”Автоматизація досліджень, проектування та випробування складних технічних систем і технологічних процесів” (м. Калуга, 1995 р.), V Українській НТК ”Пристрої перетворення інформації для контролю і управління в енергетиці” (м. Харьків, 1996 р.), Міжнародній конференції по математичному моделюванню (м. Херсон, 1998 р.).

Патент Російської Федерації №2069045 ”Промисловий лічильник витрат газу”, серед авторів якого є пошукувач, демонструвався на Міжнародній виставці винаходів та нових ідей ”Genius-96” (м. Будапешт, Угорщина), на ІІІ Міжнародній виставці винаходів, ідей та ноу-хау ”Iтveтtica-96” (м. Яси, Румунія).

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

Публікації.

Результати дисертації опубліковані у монографії, навчальному посібнику, 14 статтях у наукових журналах, 3 збірниках наукових праць, 40 матеріалах і тезах конференцій, 2 авторських свідоцтвах СРСР, 6 патентах України і трьох патентах Російської Федерації.

Обсяг і структура дисертації.

Дисертаційна робота складається із вступу і семи розділів, списку використаних джерел та трьох додатків, загальний обсяг дисертації 406 сторінок, з яких основний зміст викладений на 272 сторінках друкованого тексту, містить 58 рисунків, 3 таблиці. Список використаних джерел складається з 277 найменувань. Додатки містять виведення проміжних результатів, протоколи випробувань систем і приладів, при розробці яких використані наукові результати, отримані в дисертаційній роботі та акти впровадження результатів роботи.

2. ОСНОВНИЙ ЗМІСТ РОБОТИ

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

У першому розділі проведений огляд методів аналізу поведінки нелінійних систем із змінними параметрами.

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

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

Також показано, що при повільній зміні параметрів нелінійних систем, для оцінки стійкості такого класу систем у частотній області, можливо використовувати Фур'є - інтегральний метод ідентифікації, який був розроблений у 80-х роках професором Б.І. Мокіним.

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

У другому розділі здійснено побудову нового простору опису поведінки нелінійних нестаціонарних систем - простору приростів параметрів.

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

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

Нині задача оцінки асимптотичної стійкості динамічних систем частіше всього розв'язується в просторі змінних стану. Якщо ж первинні параметри нелінійних систем змінюються з часом, то аналітичний розв'язок задач аналізу поведінки таких систем відомими методами в просторі змінних стану отримати практично неможливо, а дослідження процесів в них проводиться, в основному, за допомогою безпосереднього моделювання на ЕОМ.

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

Причому така невизначеність часто виникає за рахунок зміни чисельних значень первинних параметрів любих систем, що обумовлена впливом неконтрольованих параметричних збурень.

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

Припустимо, що мірні вектори керуючих (вхідних) впливів X(t), первинних параметрів систем, що аналізуються, P(t) і координат,що спостерігаються, Z(t) представляються кусково-неперервними функціями, заданими на інтервалі часу [t0, T], тобто:



Природно вважати, що при виконанні умов (1)-(3) мірний вектор змінних стану Y(X(t), P(t), або t) також буде представлятися кусково-неперервною функцією, яка задана на інтервалі часу [t0,T], тобто:

При виконанні умов (1)-(4) функція Y(X(t), P(t),t) здійснює відображення Уk (оператор Немицького) кусково-неперервних функцій, визначених на інтервалі часу [t0,T] наступного вигляду:

Що задається тотожністю:

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

Відзначимо, що якщо:

- належать, відповідно, лінійним нормованим просторам X, P і Y, а U, тоді формально можна записати:



Де збіжність (7) приймається по нормі в просторі Y.

Слід підкреслити, що перша варіація може бути і не лінійна по h. Якщо ж лінійність має місце, тобто якщо:

То відображення Уk має в точці першу варіацію і при цьому існує лінійний неперервний оператор L(P,Y) такий, що:

Підкреслимо також загальновідомий факт, що диференціювання є лінійний оператор. Сукупність L(P, Y) всіх неперервних лінійних операторів, які визначені на всьому просторі P і здійснюють відображення P в Y (де P і Y - нормовані простори), утворять, по відношенню до операцій додавання та множення на число, лінійний простір. Якщо ж, крім того Y - банаховий простір, то й L(P, Y) - банаховий.

Таким чином, якщо мірний вектор змінних стану системи Y(X(t), P(t), t) здійснює деяке відображення Уk у банаховому просторі, то сукупність операторів L(P, Y) буде також банаховим простором. Проте на жаль, простори кусково-неперервних функцій (1)-(4) не є повними. Тому, як правило, розглядають їх вкладення в банаховий простір L - простір вимірювальних суттєво обмежених вектор - функцій.

Після виконання даної операції введемо у розгляд новий простір, що буде банаховим простором в силу наведеного вище, а опис поведінки систем будемо проводити у вигляді:



Де:

P0 - мірний вектор початкових значень первинних параметрів системи вимірністю (p).

Якщо ж вектор змінних стану системи задає неперервну функцію Y(X(t), P(t) або t), визначену на множині:

то при цьому функція Y(X(t), P(t), t) здійснює відображення У неперервних функцій (що є оператором Немицького), визначених на інтервалі часу [t0, T]:

Дане відображення задається тотожністю:

При виконані (11)-(14), (10) можна представити у вигляді:



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

Слід підкреслити, що опис систем у вигляді (10) і (15) проводиться за допомогою апарату теорії чутливості. Використання такого апарату дозволяє перейти від нелінійних нестаціонарних диференціальних рівнянь, що описують поведінку систем у просторі змінних стану, до лінійних нестаціонарних диференціальних рівнянь відносно зміни приростів змінних стану, які обумовлені зміною первинних параметрів системи, що полегшує моделювання поведінки класу систем, що розглядаються.

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

То для систем другого порядку у випадку їх опису у вигляді (15) аналоги фазових траєкторій являють собою залежності:



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

Підкреслимо, що якщо Y(t) визначена на множині Rm, є неперервною та має неперервні частинні похідні першого порядку, то отримання опису систем у просторі приросту параметрів (15) значно спрощується порівняно з її описом у вигляді (10).

Тому, при описі поведінки систем у просторі приростів параметрів бажано передбачити деякі операції, що без суттєвих втрат у точності розв'язання, дозволяли б використовувати опис виду (15), а не (10).

Отримати конструктивні результати при описі поведінки систем у новому просторі можливо шляхом чисельного інтегрування рівнянь чутливості динамічних систем, що аналізуються.

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

У третьому розділі розроблена методологія аналізу поведінки нелінійних нестаціонарних систем (МАП ННС) у просторі приростів параметрів.

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

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

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

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

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

Крім того, для його практичного застосування, як уже відзначалося, необхідно передбачити деякі операції, які б дозволяли використовувати опис систем у вигляді (15), а не (10) в зв'язку із значним спрощенням математичного апарату, який використовується, у цьому випадку, при розв'язку задач стійкості класу систем, що розглядається.

Тому ще однією складовою в МАП ННС повинні входити саме такі операції, які дозволяють використовувати розроблений метод оцінки стійкості в інженерній практиці.

Причому такі операції необхідно передбачити для різних характерів зміни первинних параметрів P(t), які входять у (15) та обумовлені впливом неконтрольованих параметричних збурень.

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

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

Таким чином, МАП ННС, як витікає із вищенаведеного, повинна складатися із взаємопов'язаного ланцюга проблем, які розглядаються як етапи розв'язку задачі оцінки стійкості.

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

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

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

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

Такі вирази отримані при лінійних характерах зміни первинних параметрів типових ланок з однозначними нелінійними характеристиками, а також при гармонічних характерах зміни первинних параметрів типових ланок:

- з однозначними нелінійними статичними характеристиками;

- з неоднозначними нелінійними статичними характеристиками.

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



Табл. - Описуючі функції типових неоднозначних нелінійних ланок при гармонічних характерах зміни їх параметрів:

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

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

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

У четвертому розділі з метою підтвердження можливості використання запропонованої МАП ННС розв'язано ряд задач оцінки стійкості та чутливості нелінійних інформаційно-вимірювальних систем (ІВС) із змінними параметрами.

Так, для електроприводів релейної дії розв'язана задача вибору параметрів релейних елементів, що забезпечують мінімальну чутливість вихідного сигналу до їх змін і визначення областей зміни параметрів, які забезпечують стійкість такого класу електроприводів у введеному просторі приростів параметрів. Скориставшись методом операторів, отримані нестаціонарні рівняння відносно функцій чутливості вихідного сигналу електроприводу до зміни вихідного сигналу B(t) релейної ланки (Ub(t)) та до зміни зони нечутливості C(t) релейної ланки (Uс(t)) у припущенні, що зміни параметрів B(t) та C(t) мають лінійний характер:



Тоді:

Розв'язок рівнянь (17) та (18), отриманий за допомогою пакету MathCad - PLUS 5.0, наведено на рисунках 1 та 2.

Отримані результати дозволяють зробити висновок про слабкий вплив частоти вхідного сигналу на значення відповідних функцій чутливості Ub(t) та Uс(t).

Рис. 1. - Графіки зміни функції чутливості вихідного сигналу електроприводу до зміни вихідного сигналу релейної ланки:

Як слідує з рис. 1 та рис. 2, існують такі параметри, при яких Ub=0 та Uc=0.

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



Рис. 2. - Графіки зміни функції чутливості вихідного сигналу електроприводу до зміни зони нечутливості релейної ланки:

Аналіз рівнянь (17) та (18) дозволив отримати аналітичну залежність, яка задає необхідні умови стійкості системи та визначити області стійкості у новому просторі, представлені на рисунку 3.

Рис. 3. - Область стійкості електропривода в області параметрів релейного елементу:



Використовуючи МАП ННС, для ІВС, які побудовані на основі електромеханічних компенсаційних систем, що працюють у режимі автоколивань, узагальнена структурна схема яких наведена на рисунку 4, отримано опис поведінки системи такого типу у просторі приростів параметрів.

Рис. 4. - Узагальнена структурна схема ІВС, що будуються на основі електромеханічних компенсаційних систем, які працюють у режимі автоколивань:

На рисунку 4:

ЧЕ - чутливий елемент, РС - рухома система, ДП - датчик положення, ПППр - підсилювально-перетворювальний пристрій;

ЕК - електронний ключовий пристрій;

СДЕТ(Н) - стабілізоване джерело електричного струму (напруги);

КП - перетворювач;

РВ - рухомий вузол;

НЛ - нелінійна ланка;

Fвх, Fкп вхідний вплив, відповідно;

u - вихідний сигнал НЛ.

Поведінка системи у новому просторі описується рівняннями:



Для системи, поведінка якої описується рівняннями (19) та (20) визначені достатні умови абсолютної стійкості шляхом вивчення асимптотичної поведінки рішень системи рівнянь відносно канонічних складових розв'язку системи рівнянь вільних коливань. Знаходження таких умов дозволило отримати області стійкості системи, що розглядається. Результати, що отримані, дали можливість побудувати області зміни первинних параметрів, у яких зберігається властивість стійкості, які зображені на рисунках 5 та 6.

Рис. 5. - Область ІВС, за умови, що с>0:

Рис. 6. - Область ІВС, за умови, що с<0:

Шляхом математичного моделювання за допомогою пакету прикладних програм МОДО-С, розробленому у Московському енергетичному інституті, підтверджена адекватність отриманих результатів. Для отримання вимірювальної інформації в обробляючих центрах автоматизованих виробництв часто входять ІВС, які являють собою головку дотику типу ”індикатор контакту”, структурна схема яких наведена на рис. 7.

Рис. 7. - Структурна схема ІВС обробляючих центрів типу ”індикатор контакту”:

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

Використовуючи МАП ННС, отримані рівняння динаміки ІВС у просторі приростів параметрів, які мають вигляд:



Для чисельних даних, які характеризують обробляючий центр ФП27С, отриманий розв'язок рівнянь (21), (22) за допомогою пакету MathCad PLUS 5.0, наведений на рисунках 8 та 9.

Рис. 8. - Графік функції чутливості до змін параметру b:

Рис. 9. - Графік функції чутливості до змін параметру c:

Як слідує з рис. 8, 9 вихідний сигнал ІВС обробляючих центрів типу ”індикатор контакту” більш чутливий до зміни параметру с, що характеризує зміну зони нечутливості релейного елементу, ніж до зміни параметру b. Тому елементи, що забезпечують задану зону нечутливості. Крім того, існують такі параметри релейних елементів, при яких Ub=0 та Uc=0, тобто прирости відповідних параметрів у деяких областях не впливають (або мало впливають) на вихідний сигнал системи, що розглядається.

Такі особливості поведінки функцій чутливості дадуть можливість проектувати ІВС, інваріантні до зміни параметрів релейних ланок, що входять в їх структурні схеми.

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

У п'ятому розділі, застосовуючи МАП ННС, розв'язана задача стійкості та оптимізації параметрів САР силових трансформаторів з пристроями регулювання під навантаженням (РПН), які є головним і найбільш ефективним засобом забезпечення якості напруги в електричній мережі. Побудована структурна схема САР силових трансформаторів з РПН, ланками якої є суттєво нелінійні, логічні елементи та елементи цифрової техніки. Отримана структурна схема дозволяє здійснити комп'ютерний аналіз стійкості таких САР за допомогою пакету програм МОДО-С. Моделювання САР, що розглядається, дозволило отримати перехідні функції для робочого діапазону часу затримок. Такі функції мають аперіодичний характер та дозволяють зробити висновок по стійкість таких САР при зміні її відповідного параметру в усьому робочому діапазоні. Аперіодичний характер перехідних процесів обумовлюється фільтруючими властивостями електроприводу пристрою РПН, за час переключення якого встигають згаснути всі процеси в інших елементах схеми. Головними параметрами регуляторів рівня напруги, які можливо настроювати, є зона нечутливості та час затримки. Якщо використати можливості ППП МОДО-С, то крім задачі оцінки стійкості САР, можливо розв'язати і більш складну задачу - задачу оптимізації таких її параметрів.

В якості критерію оптимізації використовується вираз:

Де:

(23) С1 - вартість збитків підприємств електричних мереж, які обумовлені відмовою трансформаторів;

С2 - вартість збитків споживача, які обумовлені поганою якістю електричної енергії (по напрузі) на його шинах;

J1 - функціонал, який визначає ймовірність відмови трансформатору з РПН при фіксованій зоні нечутливості;

J2 - функціонал, який визначає зниження якості електричної енергії (по напрузі) при такій самій зоні нечутливості.

Функціонал J1 та функціонал J2 запропоновано підрахувати за допомогою виразів:

У виразах (24) та (25):

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

f - умовна густина ймовірності відмови трансформатора, яка є функцією часової затримки при фіксованому значенні зони нечутливості регулятора;

t1-t0 - інтервал спостереження;

Kj - ”вага” споживача електричної енергії;

Ujp(t) - біжуче значення напруги на шинах j-того споживача, яке залежить від параметрів, які настроюються у регуляторів силових трансформаторів з пристроями РПН;

Ujном - номінальне значення напруги на шинах у j-го споживача;

N - кількість споживачів електроенергії, які підключені до вузла мережі, що розглядається.

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

Зрозуміло, що пошук оптимальних за критерієм (23) значень часової затримки регулятора будемо проводити саме в цьому діапазоні.

З досвіду статистичної оцінки умовної густини ймовірностей f відмови трансформаторів з РПН, яка входить у складову J1 (24) критерію (23) встановлено, що її можливо апроксимувати гамма - розподілом:

Таким чином, якщо використати умову (28), а відповідні значення частотних характеристик знайти за допомогою ФІМІ, то можна побудувати пристрій для автоматичного визначення стійкості замкнених нелінійних систем, які можливо, по частотних характеристиках розімкнених систем. При цьому підкреслимо, що якщо організувати по такту роботу пристрою, що пропонується, таким чином, щоб за такт його роботи параметри практично не змінювались, то такий пристрій.



ВИСНОВКИ

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

В дисертаційній роботі одержані такі основні наукові та практичні результати:

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

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

3. На основі отриманих результатів побудована методологія аналізу поведінки нелінійних нестаціонарних систем (МАП ННС), придатна для розв'язання широкого класу наукових та інженерних задач;

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

5. Шляхом математичного моделювання, на основі МАП ННС, розв'язана задача стійкості та запропонована нова методика параметричної оптимізації одного класу дискретних нелінійних нестаціонарних систем автоматичного регулювання в електроенергетиці - САР силових трансформаторів з пристроями регулювання під напругою, які є головним і найбільш ефективним засобом забезпечення якості напруги в електричній мережі;

6. На основі Фур'є - інтегрального методу ідентифікації (ФІМІ), розробленого професором Б.І. Мокіним, запропоновані алгоритми оцінки стійкості, діагностики стану нелінійних нестаціонарних систем. Дані алгоритми дозволили розробити функціональні схеми пристроїв, які апаратним шляхом розв'язують задачі оцінки стійкості та діагностування такого класу систем у частотні області. Розроблені пристрої захищені авторськими свідоцтвами на винаходи СРСР, патентом України та патентом Російської Федерації;

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

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



СПИСОК ОПУБЛІКОВАНИХ ПРАЦЬ

1. Юхимчук С.В. Математические модели оценки устойчивости нелинейных нестационарных систем. - Винница.: Универсум - Винница, 1997. - 141 с.

2. Мокин Б.И., Юхимчук С.В., Головатюк Н.А. Об одном способе определения оптимальных параметров настройки регуляторов напряжения силовых трансформаторов с устройствами РПН // Известия вузов. Сер. Энергетика. - 1990. - №11. - С. 26-29.

3. Мокин Б.И., Юхимчук С.В. К вопросу применения теории чувствительности для исследования поведения нелинейных САУ // Известия вузов. Серия Электромеханика. -1991. - №10. - С. 91-96.

4. Мокин Б.И., Юхимчук С.В. Определение устойчивости нелинейной САУ с переменными параметрами в пространстве приращений параметров // Автоматика. - 1993 - №4. - С. 91-94.

5. Мокин Б.И., Юхимчук С.В. Исследование влияния изменения параметров релейного звена на выходной сигнал в информационно-измерительных системах // Электронная техника. Сер.10. Микроэлектронные устройства. - 1993. - Вып. 3, 4 (97, 98). - С. 31-33.

6. Дементьев Ю.В., Мокин Б.И., Юхимчук С.В. Контроллер с расширенными функциональными возможностями по входу и выходу // Управляющие системы и машины. - №1/2. - 1994. - С. 14-16.

7. Юхимчук С.В., Батиха Х.М. Обобщение метода гармонической линеаризации для исследования нелинейных электромеханических систем с переменными параметрами // Известия вузов. Сер. Электромеханика. - 1994. - №1-2. - С. 9-15.

8. Юхимчук С.В., Батіха Х.М. Визначення області стійкості прецизійних давачів механічних ве-личин, які працюють в режимі автоколивань // Вісник Вінницького політехнічного інституту. - 1994. - №3(4). - С. 7-11.

9. Юхимчук С.В., Головченко О.М., Каплун В.А., Юношев І.М., Коханчук Ф.П. Питання методичного та програмного забезпечення комп'ютеризації навчання робітничим професіям операторів енергетичних виробництв // Вісник Вінницького політехнічного інституту. - 1996. - №1(10). - С. 50-53.

10. Юхимчук С.В., Кондратенко Н.Р. Метод розв'язку лінійних диференційних рівнянь вищих порядків за допомогою елементів комп'ютерної алгебри // Вісник Вінницького політехнічного інституту. - 1997. - №1. - С. 83-87.

11. Юхимчук С.В. Спосіб розв'язку задачі стійкості нелінійних нестаціонарних систем // Вісник Житомирського інженерно-технологічного інституту. - 1997. - №6. - С. 126-129.

12. Юхимчук С.В. Узагальнення методу описуючих функцій для неоднозначних нестаціонарних систем // Вісник Вінницького політехнічного інституту. - 1997. - №4. - С. 5-9.

13. Видмиш А.А., Юхимчук С.В. Алгоритмічна корекція похибки первинних давачів // Вісник Вінницького політехнічного інституту. - 1997. - №4.

14. Юхимчук С.В. Визначення областей стійкості систем керування релейним електроприводом // Праці 4-ої Української конференції з автоматичного управління (Автоматика - 97). - Том ІІІ. - Черкаси, ЧІТІ. - 1997. - С. 141.

15. Юхимчук С.В. Методологія аналізу поведінки нелінійних нестаціонарних систем у просторі приростів параметрів // Вісник ВПІ. - 1998. - №4.

16. Юхимчук С.В. Математические модели оценки робастной устойчивости нелинейных нестационарных систем // Нелинейные краевые задачи математической физики и их приложения: Сб. научн. тр. -К.: Институт математики НАНУ. - 1998. - С. 242-244.

17. Юхимчук С.В., Ваховська Л.М., Запорожан А.В. Визначення областей стійкості релейного електроприводу стежної дії зі змінними параметрами // Вісник Вінницького політехнічного інституту. - 1998. - №4. - С. 53-56.

18. Юхимчук С.В., Лысогор В.Н., Марущак В.Ю. Автоматизация проектирования систем автоматического управления: Учебн. пособие. - К.: УМК ВО, 1989. - 172 с.

19. Головатюк Н.А., Ткачук Б.Д., Юхимчук С.В. Моделирование систем автоматического регулирования в электрических сетях // Контроль и управление в энергетике. - Киев: УМК ВО, 1988. - С. 54-62.

20. Мокин Б.И., Юхимчук С.В. Математические модели процессов регулирования технологических объектов электрической сети // Труды ІІІ конф. ”Компле-ксная автоматизация промышленности” - Вроцлав. - 1988.

21. Дементьев Ю.В., Юхимчук С.В. Микропроцессорные системы контроля технологических параметров в теплоэнергетике //Комплексная автоматизация промышленности. - Труды четвертой международной научно-техн. конф. Секция 5. - К.: 1990. - С. 21-25.

22. Кузьмин И.В., Юхимчук С.В., Мусаев И.М. Подход к определению времени контроля и управления сложных адаптивных систем //Адаптивные системы автоматического управления. - 1987. - №15. - С. 65-68.

23. Юхимчук С.В. Метод і пристрій для діагностики лінійних динамічних об'єктів // Праці першої Всеукраїнської конференції ”Обробка сигналів і зображень, розпізнавання образів”(УкрОБРАЗ-92). - Київ: ІК ім. В.М. Глушкова. - 1992. - С. 211-214.

24. Юхимчук С.В. Оценка стабильности динамических характеристик ИИС обрабатывающих центров // Контроль и управление в технических системах. - Книга по матер. IV международной научно-техн. конф. Том 3. - Винница. - 1997. - С. 244-250. моделювання нелінійний математичний

25. Юхимчук С.В. Метод и аппаратура контроля стабильности динамических характеристик обрабатывающих центров // Современная контрольно-испытательная техника промышленных из-делий и их сертификация. - Труды первой научно-практической конф. Том2. - Киев. - 1997. - С. 158-161.

26. Позитивне рішення по заявці №93006169/2(1375) на видачу патенту України від 27.04.98 р., МПК G 05 В 23/02. Устройство для определения устойчивости линейных систем управления / С.В. Юхимчук, Х.М. Батиха, Заявлено 16.07.93. - 13 с.

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

...