Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
Скачать 25.93 Mb.
|
wt = ψ (19.6) Отсюда видно, что в общем случае все три коэффициента являются функциями трех неизвестных (19.7) В частных случаях эти зависимости могут быть более простыми. Подстановка выражений (19.4) ив заданное дифференциальное уравнение (19.1) с учетом свойства фильтра (см. § 18.2) дает Это уравнение разбивается на два (19.8) (19.9) При таком разделении сохраняются существенно нелинейные свойства и отсутствие суперпозиции решений, так как остается нелинейная взаимосвязь обоих уравнений через соотношения (19.7). Можно предложить два метода решения задачи. Первый метод состоит в следующем. Уравнение (19.9) совпадает с прежним уравнением (18.33); отличие состоит лишь в том, что теперь коэффициенты q и q' согласно (19.7) зависят не только от аи, но и от смещениях. Поэтому, написав как прежде, характеристическое уравнение (19.10) заменив р на jw и выделив вещественную и мнимую части, в отличие от (18.36), получим здесь два алгебраических уравнения стремя неизвестными (19.11) Эти уравнения дают возможность определить амплитуду а пи частоту w п автоколебаний как функции постоянной составляющей х (19.12) Для решения этой задачи можно применять любой из способов, описанных в § 18.2, в зависимости оттого, какой из них лучше подходит к условиям заданной конкретной задачи. Таким же способом можно определить зависимость аи не только от х, но и от параметров системы с целью выбора последних. Что касается тех способов § 18.2, где используются графики аи а, то здесь их необходимо строить в виде серии кривых при разных постоянных значениях х (рис. 19.3). После того как из уравнений (19.11) определены зависимости (19.12), можно, воспользовавшись первым из выражений (19.7), найти функцию смещения (19.13) Подставив ее в (19.8), получим алгебраическое уравнение (19.14) с одной неизвестной х, которая отсюда и определяется. Чаще всего это уравнение относительно х является трансцендентными решается графически. Затем согласно (19.12) определяются также амплитуда а пи частота w п См. конкретные формулы для разных нелинейностей в § 19.3. Указанную зависимость (19.12) амплитуды и частоты автоколебаний от величины смещения центра колебаний, характерную именно для нелинейных систем (в линейных системах смещение центра колебаний не играет роли, надо всегда иметь ввиду. При одних нелинейностях она может быть весьма существенной, при других — менее существенной. Второй метод решения той же задачи состоит, наоборот, в том, что сначала решается уравнение (19.8), где согласно (19.7) будет F° ха) или часто F° (ха. Решение получает вид (19.15) Это решение подставляется затем в уравнения (19.11), которые, таким образом, будут содержать только две неизвестные а пи п. Определив последние (по любому из способов § 18.2), вычисляем потом пои величину х, которая будет в результате зависеть от формы нелинейности, от параметров системы и от внешнего воздействия М. Величинах и является искомой статической или скоростной ошибкой соответственно для статической и астатической систем. В тех случаях, когда передаточная функция линейной части системы р) имеет нулевой корень в знаменателе (что часто бывает, те. когда Q(0) = 0, вместо (19.14) получаем уравнение (19.16) откуда определяется статическое отклонение или скоростная ошибках (МВ случае, когда при отсутствии внешнего воздействия (М = 0) определяются автоколебания в системе с несимметричной нелинейностью, те. нелинейностью х) или же х, рх), для которой (19.17) вместо уравнения (19.8) получаем (19.18) Оно решается любым из тех же двух методов, описанных выше для уравнения (19.8). Одновременно согласно (19.11) определяются ха п, w п Если в этом случае знаменатель р) передаточной функции линейной части системы имеет нулевой корень, то Q(0) = 0 и, следовательно, уравнение (19.18) с учетом (19.13) принимает вид (19.19) откуда определяется х. Это означает, что в указанных системах возникает такое смещение х колебаний переменной х, которое ликвидирует свойственную данной нелинейности несимметрию колебаний переменной те. обеспечивается F° = 0), как показано, например, на рис. 19.4 в отличие от рис. 19.1, б. Приведем пример исследования совместного влияния двух внешних воздействий, причем из дальнейшего будет видно, что, в отличие от линейных систем, здесь нельзя просто складывать статические ошибки от отдельно взятых воздействий. Уравнения автоматической системы заданы в виде (19.20) (19.21) (19.22) где х) — простейшая симметричная релейная характеристика, показанная на рис. 19.5: (19.23) Чтобы воспользоваться выведенными выше общими формулами, надо сначала привести заданную систему уравнений (19.20) — (19.22) к одному уравнению типа (19.1). В результате получаем (19.24) Пусть f1(t) является задающим воздействием, изменяющимся с постоянной скоростью (19.25) которое требуется воспроизвести на выходе системы в виде x 4 (t). Допустим также, что второе внешнее воздействие f2(t) является возмущающими имеет постоянную величину например, постоянная нагрузка на выходном валу системы (19.26) Его влияние требуется свести к минимуму. Найдем установившуюся ошибку на выходе системы. Правая часть уравнения (19.24) будет при этом постоянной, и установившееся решение для х с учетом автоколебаний в виде (19.27) Гармоническая линеаризация нелинейности (19.23) при этом согласно (19.6) дает (19.28) В данной задаче согласно (19.24) уравнение (19.9) для периодических составляющих получит вида уравнение (19.8) для постоянных составляющих будет (19.30) где согласно (19 24)— (19.26) (19.31) Выше были указаны два метода решения задачи. Для иллюстрации обоих методов решим данную задачу каждым из них. Согласно первому методу сначала решается уравнение (19.29) для определения зависимостей ахи х. Характеристическое уравнение здесь будет (19.32) и уравнения (19.11) поэтому примут вид (19.33) Исключая отсюда q, находим частоту автоколебаний (19.34) Частота w в данной задаче оказалась независящей от смещениях, а следовательно, и от величины внешнего воздействия. Затем, подставляя в первое из уравнений (19.33) выражение q из (19.28) и w 2 из (19.34), получаем биквадратное уравнение для отыскания зависимости амплитуды автоколебаний аи от смещениях) где величина (19.36) представляет собой амплитуду автоколебаний в данной системе при отсутствии смещения при х = 0). Отсюда (19.37) Полученное выражение можно записать также в виде (19.38) если обозначить (19.39) Результат (19.37) или (19.38) и представляет собой искомую зависимость а п (x°). Далее, согласно первому методу решения задачи подставим полученное значение амплитуды а п изв выражение (19.28) для F°, откуда с использованием (19.39) найдем функцию смещения (19.40) где А определяется через параметры системы формулой (19.36). Подставив величину (19.40) в уравнение (19.30) для постоянных составляющих, с учетом (19.31) и (19.36) получим (19.41) Сравнивая это с формулой (19.39), видим, что для искусственно введенной ранее величины а можно записать следующее выражение (19.42) Эта величина характеризует совокупность приложенных к системе внешних воздействий. Учитывая это, из формулы (19.38) находим амплитуду автоколебаний (19.43) Существенно то, что амплитуда автоколебаний зависит не только от параметров системы см. (19.36)), но еще и от величины внешних воздействий. Эта зависимость нелинейная. В данном случае при увеличении внешних воздействий амплитуда уменьшается по закону косинуса, в то время как частота не зависит от внешних воздействий. Из формулы (19.43) видно, что автоколебания существуют до тех пор, пока величины внешних воздействий удовлетворяют условию (19.44) При этом амплитуда автоколебаний изменяется в пределах А > a п > 0. Таков первый метод определения установившихся величин смещениях, амплитуды а пи частоты w п автоколебаний при наличии внешних воздействий. Проиллюстрируем также и второй метод. Согласно второму методу сначала решается уравнение (19.30). Пои первой из формул (19.28) находим или, с учетом (19.31), (19.45) Для отыскания входящей сюда амплитуды а воспользуемся уравнением (19.29). Характеристическое уравнение для него будет (19.32), и уравнения (19.11) поэтому примут вид (19.46) где согласно (19.28) и (19.45) (19.47) Исключая из уравнений (19.46) величину q, находим частоту автоколебаний (19.48) Подставив найденные выражения q и w п в первое из уравнений (19.46), найдем амплитуду автоколебаний (19.49) где величина (19.50) является амплитудой автоколебаний при отсутствии внешних. Подставив найденное выражение амплитуды (19.49) в формулу (19.45), получим окончательно величину смещения (19.51) Как видно, второй метод в данной задаче приводит к тем же самым результатам значительно более коротким путем, чем первый, что очень важно для практических расчетов (принципиально же оба метода эквивалентны друг другу. По-видимому, большая простота второго метода будет иметь место ив большинстве других задач. В этом втором методе, в отличие от первого, функция смещения Ф(х°) не определяется. Однако последняя может понадобиться в дальнейшем для других целей. Но ее тоже легко можно определить при использовании второго метода. Здесь величины ха пи их отношение выражены через величины внешних воздействий. Функция же смещения Ф хне должна содержать ни величин внешних воздействий, ни амплитуды а п, зависящей от них. Подставив значение квадратной скобки изв, получим а подставив это в функцию смещения (19.52) первую из формул (19.28), сразу получим искомую (19.53) где А выражается только через параметры системы согласно (19.50). Важно отметить, что функция смещения Ф(х°) не зависит ни от числа внешних воздействий, ни от характера их изменения (если они постоянные или медленно меняющиеся, что наиболее наглядно было видно из первого метода решения задачи. Итак, двумя разными методами определена величина смещениях автоколебаний на входе реле. Найдем теперь установившуюся ошибку на выходе системы x4. Поскольку на выходе должно воспроизводиться внешнее воздействие f 1 (t), то согласно рис. 19.5 и второму уравнению (19.20) ошибка данной системы выражается величиной х, установившееся решение для которой, следовательно, и надо искать. Выразив переменную х через х, которая уже известна, из заданных уравнений системы (19.20) и (19.21) получаем Учитывая (19.25) и (19.27), перепишем данное уравнение в виде (19.54) В соответствии с видом правой части установившееся решение этого линейного уравнения следует искать в виде (19.55) где хи с — постоянные, ах периодическая составляющая. Подставив это в (19.54), получим три уравнения для отыскания указанных величин (19.56) (19.57) (19.58) Второе из них дает (19.59) Тогда из (19.56) находим (19.60) где x° определяется формулой (19.51) через внешние воздействия. Наконец, из уравнения (19.58) получаем амплитуду автоколебаний переменной x4: (19.61) где а п определяется формулой (19.49) через внешнее воздействие, а w п — формулой (19.48). Итак, в данной системе имеются все три составляющие ошибки (19.55), зависящие от величины внешних воздействий и от параметров системы. Наиболее нежелательной из них является составляющая с, возрастающая пропорционально времени. Поэтому систему необходимо видоизменить в первую очередь так, чтобы уничтожить эту составляющую ошибки, те. сделать с = 0. Для этого можно было бы вовсе изъять дополнительную обратную связь (рис. 19.5), так как при ос = 0 согласно (19.59) будет с 0. Однако при этом существенно возрастает амплитуда автоколебаний (19.61), те. периодическая составляющая ошибки. Поэтому более целесообразной мерой будет замена жесткой обратной связи хос = k ос x 4 на гибкую хос = k ос р. Тогда в уравнении (19.54) величина k ос заменится на k ос р (19.62) Как видим, составляющая, пропорциональная времени, в правой части уравнения исчезла, вследствие чего установившееся решение для ошибки х, в отличие отбудет) При этом из (19.62) находим (19.64) (19.65) причем изменяются, конечно, и формулы для x°, а пи п (их можно получить таким же способом. Подбором параметров системы амплитуду автоколебаний ошибки a1 можно сделать весьма малой. § 19.2. Прохождение медленно меняющихся сигналов в автоколебательных системах Рассмотрим очень важный для практики случай, когда внешнее воздействие f(t), которое может быть либо возмущающим, либо управляющим (задающим, в автоколебательной системе является непостоянным, а медленно меняющимся. Медленно меняющейся будем называть такую функцию времени, которая сравнительно мало изменяется за период автоколебаний, те. соблюдается условие в виде одного из неравенств Где пап частота автоколебаний. Соответственно для астатических систем медленно меняющейся скоростью f'(t) будет такая, для которой выполняется условие Указанными свойствами почти всегда обладают полезные сигналы управления, проходящие через автоматическую автоколебательную систему (в том числе в переходных процессах. Условие медленного изменения любой функции времени можно выразить также ив частотной форме, а именно медленно меняющейся считается такая функция, возможные частоты изменения которой во времени значительно ниже возможной частоты возникающего в системе периодического решения (автоколебаний. Сделанные предположения позволят величину f(f) или, соответственно, р) считать постоянной за время каждого периода исследуемых автоколебаний и искать решение в той же форме (19.4); где t' отсчитывается отдельно внутри каждого периода, ибо теперь ха и w будут непостоянными, а переменными во времени t (от периода к периоду) вместе с изменением внешнего воздействия / (г. При этом х (t) будет медленно меняющимся сигналом на входе нелинейности. В связи со сказанным здесь остается в силе разложение (19.5) — (19.7). Но подстановка его в заданное уравнение нелинейной автоматической системы (19.1) дает При достаточно медленном изменении функции f(t) (а в астатических системах р) и величина, входящих в коэффициенты F°, q, q', данное уравнение может быть разделено на два отдельных уравнения (19.66) (19.67) соответственно для медленно меняющейся составляющей и для колебательной составляющей. При этом разделении уравнений, как и прежде, сохраняются существенно нелинейные свойства системы. Следовательно, здесь сохраняется целиком прежний (§ 19.1) первый метод решения задачи (второй здесь неприемлем, выраженный формулами (19.10) — (19.13), где в данном случае х является величиной непостоянной, а медленно меняющейся. Поэтому прежний процесс решения заканчивается определением функции смещения (19.13). Подставив (19.13) в (19.66), получим дифференциальное уравнение для определения медленно меняющегося сигнала управлениях) (на фоне автоколебаний системы) в виде (19.68) Таким образом, получается, что для определения медленно меняющихся процессов функцию смещения (19.69) следует подставить в уравнение автоматической системы (19.1) вместо заданной нелинейности Р (х, рх). Следовательно, функция смещения Ф (х) представляет собой как бы статическую характеристику (обычно криволинейную, которая определяет зависимость между выходной и входной величинами заданной нелинейности для постоянных или медленно меняющихся сигналов в автоколебательной системе. При любых нелинейностях, в том числе и скачкообразных, функция смещения Ф(х°) может получать при определенных условиях вид весьма плавной кривой. Этот эффект называется вибрационным сглаживанием нелинейностей при помощи автоколебаний, а функцию смещения Ф(х°) можно называть сглаженной нелинейной характеристикой. Так, в примере § 19.1 согласно (19.40) функция смещения будет иметь вид риса, те. для медленно меняющегося сигнала в данной релейной системе нелинейная характеристика будет в определенных пределах иметь плавный вид (риса) вместо скачкообразного (рис. 19.6, 6) — за счет сглаживающего влияния автоколебательных вибраций. Далее, например, для нелинейностей, обусловленных зоной нечувствительности (риса, а также зазором (рис. 19.7, в) и петлей, сигналы х < b при отсутствии автоколебаний не передаются (F = 0). При наличии же автоколебаний сигнал х < b передается в виде составляющей Р. Поэтому для медленно меняющегося сигнала получается плавная характеристика (функция смещения) Ф (х) без зоны нечувствительности (рис. 19.7, б. Эффект вибрационного сглаживания нелинейностей в этих примерах является положительным (ликвидация зон нечувствительности и петель. Однако в других случаях эффект вибрационного сглаживания нелинейности может оказаться и отрицательным. Возьмем, например, нелинейную характеристику с зоной насыщения (ограниченно-линейную), показанную на рис. 19.8; В этом случае за счет того, что верхушки синусоиды с одной стороны срезаются, постоянная составляющая F° будет меньше, чем само значение F, соответствующее линейному начальному участку. Поэтому постоянный или медленно меняющийся сигнал будет при наличии автоколебаний проходить через данную нелинейность с меньшим коэффициентом усиления, чем без автоколебаний, что может в известных случаях отрицательно сказаться на качестве автоматической системы в целом. Во многих случаях вычисление аи будет необходимо только сточки зрения проверки выполнения условий вибрационного сглаживания нелинейности и допустимости таких вибраций в данной конкретной автоматической системе. Основными же для качества работы автоматической системы при этом будут являться медленно меняющиеся процессы, определяемые уравнением (19.68). Сточки зрения упрощения их определения весьма важными являются следующие два обстоятельства. Во-первых, вид функции смещения Ф (х, как видно из § 19.1, не зависит ни от количества и места приложения внешних воздействий на систему, ни от характера их изменения (если только они медленно меняющиеся. Вид Ф (x°) зависит от формы нелинейности, от структуры и от параметров системы. Поэтому можно пользоваться любым методом определения Ф (х) при любых частных упрощающих предположениях относительно внешних воздействий. Можно, например, пользоваться более простым вторым методом из описанных в § 19.1 методов и проиллюстрированных там на примере, взяв любое одно постоянное по величине внешнее воздействие. Во-вторых, какова бы ни была заданная нелинейность х, рх) (скачкообразная, петлевая и т. п, обычно функция смещения Ф (х) получает вид плавной кривой. Поэтому, в отличие от первоначально заданной нелинейности, ее легко можно линеаризовать обычным способом (по касательной или но секущей вначале координат или в другом начале отсчета. Имея ввиду это свойство, часто вместо термина вибрационное сглаживание употребляют термин вибрационная линеаризация (будем придерживаться первого из них. Итак, в определенном диапазоне можно считать где (19.70) Графически k н представляет собой тангенс угла наклона прямой (касательной или секущей, рис. 19.7, б. Величина коэффициента k н зависит от соотношения, вообще говоря, всех параметров системы. Например, для системы, описываемой уравнениями (19.20) - (19.23), согласно (19.53) и (19.50) имеем (19.71) Поэтому для расчета медленно протекающих процессов в данной системе на основании (19.20) — (19.22) и (19.70) получаем линейные уравнения (19.72) или единое линейное уравнение (19.24), в котором надо заменить хна хи хна н x 0 Определение коэффициента усиления k н можно значительно упростить следующим образом. Поскольку (функция смещения Ф (х) определяется согласно (19.13) и (19.7) по выражению ха, в которое подставлена зависимость ах, то формулу для вычисления k н можно представить в виде В тех случаях, когда рассматриваются нечетно-симметричные нелинейности х, величина F° не зависит от w и, кроме того, согласно (19.6) так как производная под знаком интеграла будет четной функцией. Следовательно, для нечетно-симметричных нелинейностей х, как однозначных таки петлевых, величину k н можно вычислять по формуле непосредственно из выражения (19.7), не определяя функции смещения Ф (х, Во многих задачах это будет существенным упрощением решения. Это упрощение не относится к несимметричным нелинейностям, а также к тем случаям, когда k н приходится определять не по касательной, а по секущей. Итак, с подстановкой (19.70) уравнение (19.68) для определения медленна протекающих процессов становится обыкновенным линейным уравнением (19.73) я, как таковое, легко решается. Как видим, введенный здесь принцип разделения уравнений для колебательных и для медленно меняющихся составляющих, при котором сохраняются существенно нелинейные свойства системы, приводит к весьма важным для практических расчетов результатам. Существенным выводом является то, что медленно меняющиеся сигналы проходят через нелинейность с другим коэффициентом усиления k н, чем автоколебания. Особенно важно использовать свойство вибрационного сглаживания нелинейностей с последующей их обычной линеаризацией при расчете сложных автоматических систем. Если, например, система автоматического управления полетом самолета работает по схеме, изображенной на рис. 19.9, то часть системы, обведенную пунктиром (релейный усилитель, приводи дополнительная обратная связь, как отдельную следящую систему, можно рассчитывать изложенным выше методом с учетом автоколебательных вибраций. Частоту последних путем соответствующего выбора параметров этой части системы или введением корректирующих устройств можно сделать достаточно большой стем, чтобы амплитуда автоколебаний переменной хна выходе этой части системы была мала. Если же указанную амплитуду x2 не удается сделать малой (тогда руль будет колебаться, то необходимо, чтобы указанная частота практически не воспринималась корпусом самолета в процессе его движения вокруг центра тяжести. Тогда расчет автоматической системы будет выглядеть следующим образом. Автоколебания определяем только в обведенной пунктиром внутренней части (рис. 19.9), как в отдельной самостоятельной системе, считая х) произвольным медленно меняющимся внешним входным воздействием, ах выходной величиной. Для такой простой системы находим, как изложено выше, функцию смещения Ф (ха также частоту и амплитуду автоколебаний в зависимости от величины внешнего воздействия. Выбираем параметры данной части системы так, чтобы условия вибрационного сглаживания нелинейности соблюдались во всем практически возможном диапазоне изменения входной величины х. При этом следим затем, чтобы частота автоколебаний системы (зависящая от параметров системы) лежала за пределами частот возможных колебаний самолета (чтобы она практически не воспринималась корпусом самолета. После такого расчета внутренней части системы производим обычную линеаризацию функции смещения Ф (х, те. заменяем ее одной прямой линией F° = k н (причем можно использовать указывавшееся упрощение в определении k н. В результате получаем линейное уравнение для медленно протекающих процессов в данной части системы. К этому уравнению добавляем уравнение всей остальной части системы (в данном случае самолета, чувствительных элементов и руля, см. рис. 19.9) и рассчитываем всю систему в целом, как линейную, по любым обычным методам теории автоматического регулирования. При этом не обращаем уже внимания на автоколебания, которые локализуются в рассчитанном ранее внутреннем контуре системы. Однако их влияние не игнорируется, ибо оно было учтено при определении функции смещения Ф (хи коэффициента k н Изложенный принцип позволяет, во-первых, вести расчет автоколебаний поболее простым уравнениям (так как выделяется только внутренняя часть системы) и, во-вторых, значительно упрощает расчет всей системы в целом, сводя его к исследованию обыкновенных линейных уравнений (нос коэффициентом k н, зависящим от автоколебаний, те. от параметров внутреннего контура системы. Если необходимо учесть переменные коэффициенты и нелинейности самого самолета, то уравнения системы в целом не будут уже столь простыми. Однако ив этом случае полностью сохраняет смысл предварительный отдельный расчет внутреннего контура системы, так как влияние нелинейностей самого самолета будет распространяться обычно только нате процессы движения, которые по сравнению с автоколебаниями внутреннего контура являются медленными. Известно, что и самодвижение самолета, например по тангажу, можно разделить на два, одно более быстрое — движение около центра тяжести (угловое движение) и другое более медленное — движение центра тяжести (движение по траектории. Оба они являются медленными по сравнению с автоколебаниями внутреннего контура системы управления. Однако их в свою очередь тоже можно рассматривать отдельно. Следовательно, в этом случае, кроме обычно применяемого пространственного разбиения движения самолета по каналам (тангажа, курса, крена, расчет системы по каждому каналу (например, тангажа) разбивается еще натри этапа по степени медленности движения во времени. Аналогичное разделение расчета по крайней мерена два этапа по степени медленности во времени бывает целесообразными для многих других нелинейных автоматических систем регулирования, слежения, стабилизации и т. п. В сложных системах такие приемы, существенно упрощающие все исследование, оказываются единственными, которые могут сделать расчет системы практически осуществимым. Важно иметь ввиду, что при этом принципе разделения движений сохраняется существенная нелинейная взаимосвязь между ними. |