Исследование устойчивости - (реферат)
p>Следует отметить, что критерий Рауса неприменим в двух случаях. Во-первых, когда какой-либо элемент первого столбца, начиная со второго, равен нулю. Тогда все члены следующей строки будут равны бесконечности. Во-вторых, когда все элементы второй или любой из следующих строк равны нулю. В этих специальных случаях необходимо использовать для анализа другие методы.Для примера рассмотрим уравнение:
Сопоставляя (43) с (17), можно записать
Тогда таблица Рауса будет иметь следующий вид:
Замечаем, что знак элементов первого столбца таблицы сначала изменяется с плюса на минус, а затем–опять на минус. Это означает, что уравнение имеет два корня с положительными вещественными частями. Действительно, корнями уравнения (43) являются: Следует иметь в виду, что для упрощения вычислений можно разделить (или умножить) все элементы любой строки на положительное число, прежде чем использовать их для получения следующей строки. Очевидно, что такая операция не изменит знака членов следующей строки и не отразится на конечном результате. Например, элементы третьей строки таблицы (45) можно было бы разделить на 8 для упрощения последующих вычислений.
Анализ результатов устойчивости в нелинейных системах.
При исследовании устойчивости в цепях постоянного тока при малых возмущениях обнаружение неустойчивости возможно только при наличии элементов с отрицательным дифференциальным сопротивлением. Элементы с “падающими” участками на вольт - амперной характеристике, как известно, разделяются на две группы. К первой группе относятся элементы с S –образной ВАХ, например, электрическая дуга. У них ток относительно напряжения является неоднозначной функцией, т. е. при определенных напряжениях при их плавных изменениях теоретически возможны резкие изменения, так называемые “скачки” тока. Однако, опыт показывает, что при этом в такой цепи всегда присутствует небольшая паразитная индуктивность, которая “сглаживает” скачкообразные изменения, не допуская скачкообразного изменения энергии магнитного поля, поскольку оно в реальной системе невозможно. По этой причине при исследовании устойчивости в таких цепях нелинейное сопротивление сS –образной ВАХ представляют в виде последовательного соединения его дифференциального сопротивления и малой начальной индуктивности. У элементов с N –образной ВАХ напряжение является неоднозначной функцией тока. Здесь при определенных токах при их плавных изменениях теоретически возможны скачкообразные изменения напряжения, которые, однако, в реальных системах предовращаются наличием малой паразитной емкости. Поэтому при исследовании устойчивости в таких цепях нелинейное сопротивление с ВАХN –типа заменяют эквивалентной схемой с параллельным соединением дифференциального сопротивления и малой паразитной емкости.
Устойчивость точки равновесия электрической дуги.
В качестве примера цепи постоянного тока рассмотрим электрическую дугу с ВАХ S - типа. С учетом паразитной индуктивностиLцепь имеет эквивалентную схему, представленную на рис. ? ?? ?? ?. Уравнение цепи на основании закона Кирхгофа имеет вид:
В состоянии равновесия ток в цепи i не должен изменяться во времени t: Поэтому уравнение примет вид:
Запишем уравнение возмущенного движения, полагая при этом
где о(t) – малое возмущение; тогда получим:
Разложим функцию U(i) в ряд в точке I0 по степеням малых возмущений ои ограничимся величинами первого порядка малости, пренебрегая всеми членами более высоких порядков малости:
Обозначим величину дифференциального сопротивления нелинейного сопротивления Rd в точке i=I0:
Тогда подставляя (53), (52), (49) в (51) получим линеаризованное уравнение цепи:
где
Представляя решение (54) в виде
получим характеристическое уравнение системы
и явную зависимость о(t) в следующей форме:
Следовательно, состояние линейной системы (54) асимптотически устойчиво, а исходная нелинейная система (51) устойчива в обычном смысле. Если выполняется условие Используя обозначение (55) окончательно получаем:
Последнее выражение представляет собой известный критерий Кауфмана для устойчивой рабочей точки электрической дуги.
Устойчивость решений уравнения Дуффинга.
Запишем уравнение движения неконсервативного нелинейного осциллятора, находящегося под гармоническим внешним воздействием, для случая среды с вязким трением (7. 2), (11. 1)
Формулы с двойными номерами здесь – (7. 2), (11. 1) - и ниже – (7. 5), (3. 20), (9. 5), (11. 3), (11. 5) – цитируются по книге [4]. 2) Поскольку символ л использован везде в настоящем разделе для обозначения корней характеристических уравнений.
где символом д обозначена в соответствии с (7. 5) удельная вязкость среды; щ0, м – (3. 20), F – (9. 5). Правую часть уравнения можно представить в виде суммы синусной и косинусной компоненты:
где F1, F2 определяются выражениями (11. 3) и справедливы формулы (11. 5). При исследовании устойчивости для описания поведения рассматриваемой системы при появлении малых возмущений необходимо использовать полную подстановку Ван дер Поля:
где a(t), b(t) – медленно изменяющиеся функции. Вычислим первую и вторую производные функции y(t) по времени t:
Используя медленность изменения функции a(t), b(t) и малость параметра д, пренебрежем в формулах (64) слагаемыми вторых порядков малости: Подставив последние выражения и (63) в уравнение (61), получим: Тригонометрический двучлен третьей степени в левой части равенства без учета всех компонент, кроме колебаний с основной частотойщ может быть представлен в следующем виде: Подставив это выражение в предыдущие и сгруппировав слагаемые с одинаковыми тригонометрическими функциями получим два соотношения:
Отсюда, разрешая равенства относительно a, b, можно записать систему “укороченных” уравнений: Рассмотрим стационарное решение:
Тогда для определения амплитуд стационарных колебаний a0, b0 на основании системы (69) получаем алгебраические уравнения:
Ранее решение (11. 6) было получено в частном случае наличия одного только синусного колебания:
При этом из (71) получаем выражения:
совпадающие, как и следовало ожидать, с (11. 9). Тогда. Использую (11. 5), (73) можно записать формулу (11. 10):
Определяющую резонансную зависимость рассматриваемого осциллятора |a0|(щ) или его “управляющую” характеристику |a0|(F). Для исследования устойчивости полученных стационарных колебаний с амплитудами a0, b0 (70), (71) введем теперь в рассмотрение их малые возмущения о(t), з(t):
где
Тогда подставляя (75) в укороченные уравнения (69) и используя условия для стационарных амплитуд (71), получим нелинейную систему возмущенного движения в следующем виде:
Анализируем уравнения (77), используя малость возмущений (76): где постоянные коэффициенты aнm для частного случая, рассмотренного ранее (72), определяются так:
Записывая возмущения в экспоненциальной форме (13), получая систему алгебраических уравнений (14), приравнивая нулю определитель этой системы (15), (16), окончательно получаем характеристическое уравнение в виде: или
Следовательно, необходимые и достаточные условия устойчивости линеаризованной системы можно записать так:
или, используя обозначения (79)
Первое из этих условий в рамках рассматриваемой задачи, очевидно, выполняется всегда. Второе условие требует более детального анализа. Чтобы его осуществить, определим зависимостьF(a0) на основании соотношения (74):
и вычислим производную dF/da0:
Последнее выражение легко преобразовать к виду:
Если теперь в этой зависимости амплитуду вынужденных колебаний a0 (11. 6) заменить на модуль этой величины |a0|, как это обычно делается при построении резонансных характеристик |a0|(щ) осцилляторов и их управляющих характеристик |a0|(F), то никаких изменений в соотношении (86) не произойдет. В частности, знак производной не изменится. Таким образом, коэффициент при фигурной скобке в (86) является величиной существенно положительной.
Тогда, сравнивая выражение, заключенное в фигурные скобки в (86) с условием (83) видим, что условие устойчивости идентично условию положительности производной.
или условию положительности обратной величины
Следовательно, все точки управляющей характеристики осциллятора с положительным наклоном касательной соответствует устойчивым режимам колебаний. На рисунках? ?? ?? ?? ?? ?? ?? ? и 12. 3 эти ветви изображены сплошными кривыми. “Падающей” ветви характеристики соответствуют неустойчивые колебания. На рис? ?? ?? ?? ?? ? Эта ветвь представлена пунктиром, а на рис. 12. 3 вообще отсутствует.
В переходных точках перегиба кривой
имеет место упомянутый ранее “критический” случай, поскольку можно показать, что появляются корни характеристического уравнения с равными нулю вещественными частями. Анализа устойчивости на основе линейного приближения здесь оказывается недостаточно. Устойчивость или неустойчивость в этих точках определяют слагаемые высших порядков малости в уравнениях возмущенного движения.
