В

Рис. 2.2. Звено САР

большинстве случаев можно линеаризовать нелинейные зависимости, используя метод малых отклонений или вариаций. Для рассмотрения его обратимся к некоторому звену системы автоматического регулирования (рис. 2.2). Входная и выходная величины обозначены через X 1 иX 2 , а внешнее возмущение – через F(t).

Допустим, что звено описывается некоторым нелинейным дифференциальным уравнением вида

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

Основанием для линеаризации служит предположение о достаточной малости отклонений всех переменных, входящих в уравнение динамики звена, так как именно на достаточно малом участке криволинейную характеристику можно заменить отрезком прямой. Отклонения переменных отсчитываются при этом от их значений в установившемся процессе или в определенном равновесном состоянии системы. Пусть, например, установившийся процесс характеризуется постоянным значением переменной Х 1 , которое обозначим Х 10 . В процессе регулирования (рис. 2.3) переменная Х 1 будет иметь зна­чениягде
обозначает отклонение переменнойX 1 от установившегося значения Х 10 .

А

Рис. 2.3. Процесс регулирования в звене

налогичные соотношения вводятся для других переменных. Для рассматриваемого случая имеем: а также
.

Далее можно записать:
;
и
, так как
и

Все отклонения предполагаются достаточно малыми. Это математическое предположение не противоречит физическому смыслу задачи, так как сама идея автоматического регулирования требует, чтобы все отклонения регулируемой величины в процессе регулирования были достаточно малыми.

Установившееся состояние звена определяется значениями Х 10 , Х 20 и F 0 . Тогда уравнение (2.1) может быть записано для установившего состояния в виде

Разложим левую часть уравнения (2.1) в ряд Тейлора

где  – члены высшего порядка. Индекс 0 при частных производных означает, что после взятия производной в её выражение надо подставить установившееся значение всех переменных
.

В состав членов высшего порядка в формуле (2.3) входят высшие частные производные, умноженные на квадраты, кубы и более высокие степени отклонений, а также произведения отклонений. Они будут малыми высшего порядка по сравнению с самими отклонениями, которые являются малыми первого порядка.

Уравнение (2.3) является уравнением динамики звена, так же как (2.1), но записано в другой форме. Отбросим в этом уравнении малые высшего порядка, после чего из уравнения (2.3) вычтем уравнения установившегося состояния (2.2). В результате получим следующее приближённое уравнение динамики звена в малых отклонениях:

В это уравнение все переменные и их производные входят линейно, то есть в первой степени. Все частные производные представляют собой некоторые постоянные коэффициенты в том случае, если исследуется система с постоянными параметрами. Если же система имеет переменные параметры, то уравнение (2.4) будет иметь переменные коэффициенты. Рассмотрим только случай постоянных коэффициентов.

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

где введены следующие обозначения

. (2.6)

Кроме того, для удобства принято все дифференциальные уравнения записывать в операторной форме с обозначениями

Тогда дифференциальное уравнение (2.5) запишется в виде

Эту запись будем называть стандартной формой записи уравнения динамики звена.

Коэффициенты Т 1 и Т 2 имеют размерность времени – секунды. Это вытекает из того, что все слагаемые в уравнении (2.8) должны иметь одинаковую размерность, а например, размерность(илиpx 2) отличается от размерности х 2 на секунду в минус первой степени (
). Поэтому коэффициенты Т 1 и Т 2 называютпостоянными времени .

Коэффициент k 1 имеет размерность выходной величины, деленную на размерность входной. Он называетсякоэффициентом передачи звена. Для звеньев, у которых выходная и входная величины имеют одинаковую размерность, используются также следующие термины: коэффициент усиления – для звена, представляющего собой усилитель или имеющего в своем составе усилитель; передаточное число – для редукторов, делителей напряжения, масштабирующих устройств и т. п.

Коэффициент передачи характеризует статические свойства звена, так как в установившемся состоянии
. Следовательно, он определяет крутизну статической характеристики при малых отклонениях. Если изобразить всю реальную статическую характеристику звена
, то линеаризация дает
или
. Коэффициент передачи k 1 будет представлять собой тангенс угла наклона касательной в той точкеC(см. рис. 2.3), от которой отсчитываются малые отклонения х 1 и х 2 .

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

