УПРАВЛЕНИЕ ПЕРЕХОДНЫМИ ПРОЦЕССАМИ В СИСТЕМАХ С ОГРАНИЧЕНИЕМ КООРДИНАТ | Вестник Томского государственного университета. Математика и механика. 2008. № 3 (4).

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

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

Transitional Processes Control in the System with Limitation of .pdf Синтез оптимальных по быстродействию управлений линейными и нелинейными системами с ограничением координат традиционными методами не всегда эффективен. Следует применять новые, нестандартные подходы, в которых требуется не столько искусство математика, сколько хорошее знание рассматриваемой технической задачи и понимание, какими факторами можно пренебречь и к каким последствиям это приведет.Для технических систем характерна иерархия вышестоящего и нижестоящего уровней. Продвижение к глобальной цели, стоящей перед всей системой, обычно осуществляется за счет соответствующей координации деятельности подсистем. Глобальная цель развертывается в подцели, причем часто лишь после достижения подцели появляется возможность оценить целесообразность принятия того или иного закона управления. Такие действия следует выполнять путем предварительного прогнозирования поведения системы.Для синтеза оптимальных управлений применим метод избыточных переменных. Следует использовать принципы «перемены цели» и «ведущего слабого звена». Оптимизация управления переходными процессами в реальном масштабе времени микропроцессорными средствами требует разработки простых алгоритмов.Методика последовательного многошагового синтеза оптимальных управленийСистемный анализ позволяет создать методологию исследования и проектирования сложных систем. Во всяком действии легко увидеть его составные части, более мелкие действия. Они должны выполняться не в произвольном порядке, а в определенной последовательности. В сложных системах наблюдается иерархическая система противоречий, которую следует использовать для создания алгоритмов управления. Формы проявления единства и борьбы противоположностей неразрывно связаны со спецификой самих систем. Когда цели связаны с различием между состояниями, приходится отображать процессы изменения состояния с помощью динамических моделей. В настоящее время используется для анализа структур сложных процессов и систем имитационное моделирование, которое позволяет постичь суть явления, не прибегая к экспериментам на реальном объекте. Обычно модель будущего состояния системы выполняется с учетом того, что происходит на промежуточных этапах. Приходится оценивать текущий результат предыдущих действий и выбирать следующий шаг из числа возможных. Необходимо сравнивать последствия возможных шагов, не выполняя их реально.Системный анализ позволяет реализовать диалектический метод при рассмотрении прикладных задач. Он применим для синтеза оптимальных по быстродействию управлений как линейными, так и нелинейными системами с ограничением координат. Новые методики синтеза оптимальных управлений могут быть созданы на основе методов динамического программирования [1] и имитационного моделирования, принципов «перемены цели» и «ведущего слабого звена» [2]. На моделях возможен поиск оптимальных управлений путем исследования переходных процессов в системе, получаемых в результате приложения пробных воздействий.Принцип «перемены цели» служит в качестве средства приспособления системы к изменению параметров, фазовых координат, ограничений и требований к переходным процессам. Принцип «ведущего слабого звена» подразумевает объединение слабых и сильных звеньев для достижения цели, причем сильные звенья, имеющие большие адаптационные возможности, подстраиваются под слабые звенья. В течение переходного процесса главными на разных этапах становятся разные ограничения и требования, причем некоторые ограничения могут и не стать главными для конкретного процесса и не участвовать в формировании оптимального управления. Все зависит от величины задания по главной координате и начального состояния системы. Наличие в любой момент времени переходного процесса одного «слабого звена» существенно упрощает синтез управления.Разработана методика синтеза оптимальных по быстродействию управлений линейными и нелинейными системами путем многократного численного решения обыкновенных дифференциальных уравнений. Оптимальный закон управления системой составляется из управлений, найденных во время переходного процесса для малых интервалов. Расчет этих управлений осуществляется с учетом технологических требований, ограничений фазовых координат, принятого критерия оптимальности, заданного конечного состояния системы и начальных состояний системы на малых интервалах. Методика предполагает действия над разностными уравнениями, которыми описывается поведение системы. В общем случае она представляет собой методику составления программы для численного решения задачи на вычислительных машинах.Поиск оптимальных управлений на малых интервалах ведется последовательно по шагам с учетом координат системы, полученных при оптимальном управлении на предыдущих шагах. Синтез оптимального управления системой сводится к последовательной оптимизации более простых процессов управления. Характерно, что каждый последующий участок управления усложняется, но формируется всегда путем определенной последовательности включения более простых участков управления.На первом этапе методом динамического программирования с учетом принятых ограничений рассчитывается прогнозируемое оптимальное управление для очередного шага. Это управление в дальнейшем может быть скорректировано после проведения проверок на отсутствие нарушений ограничений координат во время переходного процесса. На втором этапе определяются координаты системы в результате выполнения пробного шага с найденным прогнозируемым оптималь-ным управлением. Расчеты ведутся последовательно от входа к выходу системы. На третьем этапе методом имитационного моделирования выполняется перевод системы по оптимальному закону с учетом принятых ограничений из состояния, полученного в результате выполнения пробного шага, в равновесное состояние. Под равновесным состоянием понимается состояние системы, в котором она может оставаться длительное время без изменения координат. На четвертом этапе сравниваются значения координат системы при переводе ее по оптимальному закону в равновесное состояние с допустимыми значениями координат. Если нет нарушений принятых ограничений, то использованное на пробном шаге управление считается оптимальным и его можно использовать для расчета реальных координат системы на очередном шаге интегрирования. Эти координаты используются в качестве начальных условий для определения оптимального управления на следующем шаге. Если наблюдаются нарушения принятых ограничений, то использованное на пробном шаге управление не является оптимальным, его следует скорректировать и повторить расчеты по описанному циклу, начиная со второго этапа.Перевод системы в равновесное состояние выполняется методом имитационного моделирования путем изменения в иерархической последовательности всех координат до установившихся значений. Под установившимся значением фазовой координаты подразумевается такое ее значение, которого она может достичь при изменении по оптимальным законам всех предшествующих фазовых координат до уровней, соответствующих прекращению дальнейшего изменения анализируемой координаты. При изменении координат до установившихся значений по оптимальным законам могут формироваться различные цели, отличные от использованной при расчете пробного шага цели, однако всегда используется принцип «ведущего слабого звена» и идет подстройка под самое «сильное» в данный момент ограничение. Сложность состоит в необходимости одновременного выхода на установившееся значение как анализируемой координаты, так и всех предшествующих координат. Задача усложняется с повышением порядка системы, приходится использовать большее количество пробных шагов, формируя различные цели. Все расчеты выполняются по циклическим алгоритмам, для одной и той же координаты могут формироваться цели выхода по оптимальным законам на разные заданные установившиеся значения. Изменение целей связано с необходимостью выполнения предъявляемых к системе противоречивых требований.Особенностью предложенной методики, в отличие от других методов решения многошаговых задач, в которых приходится анализировать на каждом шаге все возможные варианты управления, является использование промежуточных критериев, позволяющих сразу отсечь заведомо неприемлемые управления и тем самым сократить объем вычислений. Вычисления выполняются по циклическим алгоритмам, обеспечивающим перевод координат в установившиеся состояния, поэтому для решения частных задач в процессе поиска оптимального управления можно использовать подпрограммы.В некоторых случаях удается получить простые аналитические выражения для расчета процесса перевода координат системы в установившиеся состояния после выполнения пробного шага, что открывает широкие перспективы по разработке алгоритмов синтеза в реальном масштабе времени микропроцессорными средствами оптимальных управлений системами высоких порядков. Методика позволяет синтезировать оптимальные управления при переводе объекта в любое требуемое состояние. При этом изменяется только понятие «установившегося значения фазовой координаты» и промежуточные критерии.Оптимальное управление системой третьего порядкаУправляемый линейный объект описывается системой уравнений\ dXi = X XM) следует изменить прогнозируемое оптимальное управление U(n+1) и повторить расчеты по описанному циклу, начиная со второго этапа. Синтез оптимального управления для одного шага выполняется по методике, которая применима и для других шагов.Определение оптимального управления для (n + 1) шага интегрирования выполняется в следующей последовательности. Методом динамического программирования рассчитывается управление, обеспечивающее максимальное приращение координаты Х1 на (n+1) шаге интегрирования.Определяется требуемое приращение по координате Х1АХ1(п+1) = XM - X1(n),(5)затем вычисляется значение координаты Х2, способное обеспечить это приращениеX2(n+1) = AX1(n+1) / At.(6)Определяется требуемое приращение по координате Х2:AX2(n+1) = X2(n+1) - X2(n) (7)и значение координаты Х3, способное обеспечить это приращение:X3(n+1)= AX2(n+1) / At.(8)Вычисляется требуемое приращение по координате Х3:AX3(n+1) = X3(n+1) - X3(n) (9)и обеспечивающее его управлениеU(„+i) = AX3(„+i) / At.(10)Это управление ограничивается, при необходимости, на уровне UM с соответствующим знаком. Рассчитываются координаты объекта после выполнения первого пробного шага (X11(„+1), X21(„+1), X31(„+1)) с найденным управлениемX31(n+1) = X3(n) + U(n+1) At;(1) X21(n+1) = X2(n) + X31(n+1) At; (12) X11(n+1) = X1(n) + X21(n+1) At. (13)Они используются в качестве начальных условий для перевода объекта в равновесное состояние. Это необходимо сделать для проверки выполнения требования к координате Х1, которая не должна превышать значения Хм.Перевод объекта в равновесное состояние начинается с расчета второго пробного шага (n+2), выполняемого с целью скорейшего достижения координатой Х2 значения Х2 = 0. В качестве начальных условий используются координаты системы X11(„+1), X21(„+1), X31(„+1) , полученные в результате расчета первого пробного шага. Расчет ведется методом динамического программирования, аналогично первому пробному шагу, только теперь изменена цель управления. По аналогичной методике определяются координаты системы в результате выполнения второго пробного шага Х12(и+2), Х22(и+2), Х32(и+2) , В качестве «ведущего слабого звена» вновь выступает ограниченное управление.Новые координаты объекта Х32(и+2), Х22(и+2), Х12(и+2) используются в качестве начальных условий в циклическом алгоритме, обеспечивающем изменение координаты X3 до значения Х3 = 0, соответствующего установившемуся значению координаты Х2. Расчеты выполняются методом динамического программирования.После достижения координатой Х3 значения Х3 = 0 оценивается значение координаты Х2. Если оно оказывается больше Х2(Т), то рассчитывается новый пробный шаг по скорейшему достижению координатой Х2 значения Х2(Т), только в качестве начальных условий используются координаты объекта Х32(и+2), Х22(и+2), Х12(и+2) , полученные в результате выполнения предыдущего второго пробного шага. Вновь используется циклический алгоритм, обеспечивающий изменение координаты Х3 до значения Х3 = 0, и оценивается значение координаты Х2.Расчеты по такому циклу продолжаются до тех пор, пока координата Х2 не достигнет значения Х2 = 0. Иногда для одновременного достижения координатами значений Х2 = 0 и Х3 = 0 приходится использовать метод последовательных приближений. Полученные значения Х2 = 0 и Х3 = 0 соответствуют установившемуся состоянию объекта, оценивается значение координаты Х1. Если оно не превышает значение Хм, то найденное на первом пробном шаге управление U(„+1) считается оптимальным. В противном случае управление на (n+1) шаге определяется методом последовательных приближений из диапазона (- UM - + UM).В отдельных случаях могут быть получены аналитические выражения, позволяющие упростить расчеты по переводу объекта в установившееся состояние, однако упрощений чаще всего не удается добиться в системах с нелинейностями и многими ограничениями. Поиск оптимального управления на любом шаге интегрирования выполняется по одному алгоритму, даже при наличии нелинейных ограничений.Синтез оптимального управления глубиной погруженияподводного аппаратаНеобходимо синтезировать закон изменения управляющего воздействия, обеспечивающего оптимальный по быстродействию переход с ограниченной скоростью без перерегулирования по положению подводного аппарата на заданную глубину.Модель объекта представляется обыкновенным дифференциальным уравнением [3]dV dtгде F - сила, действующая на подводный аппарат; m - масса аппарата; V - скорость движения аппарата; ^ - коэффициент сопротивления.Полагаем, что силовые воздействия на аппарат F1 и F2 имеют форму прямоугольных импульсов, воздействия внешней среды от изменения плотности и солености воды и сжатия корпуса незначительны, а коэффициент сопротивления не меняется во время переходного процесса. Импульс, создаваемый силой F1, разгоняет аппарат на погружение или всплытие, а импульс, создаваемый силой F2, тормозит движение аппарата, стабилизируя его на определенной глубине. Знаки Fi и F2 определяют направление перемещения аппарата.Рассмотрим алгоритм синтеза оптимального по быстродействию управления U подводным аппаратом при перемещении по глубине. Для поиска управления воспользуемся расчетом пробных шагов, выполняемых с предельными динамическими возможностями в сторону увеличения скорости аппарата и последующими переводами его в равновесные состояния. Анализ текущих значений скоростей и положений аппарата на модели позволит использовать рассчитанные для пробных шагов управления на реальном объекте.Аппарат перемещается из исходного положения L в заданное положение Lz при ограничении скорости движения на уровне Vm. Начальное значение пути торможения LT аппарата до останова после его разгона принимается равным нулю. Расчет управления U на очередном шаге At начинается с определения знака силового воздействия Fi, которое задает направление перемещения аппарата. При погружении аппарата начальное значение U принимает отрицательное значение Fi. Выполняется расчет пробного шага с целью получения максимального приращения скорости аппарата на очередном шаге. По имитационной модели подводного аппарата (14) определяется прогнозируемое значение скорости Vn после выполнения пробного шага с предельным управляющим воздействием:U-iy\v\Vn = V + ицчп At .(15) mЭта скорость сравнивается с допустимой скоростью движения аппарата Vm. В случае нарушения ограничения выполняется повторный расчет пробного шага по выражению (15) с управляющим воздействием U = 0. Таким способом реализуется ограничение по скорости.Определяется прогнозируемое перемещение аппарата Ln с начала движения при учете выполнения пробного шага с найденной скоростью Vn:Ln = L + \Vn\At .(16)Рассчитывается прогнозируемый путь торможения аппарата LTn со скорости Vn до останова после выполнения пробного шагаU-iy\v\Ьтп = LT + VnAt .(17)Т П F24Vn |Vn|Если Lz превышает сумму прогнозируемого перемещения аппарата Ln с начала движения при учете выполнения пробного шага с найденной скоростью Vn и прогнозируемого пути торможения аппарата LTn со скорости Vn до останова после выполнения пробного шага, то ранее полученное значение U используется для управления подводным аппаратом. После выполнения этого шага переменным аппарата V, L и LT присваиваются соответственно ранее найденные значения Vn, Ln и LTn. По описанной методике, начиная с выражения (15), рассчитываются управления для следующих шагов. Действия по такому циклу (15) - (17) продолжаются до тех пор, пока не будет выполнено условие начала торможения аппарата:Lz < (Ln + LTJI) .(18)В этом случае на подводный аппарат подается управление U, равное F2, и определяются значения скорости V и перемещения L после выполнения очередного шага:V + UAt:(19)mL = L + IV At .(20)Оценивается скорость перемещения аппарата V. Если продолжается движение, то на него подается управление U, равное F2, значения скорости и перемещения после выполнения очередного шага определяются по выражениям (19) и (20). После прекращения движения расчет скорости и перемещения аппарата ведется по выражениям (19) и (20) с управляющим воздействием U = 0. Таким способом исключается перерегулирование по положению.U 1000-100V 0,40-0,4-0,8-1,2 L40 0 -40 -80U 1000-100V 0,40-0,4-0,8-1,2 L 400-40-80абРис. 1. Переходные процессы при перемещении подводного аппарата с глубины 20 м на глубину: а - 80 м; б - 70 мПосле выхода на заданную глубину силовое воздействие к аппарату не прикладывается. Алгоритм позволяет определять для любого перемещения длительность приложения силового воздействия Fi, момент приложения силового воздействия F2 и начала торможения, момент снятия управляющего воздействия с подводного аппарата. После окончания действия первого импульса Fi скорость погружения аппарата снижается за счет сил сопротивления движению, а после окончания второго импульса F2 движение аппарата прекращается и он стабилизируется на глубине. Машинными экспериментами доказана работоспособность предложенного алгоритма управления подводным аппаратом с параметрами m = 1800, F1 = -120, F2 = 100 и ^ = 20 при различных глубинах погружения и ограничениях скорости.На рис. 1 приведены переходные процессы при перемещении подводного аппарата с глубины 20 м на глубину 80 м и 70 м при ограничении максимальной скорости на уровне 1,2 м/с. Алгоритм позволяет определять время погружения подводного аппарата при любых коэффициентах сопротивления, изменять во время переходных процессов ограничение на скорость перемещения и задания на глубину погружения. Он обеспечивает синтез оптимальных управляющих воздействий за малый интервал времени микропроцессорными средствами с учетом ограничения скорости и требования выхода подводного аппарата на заданную глубину без перерегулирования по положению.Оптимальное управление электродвигателем постоянного токаВ большинстве электромеханических устройств легко прослеживается взаимосвязанность составных частей, структурированность и подчиненность отдельных координат. Процессы протекают не в произвольном порядке, а в определенной последовательности, наблюдается строгая взаимозависимость координат. Синтезируем оптимальное по быстродействию управление двигателем постоянного тока независимого возбуждения по методике последовательного многошагового синтеза. Электродвигатель описывается уравнениямиU = iR + C со + L- ;dt(21)(22)где R, L - активное сопротивление и индуктивность якорной цепи двигателя; J, C - момент инерции и конструктивный коэффициент двигателя; со - частота вращения; U, i - напряжение и ток якоря; ic - ток статической нагрузки; t - время.Параметры двигателя считаются постоянными, а момент статической нагрузки - реактивным, т.е. вращение двигателя из заторможенного состояния может начаться только после превышения током якорной цепи статического тока. В качестве управляющего воздействия используется напряжение на якоре двигателя. Численным методом определяется управление U(t), обеспечивающее минимум времени перевода двигателя из исходного состояния с координатами U(0), i(0) и со(0) в новое состояние с координатами U(T), i(T) = ic и со(Г) = Qm. При этом значения частоты вращения со(г), тока i(T), напряжения U(T) и скорости измененияСоставлен алгоритм синтеза оптимального по быстродействию управления двигателем с учетом принятых ограничений. На начальном этапе определяется управление, необходимое для перевода объекта за один шаг из исходного состояния в заданное. При необходимости управление ограничивается на максимальном уровне. Определяются координаты объекта после выполнения пробного шага с найденным управлением. Затем выполняется расчет процесса изменения управления с целью выхода на установившееся значение тока. Производится проверка выполнения ограничения по току. В случае нарушения ограничения по току управление на пробном шаге корректируется. Рассчитывается процесс изменения с предельными возможностями тока двигателя до статического тока без перерегулирования. Таким способом объект переводится в равновесное состояние после выполнения пробного шага. Сравнивается полученное значение частоты вращения двигателя в равновесном состоянии с заданной частотой com . Нарушение ограничения по частоте вращения указывает на недопустимость использования управления пробного шага на реальном объекте. Его следует скорректировать и повторить расчеты по описанному циклу. Если отсутствует нарушение ограничения по частоте вращения, то найденное на пробном шаге управление считается оптимальным и используется на реальном объекте. На следующих шагах синтез управления выполняется по аналогичной методике. Алгоритм позволяет отрабатывать с предельным быстродействием возмущения, вводить ограничения на ускорение, учитывать изменения момента инерции и нагрузки.На рис. 3 представлен оптимальный по быстродействию процесс разгона электродвигателя с учетом принятых ограничений до частоты вращения 100 рад/с. Параметры системы: R = 1,5 Ом; L = 0,09 Гн; J = 0,0522 кг-м2; С = 1,28 В-с/рад;Изменение в ходе переходного процесса ограничений, параметров объекта или возмущающего воздействия просто учитываются при синтезе оптимального управления предложенным методом, так как поиск управления на каждом шаге ведется с использованием имитационного моделирования.Разработан алгоритм синтеза оптимального по быстродействию управления электродвигателем постоянного тока независимого возбуждения при разгоне без перерегулирования до заданной частоты вращения и питании его от реверсивного вентильного преобразователя с раздельным управлением и ограничением тока двигателя. Исследованы его возможности по воспроизведению заданных траекторий движения с ограничениями рывка и ускорения.ЗаключениеПредставленная методика формирования оптимальных законов управления применима в сложных нелинейных динамических системах с большим количеством ограничений. Использование дополнительных критериев позволяет выбирать управление на очередном шаге вычисления из ограниченного диапазона, увеличение числа нелинейностей и ограничений координат системы существенно не влияет на сложность алгоритма управления. Повышение порядка системы приводит к увеличению объема вычислений по циклическим алгоритмам, но не вызывает принципиальных трудностей. Для линейных систем с ограничением координат составлены простые алгоритмы, позволяющие переводить системы с предельным быстродействием в равновесное состояние и создавать на их основе системы супервизорного управления электромеханическими устройствами и технологическими процессами. Методика позволяет реализовывать предельные динамические возможности систем с учетом технологических требований и ограничений координат, автоматизировать сложные производственные процессы. Разработаны простые алгоритмы синтеза в реальном масштабе времени микропроцессорными средствами оптимальных по быстродействию управлений позиционными и следящими электроприводами с ограничением рывка, ускорения и скорости. Методика последовательного многошагового синтеза оптимальных управлений позволяет разрабатывать алгоритмы синтеза управлений электромеханическими устройствами, энергетическими системами и высокоскоростными подвижными объектами.ЛИТЕРАТУРА

