Главная Контакты Добавить в избранное Авторы Вопросы и ответы
,

УДК 681.511.4: 681.5.037

МОДЕЛИРОВАНИЕ  ДИНАМИКИ ВЫХОДА ТЕХНОЛОГИЧЕСКОГО АППАРАТА ИЗ КРИТИЧЕСКОГО РЕЖИМА ЭКСПЛУАТАЦИИ

Завальнюк И.П., Бражник А.М., Завальнюк О.П.

Введение. Современное производство характеризуются напряженным режимом работы оборудования, что обусловлено требованиями к эффективности технологических процессов. При этом увеличивается  вероятность переход системы в критический режим, в котором нагрузка выше номинальной [1, 2, 3], с чем связано возникновение аварийных ситуаций. В критическом режиме контур диссипации системы перегружен и теряет способность рассеивать энергию во внешнюю среду, т.е. коэффициент диссипации падает и появляется возможность накопления энергии. При неправильном управлении (естественной реакцией в такой ситуации является резкое снятие управления) энергия высвобождается в виде динамического удара в контуре диссипации создавая угрозу аварии. 

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

Анализ последних исследований и публикаций. Вопросы, связанные с особыми режимами нелинейных систем широко освещены в теории катастроф [4, 5, 6], однако в данном случае система не изменяет атрактора и не теряет устойчивости, так как при снижении диссипации просто происходит накопление энергии в системе и резко увеличивается колебательность. Так в работе [7] рассмотрено, что наряду с возмущающей силой на систему в процессе вынужденных колебаний действуют и внутренние силы колебательной системы – силы инерции, силы демпфирования и восстанавливающие силы связанные с рассеиванием энергии - диссепацией. На основе сравнения энергии, поступающей в систему за счет внешних сил  с энергией, расходуемой на преодоление демпфирования  с использованием  энергетической диаграммы показано, что если работа демпфирования меньше работы системы , то в систему посту­пает больше энергии, нежели расходуется на преодоление сил демп­фирования, и, следовательно, амплитуда увеличивается. При  потребляемая демпфированием энергия больше энергии, производимой внешней силой, и вследствие этого происходит умень­шение амплитуды. При  обе работы равны между собой: , где  - частота колебаний.

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

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

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

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

 

.

(1)

где  ,  - функция диссипации энергии в системе.

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

 

,

(2)

 

где энергию системы оценим по фазовым переменным .

Исследуем ситуацию, когда исходное значение диссипации исключает колебательность (), а при потере диссипации, её остаточное значение очень мало (). Подчеркивая потерю эффективности контура диссипации  примем значения коэффициентов c и d в уравнении (3) равными 8 и 0.1 (соответственно).

Движение системы (1) при отсутствии управления описывается уравнением [8]:

 

(3)

где , ,  - постоянная (начальная) фаза.

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

 

Рис.1 Вид функции диссипации энергии в системе при с = 8, d = 0.1

 

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

 

(4)

 

Амплитуда управления  обеспечивает выход системы в область малой диссипации, а постоянная времени T – больше, чем собственная потоянная времени объекта. Таким образом, выходной сигнал системы:

 

(5)

 

После завершения переходного процесса система находится в критическом состоянии и благодаря плавному разгону и выполнению условия d>0:

 

(6)

 

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

 

(7)

 

Проведем анализ поведения системы при резком снятии нагрузки – естественной реакции на возникновение критического режима. Тогда управление изменяется следующим образом:

 

(8)

где 1(t) – ступенчатая функция.

В структурной схеме моделирования (рис.2) в стандартную модель объекта 2-го порядка введен контур диссипации (Product1, Product2, Product3, Sum2, Fcn, блок а22); чтобы организовать требуемое управление (4) использованы блоки Constant1 и Transfer Fcn.

 

Рис.2 Структурная схема моделирования системы при резком снятии нагрузки

 

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

Мощность удара определяется изменением энергии в системе ΔE за время Δt* прохождения участка с изменением диссипации:

 

,

(9)

 

где ,  - запас энергии в системе до и после снятия нагрузки

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

 

(10)

 

На рис.3 приведена структурная схема моделирования системы с плавным снятием управления, для чего введены блоки: Integrator1, Saturation и блок b. Важно, что длительность перехода системы в начальное состояние, длительность рассеяния энергии можно регулировать, изменяя значение коэффициента b.

 

Рис. 3 Структурная схема моделирования системы при постепенном снятии нагрузки

 

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

а)

б)

Рис. 4 Осциллограммы сигналов при резком (а) и постепенном (б) снятии нагрузки: 1,2 – сигналы  и ; 3 – оценка количества энергии в системе; 4 – характер диссипации энергии; 5 - изменение энергии после диссипирующей цепи.

 

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

а)

б)