Кроме того, отсюда вытекает другой, графический способ линеаризации. Если известна статическая характеристика и точка C, определяющая установившееся состояние, около которого происходит процесс регулирования, то коэффициент передачи в уравнении звена определяется графически из чертежа по зависимости k 1 = tgcучетом масштабов чертежа и размерностиx 2 . Во многих случаяхграфический метод линеаризации оказывается более удобным и быстрее приводит к цели.

Размерность коэффициента k 2 равна размерности коэффициента передачи k 1 , умноженной на время. Поэтому часто уравнение (2.8) записывают в виде

где
– постоянная времени.

П

Рис. 2.4. Двигатель независимого возбуждения

остоянные времени Т 1 , Т 2 и Т 3 определяют динамические свойства звена. Этот вопрос будет рассмотрен подробно ниже.

Коэффициент k 3 представляет собой коэффициент передачи по внешнему возмущению.

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

Для нахождения дифференциального уравнения, связывающего приращение скорости с приращением напряжения на обмотке возбуждения, запишем закон равновесия электродвижущих сил (эдс) в цепи возбуждения, закон равновесия эдс в цепи якоря и закон равновесия моментов на валу двигателя:

;

.

Во втором уравнении для упрощения опущен член, соответствующий эдс самоиндукции в цепи якоря.

В этих формулах R В и R Я – сопротивления цепи возбуждения и цепи якоря; І В и І Я – токи в этих цепях; U В и U Я – напряжения, приложенные к этим цепям; В – число витков обмотки возбуждения; Ф – магнитный поток; Ω – угловая скорость вращения вала двигателя; М – момент сопротивления от внешних сил;J– приведенный момент инерции двигателя; С Е и С М – коэффициенты пропорциональности.

Допустим, что до появления приращения напряжения, приложенного к обмотке возбуждения, существовал установившийся режим, для которого уравнения (2.10) запишутся следующим образом:

(2.11)

Если теперь напряжение возбуждения получит приращение U В = U В0 + ΔU В, то все переменные, определяющие состояние системы, также получат приращения. В результате будем иметь: І В = І В0 + ΔІ В; Ф = Ф 0 + ΔФ; I Я = I Я0 + ΔІ Я; Ω = Ω 0 + ΔΩ.

Подставляем эти значения в (2.10), отбрасываем малые высшего порядка и получаем:

(2.12)

Вычитая из уравнений (2.12) уравнения (2.11), получим систему уравнений для отклонений:

(2.13)

В

Рис. 2.5. Кривая намагничивания

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

Совместное решение системы (2.13) даёт

где коэффициент передачи, ,

; (2.15)

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

(2.16)

где L B = a B – динамический коэффициент самоиндукции цепи возбуждения; электромагнитная постоянная времени двигателя, с,

. (2.17)

Из выражений (2.15) – (2.17) видно, что рассматриваемая система является по существу нелинейной, так как коэффициент передачи и «постоянные» времени, на самом деле – не постоянны. Их можно считать постоянными только приближенно для какого-то определенного режима при условии малости отклонений всех переменных от установившихся значений.

Интересным является частный случай, когда в установившемся режиме U B0 = 0; І B0 = 0; Ф 0 = 0 и Ω 0 = 0. Тогда формула (2.14) приобретает вид

. (2.18)

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

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

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

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

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

Подставляя сюда вместо случайной функции какое-либо ее каноническое разложение, получим:

Это равенство представляет случайную функцию как некоторую, вообще говоря нелинейную, функцию случайных величин в которую аргумент 5 входит как параметр:

Линеаризуя эту функцию обычным в теории вероятностей способом (см. § 31) и принимая во внимтние, что математические ожидания величин равны нулю, будем иметь:

есть значение производной функции по случайной величине при нулевых значениях всех величин что и отмечено нуликом внизу у квадратной скобки. Формула (100.5) дает приближенное каноническое разложение случайной функции с коэффициентами разложения и координатными функциями

Принимая во внимание, что математические ожидания всех величин равны нулю, получим из (100.5) следующую приближенную формулу для математического ожидания случайной функции

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

Корреляционная функция случайной функции на основании общей формулы (56.2), выразится приближенной формулой