Ключевые слова

limited coordinates , systems , synthesis , control , Method , системы , ограничения координат , управления , синтез , методика

Авторы

ФИООрганизацияДополнительноE-mail
Яковенко Павел Георгиевич ТГУ доктор технических наук pgj75@yandex.ru
Всего: 1

Ссылки

Дорри М.Х. Импульсное управление глубиной погружения подводного аппарата // Сб. трудов XXXIV Всерос. конф. «Управление движением морскими судами и специальными аппаратами». М.: ИПУ РАН, 2007. С. 20 - 25.
Мясников В.А., Игнатьев М.Б., Покровский А.М. Программное управление оборудованием. Л.: Машиностроение, 1974. 243 с.
Беллман Р. Динамическое программирование. М.: ИЛ, 1960. 400 с.
 УПРАВЛЕНИЕ ПЕРЕХОДНЫМИ ПРОЦЕССАМИ В СИСТЕМАХ С ОГРАНИЧЕНИЕМ КООРДИНАТ             | Вестник Томского государственного университета. Математика и механика. 2008. № 3 (4).

УПРАВЛЕНИЕ ПЕРЕХОДНЫМИ ПРОЦЕССАМИ В СИСТЕМАХ С ОГРАНИЧЕНИЕМ КООРДИНАТ | Вестник Томского государственного университета. Математика и механика. 2008. № 3 (4).

Полнотекстовая версия