Рис. 5 Фазовые портреты системы.

 

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

 

Выводы:

                  1.     Моделирование подтверждает гипотезу о возможности возникновения ударных нагрузок в контуре диссипации при резком снятии управления, что существенно для систем с ограниченными ресурсами контура диссипации.

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

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

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

 

The behaviour of the dynamic system is сonsidered with small dissipation of the energy under different type of the removing the load. It is installed that in system can occur the powerful impact because of cutting the removing of control. The certain law of the change of control, guaranteeing absence of the striking loads in system when leaving her from critical mode.

 

                  1. Андронов А.А, Витт А.А, Хайкин С.Э. Теория колебаний.- 3-е изд., испр.- М.: Наука. Гл. ред. физ.-мат. лит., 1989. – 916 с.

                  2. Анищенко В.С. Сложные колебания в простых системах: Механизмы возникновения, структура и свойства динамического хаоса в радиофизических системах.- М.: Наука. Гл. ред. физ.-мат. лит., 1990. – 312 с.

                  3. Джакалья Г.Е.О. Методы теории возмущений для нелинейных систем: Пер с англ./ Под ред. А.П. Маркеева.- М.: Наука. Гл. ред. физ.-мат. лит., 1988. – 320 с.

                  4. Дж.М.Т.Томпсон Неустойчивость и катастрофы в науке и технике. Перевод с английского В.Л. Бердичевского. –М.: Мир., 1985.-254С.ил.

                  5. Тим Пост, Иэн Стюарт Теория катастроф и ее приложения. Перевод с английского А.В. Чернавского. – М.: Мир. 1980. 606С. Ил.

                  6. В.И. Арнольд Теория катастроф. -3-е изд.,доп.-М.: Наука гл.ред.физ.-мат. лит., 1990.-С.128.

                  7. Магнус К. Колебания: Введение в исследование колебательных систем: Пер. с нем. – М.: Мир, 1982. – 304 с., ил.

                  8. Теория систем автоматического регулирования, издание третье, исправленное / Бесекерский В. А., Попов Е. П. – М.: Наука, Гл. ред. физ.- мат. лит-ры, 1975. – 768 с.

 





Ответы на вопросы [_Задать вопроос_]

Моделирование объектов и систем управления

Соколов А.Е., Махова Е.О. Моделирование процесса принятия педагогического решения при компьютеризированном обучении

Славко О.Г. Порівняльний аналіз керування регулятором на основі локальної моделі керованого процесу та П-регулятором

Войтенко В.В., Дикусар Е.В, Ситников В.С. Определение частоты среза устройства сглаживания данных на основе метода скользящего среднего

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

Ляшенко С.А, Ляшенко А.С. Оценка модели псевдолинейной регрессии

Ладієва Л.Р. Математична модель процесу газової мембранної дистиляції

Носов П.С., Косенко Ю.І. Нечіткі моделі і методи ідентифікації та прогнозу стану інформаційної моделі студента

Китаев А.В., Глухова В.И. Анализ работы синхронного двигателя с неявнополюсным ротором по данным каталога

Дорошкевич В.К., Пироженко А.В., Хитько А.В., Хорольский П.Г. К определению требований к системам увода космических объектов

Голінко І.М., Ковриго Ю.М., Кубрак А.І. Настройка системи керування за імпульсною характеристикою об’єкта

Яшина К.В., Садовой А.В. Комплексная математическая модель тепловых процессов, происходящих в дуговых электросталеплавильных печах

Шейник С.П., Рудакова А.В. Использование функций принадлежности для моделирования параметров распределенных объектов

Хомченко А.Н., Литвиненко Е.И. Метод барицентрического усреднения граничных потенциалов электростатического поля

Селяков Е. Б. Моделирование требований к техническим системам методами математической логики

Тодорцев Ю.К., Ларіонова О.С., Бундюк А.М. Математична модель контура теплопостачання когенераційної енергетичної установки

Кириллов О.Л. , Якимчук Г.С. Моделирование процесса управления системой перегрузки углеводородных жидких топлив

Шеховцов А.Н., Козел В.Н. Построение математической модели формирования распределенных систем

Китаев А.В., Глухова В.И. Анализ поведения генератора постоянного тока по данным каталога

Хомченко А.Н., Козуб Н.О. Задачі наближення функцій: від лагранжевих до серендипових поліномів

Хобин В.А., Титлова О.А. Определение температуры парожидкостной смеси в дефлегматоре АДХМ по результатам измерений температуры его поверхности

Григорова Т.М., Усов А.В. Вероятностно-статистическое моделирование маршрутизированных пассажиропотоков в крупных городах

Горач О.О., Тернова Т.І. Моделювання технологічного процесу одержання трести при використані штучного зволоження з урахуванням складу мікрофлори