Большинство реальных систем нелинейны, т.е. поведение системы описывается уравнениями:

Часто на практике нелинейные системы можно аппроксимировать линейной в некоторой ограниченной области.

Предположим, что
для уравнения (1) известно. Заменим систему (1,2) подставив начальные условия

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

Выход
найдем в результате решения возмущенных уравнений.

Разложим правую часть в ряд Тейлора.

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

Вычитая исходное решение из разложений, получаем следующие линеаризованные уравнения:

.

Частные производные обозначим как коэффициенты зависящие от времени

Эти выражения можно переписать в виде

Получим линеаризованные уравнения в точках равновесия
.

. В точке

Решение этого уравнения

Продифференцируем правую часть исходного уравнения по x , получим

.

Выполним линеаризацию уравнения для произвольного начального значения
.

Получаем линеаризованную систему в виде нестационарного уравнения

Решение линеаризованной системы имеет вид:

.

1.7. Типовые возмущающие воздействия

Внешние возмущающие воздействия могут иметь различный характер:

мгновенного действия виде импульса и постоянного действия.

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

(t)- функция при интегрировании обладает следующими фильтрующими свойствами:

Интегрируемое произведение произвольной функции
и(t)- функции отфильтровывает из всех значений
только то, которое соответствует моменту приложение мгновенного единичного импульса.

Линейное возмущение

Гармоническое возмущение

2 U. Системы второго порядка

2.1.Приведение уравнений второго порядка к системам уравнений первого порядка

Пример линейной стационарной системы.

Другое описание этой же системы второго порядка дается парой связанных дифференциальных уравнений первого порядка

(2)

где связь между коэффициентами этих уравнений определяется следующими соотношениями

2.2. Решение уравнений второго порядка

Применяя дифференциальный оператор
уравнение можно представить в более компактном виде

Решается уравнение (1) в 3 этапа:

1) находим общее решение однородного уравнения;

2) находим частное решение ;

3) полное решение есть сумма этих двух решений
.

Рассматриваем однородное уравнение

будем искать решение в форме

(5)

где
действительная или комплексная величина. При подстановке (5) в (4) получаем

(6)

Это выражение является решением однородного уравнения, если s удовлетворяет характеристическому уравнению

При s 1  s 2 решение однородного уравнения имеет вид

Тогда ищем решение в виде
и подставляя его в исходное уравнение

Откуда следует, что
.

Если выбрать

. (8)

Частное решение исходного уравнения (1) ищем методом вариации
в форме

исходя из (11), (13) получаем систему

Полное решение уравнения.

Заменой переменных получим уравнение второго порядка:

      ФАЗОВАЯ ПЛОСКОСТЬ

Двумерным пространственным состоянием или фазовой плоскостью называется плоскость, в которой две переменные состояния рассматриваются в прямоугольной системе координат

- эти переменные состояния образуют вектор
.

График изменения
образует траекторию движения. Необходимо указать направление движения траектории.

Состояние равновесия называется такое состояние , в котором система остается при условии, что
Состояние равновесия можно определить (если оно существует) из соотношений

при любом t .

Состояния равновесия иногда называются критическими, основными или нулевыми точками.

Траектории системы не могут пересекаться друг с другом в пространстве, что вытекает и единственности решения дифференциального уравнения.

Ни одна траектория не проходит через состояние равновесия хотя и могут сколь угодно близко приближаться к особым точкам (при
) .

Типы точек

1 Регулярная точка есть любая точка, через которую может проходить траектория, точка равновесия не является регулярной.

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

Рассмотрим систему

Для определения состояния равновесия решим следующую систему уравнений

.

Получаем зависимость между переменными состояния
.

любая точка которой есть состояние равновесия. Эти точки не является изолированными.

Заметим, что для линейной стационарной системы

начальное состояние оказывается состоянием равновесия и изолированным, если детерминант матрицы коэффициентов
, тогда
есть состояние равновесия.

Для нелинейной системы второго порядка состояние равновесия называется простым , если соответствующая матрица Якоби не равна 0.

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

Рассмотрим решение уравнения состояния для линейной системи второго порядка:
.

Эту систему можно представить двумя уравнениями первого порядка,

