Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
 Скачать 25.93 Mb.
|
|
(8.83) где (8.84) (8.85) Это есть уравнение окружности с радиусом R с центром, смещенным влево от начала координат на величину С. Задаваясь различными значениями M от 1 до ∞, можно построить семейство таких окружностей (рис. 8.26). На каждой окружности написано значение ординаты амплитудной частотной характеристики. При М = 1 окружность вырождается впрямую линию, параллельную оси ординат и проходящую слева от нее на расстоянии 0,5. При М окружность вырождается в точку, совпадающую сточкой. Для значений ординат амплитудной характеристики, лежащих в пределах 0 < М < 1, получается семейство окружностей, расположенных справа от линии М = 1, симметрично с первым семейством. При М = 0 окружность вырождается в точку, совпадающую с началом координат. Для построения амплитудной характеристики (рис. 8.25) достаточно в тех же координатах, где построены окружности М = const, нанести амплитудно-фазовую характеристику разомкнутой системы. Точки пересечения этой характеристики с окружностями будут определять точки амплитудной характеристики с соответствующими значениями ординат, равными М. Для определения показателя колебательности можно не строить амплитудную характеристику, так как достаточно знать одно максимальное значение ординаты M max , определяемое по наименьшей окружности М = const, которой коснется амплитудно-фазовая характеристика. Если при проектировании системы ставится условие, чтобы ее показатель колебательности был не больше некоторого заданного значения, например М = 1,5, то для выполнения этого необходимо, чтобы амплитудно-фазовая характеристика не заходила внутрь окружности, соответствующей этому значению М (рис. 8.27). Амплитудно-фазовая характеристика может только коснуться этой окружности. В этом случае показатель колебательности будет как раз равен заданному значению М max Таким образом, окружность Мшах является запретной зоной для амплитудно-фазовой характеристики разомкнутой системы. Эта зона охватывает точку (-1, j0) и обеспечивает получение заданного запаса устойчивости. Величина показателя колебательности может быть определена ив случае использования логарифмических частотных характеристик. Для этого отобразим запретную зону (рис. 8.27) на логарифмическую сетку. Рассмотрим отдельно окружность заданного показателя колебательности (рис. 8.28). На окружности возьмем произвольную точку В и построим вектор, соединяющий эту точку с началом координат. Установим для этого вектора связь между его модулем Аи запасом по фазе µ. Из треугольника ОВО 1 по теореме косинусов находим Далее можно найти и окончательно (8.86) Из рис. 8.28 нетрудно видеть, что зависимость (8.86) существует только для модулей, лежащих в пределах (8.87) В случае, когда запас по фазе может быть любым, так как в этом случае конец вектора не может попасть в запретную зону (рис. 8.28). Задаваясь различными значениями показателя М = const, а следовательно и С = const (8.84), по выражению (8.86) можно построить графики µ = f(A), которые носят название µ -кривых. Эти графики строятся обычно таким образом, что модуль А откладывается в децибелах (рис. 8.29). , Из выражения (8.86) можно найти, в частности, максимальный запас по фазе обычным методом отыскания максимума (8.88) Этот максимум получается, когда модуль . Если имеется построенная л.а.х. (рис. 8.30), то по имеющимся (µ -кривыми при заданном значении М можно построить требуемое значение запаса по фазе для каждого значения модуля. Это -построение должно делаться для модулей, лежащих в пределах (8.87). В результате будет получена запретная область для фазовой характеристики. Чтобы показатель колебательности был не больше заданного значения, фазовая характеристика не должна заходить в эту область. Нетрудно видеть, что определение качественного показателя, характеризующего запас устойчивости, делается здесь одновременно с определением устойчивости. Удобство показателя колебательности определяется также тем, что запас" устойчивости характеризуется здесь одним числом, имеющим для сравнительно широкого класса систем регулирования сравнительно узкие пределы (1,1-1,5). Оценка быстродействия может производиться по частотным характеристикам замкнутой и разомкнутой системы. При рассмотрении замкнутой системы обычно используется амплитудная частотная характеристика (рис. 8.25) или вещественная характеристика (рис. 8.6). Использование вещественной характеристики было рассмотрено выше (см. § 8.5). Для оценки быстродействия по амплитудной частотной характеристике (рис. 8.25) могут использоваться следующие величины р — резонансная частота, соответствующая пику а.ч.х.; п — частота, соответствующая полосе пропускания замкнутой системы и определяемая из условия А(ω п ) = 0,707; , с — частота среза, соответствующая условию Ас) = 1; э — эквивалентная полоса пропускания замкнутой системы, определяемая по выражению (8.89). где Эквивалентная полоса пропускания представляет собой основание прямоугольника (рис. 8.25), высота которого равна единице, а площадь равна площади под кривой квадратов модуля Ф. Понятие эквивалентной полосы пропускания тесно связано с вопросом пропускания системой помех, что будет рассмотрено в главе 11. В отличие от показателя колебательности, который является некоторой безразмерной характеристикой и лежит в сравнительно узких пределах, приведенные выше характерные частоты, определяющие быстродействие-системы, имеют размерность и их допустимые значения могут сильно меняться в зависимости от типа и назначения системы регулирования. Здесь наблюдается полная аналогия с критериями качества, основанными на рассмотрении кривых переходного процесса. Допустимое значение перерегулирования σ% (рис. 8.3) лежит в сравнительно узких пределах для систем самого различного назначения, а допустимое время переходного процесса t п может меняться от долей секунды до нескольких часов и более. Допустимые для данной системы регулирования значения р, пс или э должны устанавливаться для каждой конкретной системы на основе изучения условий ее эксплуатации. При этом характеризовать быстродействие системы может как вся совокупность указанных выше величин, таки каждая из них в отдельности. При определении быстродействия по частотной передаточной функции W(jω) разомкнутой системы может использоваться частота среза ω ср , которая определяется из условия равенства модуля единице или . Эта частота показана, например, на рис. 8.2 и 8.30. Определение частоты среза разомкнутой системы может быть сделано-на диаграмме, изображенной на рис. 8.26, по точке пересечения а.ф.х. с окружностью единичного радиуса, центр которой расположен вначале координат. Резонансная частота замкнутой системы р близка к частоте колебаний, системы в переходном процессе. Значение р может быть приближенно определено по точке а. ф. х. (рис. 8.26), которая ближе всего расположена' к точке (—1, j0). Частота среза ω ср во многих случаях близка к резонансной частоте системы ω р Удобной и наглядной мерой быстродействия системы является также частота сак (рис. 8.2), при которой задающее воздействие вида g = g max ·sinω k t отрабатывается системой с амплитудой ошибки не более х mах Хотя приведенные выше частотные критерии запаса устойчивости и быстродействия могут рассматриваться независимо от свойств системы регулирования во временной области, представляется полезным провести некоторое приближенное сопоставление частотных и временных характеристик. Если показатель колебательности М > 1, то замкнутую систему регулирования можно аппроксимировать колебательным звеном (см. § 4.5). Тогда передаточная функция замкнутой системы может быть представлена в виде (8.90) Для этой передаточной функции сравнительно просто найти, как зависят величины, которые определяют запас устойчивости перерегулирование σ%, показатель колебательности Ми запас устойчивости по фазе µ 1 , от параметра затухания ξ. Соответствующие кривые приведены на риса. На рис. 8.31, б дается зависимость между перерегулированием σ % и показателем колебательности М для той же передаточной функции (8.90). Кривые, приведенные на рис. 8.31, в некоторой мере характеризуют связь между показателями качества ив более сложных случаях, чем выражение (8.90). Так как резонансная частота р приблизительно соответствует частоте колебаний замкнутой системы в переходном процессе, то время достижения первого максимумам на переходной характеристике (рис. 8.3) может быть определено по приближенной зависимости (8.91) Если переходный процесс в системе заканчивается за 1—2 колебания, то время переходного процесса можно определить по приближенной зависимости (8.92) Сравнение формул (8.71) и (8.89) показывает, что эквивалентная полоса пропускания э совпадает с точностью до постоянного множителя с интегральной квадратичной оценкой /', определяемой формулами) и (8.68). Совпадение будет полным, если рассматривать всю эквивалентную полосу пропускания от и измерять ее в герцах. Тогда получаем (8.93) § 8.10. Чувствительность систем регулирования Действительные значения параметров системы регулирования практически всегда отличаются от расчетных. Это может вызываться неточностью изготовления отдельных элементов, изменением параметров в процессе хранения и эксплуатации, изменением внешних условий и т. д. Изменение параметров может привести к изменению статических и динамических свойств системы регулирования. Это обстоятельство желательно учесть заранее в процессе проектирования и настройки системы. Степень влияния изменения отдельных параметров на различные характеристики системы оценивается посредством чувствительности. Чувствительностью называется некоторый показатель, характеризующий свойство системы изменять режим работы при отклонении того или иного ее параметра от номинального,'или исходного, значения. В качестве оценки чувствительности используются так называемые функции чувствительности, представляющие собой частные производные й координаты системы по вариации го параметра, (8.94) или частные производные от используемого критерия качества / по j-му параметру, (8.95) Нулевым индексом сверху отмечено то обстоятельство, что частные производные должны приниматься равными значениям, соответствующим номинальным (расчетным) параметрам. Функции чувствительности временных характеристик. Посредством этих функций чувствительности оценивается влияние малых отклонений параметров системы от расчетных значений на временные характеристики системы (переходную функцию, функцию веса и др. Исходной системой называют систему, у которой все параметры равны расчетным значениями не имеют вариаций. Этой системе соответствует так называемое основное движение. Варьированной системой называют такую систему, у которой произошли вариации параметров. Движение ее называют варьированным движением. Дополнительным движением называют разность между варьированными основным движением. Пусть исходная система описывается совокупностью нелинейных уравнений первого порядка (8.96) Рассмотрим мгновенные вариации параметров , так что параметры приняли значения . Если изменения параметров не вызывают изменения порядка дифференциального уравнения, то варьированное движение будет описываться совокупностью уравнений (8.97) Для дополнительного движения можно записать (8.98) При условии дифференцируемости по параметрам дополнительное движение можно разложить вряд Тейлора. Для малых вариаций параметров допустимо ограничиться линейными членами разложения. Тогда получим уравнения первого приближения для дополнительного движения (8.99) Частные производные, находящиеся в скобках, должны быть равны их значениям при Таким образом, первое приближение для дополнительного движения может быть найдено при известных функциях чувствительности. Заметим, что использование функций чувствительности удобнее для нахождения дополнительного движения по сравнению с прямой формулой (8.98), так как последняя во многих случаях может дать большие ошибки вследствие необходимости вычитать две близкие величины. При значительных вариациях может оказаться необходимым использование второго приближения с удерживанием в ряде Тейлора, кроме линейных, также и квадратичных членов. Дифференцирование исходных уравнений (8.96) по a j приводит к так называемым уравнениям чувствительности Решение этих уравнений дает функции чувствительности u ij . Однако уравнения (8.100) оказываются сложными и решение их затруднительно, Более целесообразен путь структурного построения модели, используемой для нахождения функций чувствительности [21, 53, 111] Обратимся теперь к линейным системам. Не снижая общности рассуждений, можно рассматривать случай изменения одного го параметра. В некоторых случаях функции чувствительности получаются дифференцированием известной функции времени на выходе системы. Так, если передаточная функция системы соответствует апериодическому звену второго порядка, то (см. табл. 4.2) При поступлении на вход ступенчатой функции на выходе будет Пусть, например, вариацию претерпевает постоянная времени Т. Тогда дифференцирование последнего выражения по Т даст функцию чувствительности поэтому параметру Дополнительное движение при этом будет - вариация постоянной времени Т 3 Пусть рассматриваемая система описывается совокупностью уравнений первого порядка (8.101) где - постоянные коэффициенты, x i — фазовые координаты, а f q (t) — внешние воздействия. Начальные условия в системе при t=0 . Уравнения чувствительности получаются из (8.1.01) дифференцированием по варьируемому параметру а j , от которого могут зависеть коэффициенты аи) где - частные производные от коэффициентов системы уравнений (8.101) по варьируемому параметру а j . Уравнениям (8.102) соответствуют начальные условия . Если начальные условия не зависят от параметра a j , то уравнениям (8.102) соответствуют нулевые начальные условия. Для решения (8.102) необходимо предварительно решить совокупность уравнений (8.101) и определить исходное движение Для нахождения функций чувствительности и дополнительного движения удобно использовать передаточные функции системы. Пусть, например, регулируемая величина связана с задающим воздействием зависимостью (8.103) где р) — изображение задающего воздействия. Функция чувствительности может быть получена из (8.103) дифференцированием по параметру a j : (8.104) Здесь введена функция чувствительности передаточной функции (8.105) которая определяет первое приближение дополнительной-передаточной функции, равной разности варьируемой и исходной передаточных функций при вариации параметра a j (8.106) Эти зависимости справедливы в том случае, когда вариация параметра a j - не меняет порядка характеристического уравнения системы. Может также использоваться так называемая логарифмическая функция чувствительности Формула (8.107), строго говоря, может использоваться в тех случаях, когда Фра) и a представляют собой безразмерные величины. Если эти величины размерны, то их логарифмирование возможно, если использовать прием, указанный в § 4.4. Найдем дополнительную передаточную функцию для случая, когда исходная передаточная функция может быть представлена в виде отношения двух полиномов (8.108) где - вариации полиномов числителя и знаменателя передаточной функции. Формула (8.108) позволяет составить структурную схему модели чувствительности в виде, изображенном на рис. 8.32. Эта схема может быть использована для нахождения функции дополнительного движения или функции чувствительности и ; расчетным путем или моделированием на ЭВМ. Составим, например, модель чувствительности для передаточной функции замкнутой системы (8.109) при вариации параметра τ. В соответствии с изложенным находим Равенство приращений числителя и знаменателя Фр) позволяет упростить схему модели. Она изображена на риса в исходном, а на рис. 8.33, б — в преобразованном виде. В общем случае, когда передаточная функция зависит отряда варьируемых параметров, дополнительная передаточная функция (8.110) Если к системе приложено несколько внешних воздействий , то следует найти дополнительные передаточные функции для всех исходных передаточных функций, определенных для каждого внешнего воздействия. Функции чувствительности критериев качества. Если в системе произошли изменения ряда параметров , то результирующее изменение некоторой используемой оценки качества (8.111) где — варьированное значение оценки качества, а I — ее исходное значение, можно подсчитать по формуле полного дифференциала (8-112) Так как в большинстве случаев известны только вероятностные оценки вариаций , то целесообразно использование вероятностных методов. Так, если известны максимальные возможные отклонения , то при их независимости друг от друга можно найти среднеквадратичный максимум отклонения оценки качества (8.113) и среднеквадратичный относительный максимум (8.114) Если заданы дисперсии отклонений параметров и отклонения независимы, то можно найти дисперсию оценки качества (8.115). В качестве критериев оценки качества системы могут использоваться, например, максимум ошибки, коэффициенты ошибок, оценки запаса устойчивости и быстродействия, интегральные оценки и т. п. Пример. Пусть передаточная функция разомкнутой системы имеет вид Требуется определить среднеквадратичный максимум отклонения показателя колебательности, если сек, причем изменения параметров независимы. Определим вначале исходное значение показателя колебательности. Для этого необходимо найти максимум модуля частотной передаточной функции замкнутой системы Исследование на максимум дает при показатель колебательности M= 1, при КТ > 2 показатель колебательности Функции чувствительности, если , Среднеквадратичный максимум отклонения (8.113) Таким образом, в рассматриваемой системе показатель .колебательности ГЛАВА 9 ПОВЫШЕНИЕ ТОЧНОСТИ СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ § 9.1. Общие методы К числу общих методов повышения точности систем автоматического регулирования относятся 1) увеличение коэффициента усиления разомкнутой цепи 2) повышение степени астатизма 3) применение регулирования по производным от ошибки. Увеличение общего коэффициента усиления разомкнутой цепи является наиболее универсальными эффективным методом. Увеличить общий коэффициент усиления можно обычно за счет введения в систему регулирования усилителей. Однако в некоторых случаях удается достичь этого увеличения за счет повышения коэффициентов передачи отдельных звеньев, например чувствительных элементов, редукторов и т. д. Увеличение общего коэффициента усиления благоприятно сказывается в смысле уменьшения ошибок практически во всех типовых режимах. Это вытекает, в частности, из того, что общий коэффициент усиления разомкнутой цепи входит в качестве делителя вовсе коэффициенты ошибок (см. пример, рассмотренный в § 8.3). Однако увеличение общего коэффициента усиления ограничивается устойчивостью системы регулирования. При повышении коэффициента усиления, как правило, система приближается к колебательной границе устойчивости. При некотором предельном его значении в системе возникают незатухающие колебания. В этом сказывается противоречие между требованиями к точности и требованиями к устойчивости системы регулирования. В связи с этим повышение общего коэффициента усиления до значения, при котором обеспечивается выполнение требований к точности, обычно может производиться только при одновременном повышении запаса устойчивости системы, что осуществляется при помощи так называемых корректирующих средств, рассматриваемых в следующей главе. Повышение порядка астатизма. Повышение порядка астатизма используется для устранения установившихся ошибок в различных типовых режимах в неподвижном положении, при движении с постоянной скоростью, при движении с постоянным ускорением и т. д. Формально это сводится к тому, чтобы сделать равными нулю первые коэффициенты ошибки системы, например, с = 0 при астатизме первого порядка, или с = с = 0 при астатизме второго порядка, или с = с = с = 0 при астатизме третьего порядка и т. д. Физически повышение порядка астатизма осуществляется за счет введения в канал регулирования интегрирующих звеньев. В качестве таких звеньев могут, например, использоваться звенья, изображенные на рис. 4.21. Структурная схема системы регулирования с введенным интегрирующим звеном изображена на рис. 9.1. Передаточная функция интегрирующего звена где |