Дубік Р.М., Ладієва Л.Р. Математична модель розділення неоднорідних рідких систем

Казак В.М, Лейва Каналес Родриго, Яковицкая Е.Ю. Моделирование динамики полета магистрального самолета на исследовательском стенде

Завальнюк И.П. Исследование процесса торможения автомобиля как критического режима динамической системы

Дмитриев С.А., Попов А.В. Построение портрета неисправностей проточной части газотурбинного двигателя на примере АИ-25

Русанов С.А., Луняка К.В., Клюєв О.І., Глухов Г.М. Математичне моделювання робочого процесу в апаратах з віброкиплячим шаром та розробка систем автоматизованого моделювання гідродинаміки віброкиплячих шарів

Боярчук В.П., Сыс В.Б. Экспериментальные исследования влияния технологии шлихтования на изменение жесткости текстильных нитей

Селін Ю.М. Використовування контекстних марківських моделей для аналізу дії промислових вибухів на будівельні конструкції

Рудакова А.В. Проблемы интеграции сложных систем

Передерій В.І., Касап А.М. Математична модель та алгоритм автоматизації розрахунку параметрів комп’ютеризованих систем працюючих у реальному часі

Передерий В.И., Еременко А.П. Математические модели и алгоритмы принятия релевантных решений пользователями автоматизированных систем с учетом личностных и внешних факторов на базе генетических алгоритмов

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

Хомченко А.Н., Колесникова Н.В. Явление «сверхсходимости» в задаче Прандтля для уравнения Пуассона

Китаев А.В., Глухова В.И. Анализ работы трансформатора по данным каталога

Квасницкий В.В., Ермолаев Г.В., Матвиенко М. В., Бугаенко Б.В., Квасницкий В.Ф. Оценка применимости метода компьютерного моделирования к исследованию напряженно-деформиррованного состояния цилиндрических узлов

Китаев А.И., Глухова В.И. Анализ работы асинхронного двигателя по данным каталога

Шелестов А.Ю Имитационная модель взаимодействия GRID-узлов с очередью доступа к общей памяти

Chizhenkova R.A. Mathematical Aspects of Bibliometrical Analysis of Neurophysiological Investigations of Action of Non-ionized Radiation (Medline-Internet)

Хомченко А.Н., Козуб Н.А. Геометрическое моделирование дискретных элементов с криволинейными границами

Славич В.П. Модель автоматизованої системи управління потоками транспортних засобів

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

Степанкова Г.А., Баклан І.В. Побудова гібридних моделей на основі прихованих марківських моделей та нейронних мереж

Бакшанська Т.Д., Рижиков Ю.Г., Тодорцев Ю.К. Математична модель процесу горіння природного газу з рециркуляцією продуктів згорання для цілей управління

Хомченко А.Н. Новые решения обобщенной задачи Бюффона

Передерий В.И., Еременко А.П. Математические модели и алгоритмы определения релевантности принимаемых решений с учетом психофункциональных характеристик пользователей при управлении автоматизированными динамическими системами

Ложечников В.Ф., Михайленко В.С., Максименко И.Н. Аналитическая много режимная математическая модель динамики газовоздушного тракта барабанного котла средней мощности

Ковриго Ю.М., Фоменко Б.В., Полищук И.А. Математическое моделирование систем автоматического регулирования с учетом ограничений на управление в пакете Matlab

Исаев Е.А., Наговский Д.А. Математическое описание влияния кривизны контактирующих тел на угол смачивания жидкости в межчастичном пространстве

Бідюк П.І., Литвиненко В.І., Кроптя А.В. Аналіз ефективності функціонування мережі Байєса

Тищенко И.А., Лубяный В.З. Математическое моделирование вокодера для определения оптимальной формы импульса сигнала возбуждения.

Николаенко Ю.И., Моисеенко С.В. Моделирование гармонического полиномиального базиса гексагона.

Козуб Н.А., Манойленко Е.С., Хомченко А.Н. Температурный тест для модифицированных базисов бикубической интерполяции.

Клименко А.К. Об упрощенном численном конструировании обратной модели динамического объекта.

Китаев А.В., Сушич Е.Ф. Расчет погрешностей измерительных трансформаторов.

Передерій В.І.,Касап А.М. Математична модель та алгоритм автоматизації розрахунку параметрів комп’ютеризованих систем працюючих у реальному часі

Шпильовий Л.В. Математична модель та алгоритм екстремального управління процесом осадження дисперсної фази суспензії.

Тулученко Г.Я. Інформаційний модуль експрес-пошуку точок еквівалентності процесу нейтралізації.

Тернова Т.І. Урахування морфогенетичного рівняння в математичній моделі тканини.

Попруга А.Г. Теоретические и экспериментальные исследования электрических нагревателей по критерию экономии энергии.