обозначим
,

Характеристическое уравнение
и решение будет следующим:

Решение уравнения записывается в виде

25. Методы линеаризации нелинейных САУ.

С т. зрения передачи и преобразования сигнала НЛ отлич. от линейных систем тем, что мгновенный коэфффициент передачи зависит от значения входного сигнала. САУ, содержащие звенья, динамика которых описывается НЛ дифференц. уравнениями относят к НЛ системам .

НС-динамика к-х описывается нелин-ми диф ур-ми, это сис-мы, имеющие нелинейную стст-ю хар-ку.

Систему можно представить в виде соединения из 2-х элементов:

можно свести к:

ЛЧ

ЛЧ описывается обычными диф ур-ми с пост-ми коэфф-ми.

НЭ является безинерционным и его выходная величина и вход. величина связаны связаны между собой НЛ алгебраическим уравнением. Нелинейность обусловлена нелинейностью статической характеристики одного из элементов системы.

Методы линеаризации нелинейных САУ.

метод гармонической линеаризации

статическая линеаризация

совместная стат и гармон линеаризация

вибролинеаризация

Метод гармонической линеаризации.

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

В рез-те линеаризации нелин стат хар-ку заменяют эквивалентным линейным звеном с коэффициентами

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

k / Г и k Г – гарм-ие коэф-ты усиления нелин звена.

Отличия гарм-ой линеар-ии от обычной:

При гарм-ой линеаризации нелин хар-ку заменят прямой, крутизна которой зависит от амплитуды входного сигнала.

Гарм-ая линеаризация позволяет вместо нелин звена получить линейное, к-т усиления которого зависит от а.

Гарм-ая линеар-ия дает возможность опредилить св-ва нелин САУ методами линейной теории автом-х сис-м.

Статическая линеаризация.

Этот метод приближенного исследования точности нелин сис в стационарных случ реж-ах.

В качестве примера возьмем нелин звено со стат хар-ой типа насыщение.

Пусть на входе стационарный случ. Сигнал.

X (t )= m x + x 0 (t )

Y(t)=m y +y 0 (t)

Задача стат лин-ии закл-ся в том чтобы найти линейное звено дающее при том же вх сигнале x (t ) вых сигнал = эквивалентному вых сигналу нелин звена при этом надо чтобы эквив-й сигнал максимально приближался к y (t ).

Точность линеариз зависит от того, какой критерий выбран для сравнения y экв и y .

Сущ 2 критерия сравнения y экв и y :

1. по первому способу линеаризация осущ-ся исходя из след условий

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

2. основан на линеаризации разности

К-ты стат линеар-ии:

Совместная статическая и гармоническая линеаризация.

Случай когда в сис присутствуют автоколебания и на вх сис подаются случ воздействия:

f(t)=m f +f 0 (t)

x(t)=m x +x 0 (t)+a*sin w a t

Из-за неприменимости принципа суперпозиции необходимо учитывать наличие всех 3-х составляющих для этого надо осущ-ть совместную стат и гарм линеа-ию, в рез-те этого сигнал на выходе:

в случ симметр-ой нелин стат хар-ки пост состав-ую

m y = y 0 = k сг0 m x

эти 4 к-та опред-ся по фор-ам для гарм-ой и стат линеар-ии. Эти к-ты уже будут зависеть от 4-х составляющих (m x , s x , a , w a )

При исследовании сис m x , s x , a , w a - определяются совместным решением ур-ий для колебательной составляющей и для случ состав-ей.

Применяя совместно стат и гармонич линеаризацию можно решать две задачи:

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

можно исследовать точность сис в случ режимах при наличии сис гармонических колебаний.

Вибролинеаризация.

Испол-ся для исключения эффекта наличия нелин-х хар-к (люфт и зона нечувст-ти).

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

Рассмотрим метод виб-ой лин-ии на примере релейной сис:

зависимость y 0 = f (x 0 ) ,где y 0 зависит от x 0 и от формы нелин-ой стат хар-ки, т.о. при наличии переменного воздействия, этот элемент пропускает пост воздействие x 0 как звено непрерывного действия.

Сам процесс виб-й лин-ии можно трактовать как процесс модуляции, в данном примере реле явл-ся модулятором высокочас-ое воздействие - сигнал несущей частоты, а НЧ вх сигнал x 0 явл-ся модулирующим сигналом. В данном случае осущ-ся ШИМ и ф-ей модулир-го сигнала явл-ся ширина вых имп-са и условие неискаженной передачи НЧ-составляющей явл-ся f ВЧ / f НЧ >=3

Когда реле работает в составе САУ обычно НЧ сигнал x 0 представляет собой сигнал управления и изменения во времени x 0 и есть перех-ой процесс в сис.

ВЧ воздействие осущ виб-ой лин-ей м.б. получено 3-я способами:

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

Путем создания автоколебаний в самой САУ.

Путем создания скользящего режима.

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

Нелинейная система представляется в виде соединения линейной и нелинейной части (рис. 5).

Рис. 5 Схема нелинейной системы

Выходной сигнал нелинейной части системы в общем случае определяется выражением

Обозначим как передаточную функцию линейной части. Система уравнений примет вид

Найдем условия, при которых на выходе линейной части системы возникают гармонические колебания вида

В этом случае сигнал y(t) нелинейной части будет представлять собой также периодическую функцию, но отличную от синусоиды. Эту функцию можно разложить в ряд Фурье

В этом выражении a i и b i - коэффициенты Фурье. Для симметричных нелинейностей F 0 =0.

Основным условием, которое накладывает метод на линейную часть системы, является условие фильтра нижних частот. Считается, что линейная часть пропускает только первую гармонику колебаний. Данное допущение позволяет считать высшие гармоники в (7.19) несущественными и ограничиться рассмотрением только первой гармоники сигнала y(t).

то выражение (7.20) можно переписать в виде

Первое уравнение системы (7.17) примет вид

В этом выражении


Результат замены нелинейности F(x,sx) выражением

и называется гармонической линеаризацией. Величины q и q 1 называются коэффициентами гармонической линеаризации или просто гармоническими коэффициентами. Для однозначных нелинейностей обычно q 1 =0 . Формулы для гармонических коэффициентов, соответствующих типовым нелинейностям, приводятся в приложениях.

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

Рассмотрим методику определения амплитуды и частоты автоколебаний.

1). В характеристическом уравнении системы, полученном из (7.22) делаем замену s=j и получим

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

  • 3).Решение этой системы дает частоту и значения гармонических коэффициентов. Если эти значения вещественны и положительны, то в системе существует предельный цикл. По значениям гармонических коэффициентов можно определить амплитуду предельного цикла.
  • 4). Общим признаком устойчивости предельного цикла, т.е. существования автоколебаний, является равенство нулю предпоследнего определителя Гурвица при полученных значениях амплитуды и частоты предельного цикла. Часто более удобно использовать условие устойчивости предельного цикла, в основе которого лежит критерий устойчивости Михайлова.

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

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

Рис.6 Схема нелинейного рулевого привода

Зададим следующие параметры нелинейности скоростной характеристикм рулевого привода: b = 0.12, k 1 = tg =c/b = 6.7. Коэффициенты гармонической линеаризации этой нелинейности определяются выражениями

Заменив в схеме нелинейную характеристику гармоническим коэффициентом, получим передаточную функцию рулевого привода

Подставим эту передаточную функцию в структурную схему системы стабилизации угла тангажа и определим передаточную функцию замкнутой системы

В характеристическом уравнении замкнутой системы сделаем замену s = j и выделим вещественную и мнимую части.

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

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

По этим значениям можно вычислить для двух случаев все коэффициенты характеристического уравнения и определить частоты, соответствующие каждому значению q(А). Получим:

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

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

Производная от коэффициента гармонической линеаризации, входящая в полученные выражения, вычисляется по формуле


Расчеты по выше приведенным формулам показывают, что первый предельный цикл не устойчив и возникает он при (0) 0.1166(6.7 0 ). Если начальное отклонение меньше указанного, то процесс на входе нелинейного элемента затухает (рис.7. 7) и система устойчива.


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


Рис. 8

Путем моделирования определено, что область притяжения устойчивого предельного цикла лежит приблизительно в пределах (0) 0.1167 - 1.4 (6.71 0 - 80.2 0 ).