Постановка задачи об устойчивости движения. Определение устойчивости по Ляпунову


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


(1.1)

Здесь T =1/2q•TA(q)q, П = П(q) — кинетическая и потенциальная энергия системы, q = (q1, …, qn)T, q = (q1, …, qn)T — n-мерные векторы-столбцы ее обобщенных координат qi и скоростей
qi, QН = QН(q,q) = (QН1, …, QНn)T — вектор-столбец обобщенных неконсервативных сил.

После дифференцирования уравнения (1.1) можно записать в виде
A(q)q•• + b(q,q)+ c(q) = QН(q,q), (1.2)
где b(q,q), c(q) – некоторые нелинейные вектор-функции, которые получаются после дифференцирования кинетической и потенциальной энергий.

Пусть уравнения (1.1) и, соответственно, уравнения (1.2) имеют частное решение
qp = ξ(t), qp = ξ(t) , (1.3)
удовлетворяющее начальным условиям
qp(t0) = ξ(t0) = ξ0, qp(t0) = ξ(t0) =ξ0. (1.4)
Если решение (1.3) подставить в уравнения (1.2), то получим тождество

A(ξ)ξ•• + b(ξ,ξ)+ c(ξ) = QН(ξ,ξ).
Решение (1.3) — (1.4) описывает некоторое движение системы (положение равновесия, стационарное движение, программное движение и т.д.), которое нас интересует и которое
мы хотим исследовать.

Сделаем постановку задачи об исследовании устойчивости движения (1.3) — (1.4).

Следуя А.М. Ляпунову введем следующие определения.

Определение 1.1. Движение (1.3) — (1.4) системы (1.1), устойчивость которого исследуется, называется невозмущенным движением.

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

Для исследования устойчивости невозмущенное движение (1.3) — (1.4) необходимо сравнить с другими движениями (возмущенными) системы (1.1), которые запишем в виде
q = q (t), q = q(t), (1.5)
существующими у системы (1.1) при тех же обобщенных силах Q, но удовлетворяющими другим начальным условиям, которые можно записать в виде
q(t0) = ξ0+ x10, q(t0) = ξ0+ x20. (1.6)

Величины x10, x20 — некоторые вещественные постоянные, которые называются начальными возмущениями (их еще называют мгновенными возмущениями).

Будем считать, что начальные возмущения x10, x20 численно достаточно малы.

Определение 1.2. Движения (1.5) — (1.6) системы (1.1), с которыми сравнивается невозмущенное движение, называются возмущенными движениями.

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

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

Интуитивно ясно, что, если любое возмущенное движение, близкое к невозмущенному движению в начальный момент времени, будет оставаться близким к нему и в остальные моменты времени, то невозмущенное движение можно считать устойчивым.

Если же «расстояние» от возмущенного движения до невозмущенного движения с течением времени будет увеличиваться и станет не малой величиной, то невозмущенное движение можно считать неустойчивым.

Введем точное определение устойчивости по Ляпунову, которое принято за основное в различных областях науки и техники.

Определение 1.3 (устойчивости по Ляпунову) Невозмущенное движение qp=ξ(t),
qp(t) системы (1.1) называется устойчивым по Ляпунову относительно величин q = (q1,…, qn)T, q = (q1, …, qn)T, если для любых чисел ε1, ε2 > 0, сколь бы малы они ни были,
существуют числа δ1, δ2 > 0, зависящие от ε1, ε2 и t0, такие, что для всех решений q(t), q(t) уравнений (1.1), начальные условия которых q(t0), q(t0) удовлетворяют неравенствам
||q(t0) — ξ0|| < δ1, ||q(t0) — ξ0|| < δ2, (1.7) при всех t ≥ t0 будут выполняться неравенства ||q(t) — ξ(t)|| < ε1, ||q(t) — ξ(t)|| < ε2. (1.8)
Здесь ||q(t) — ξ(t)|| — евклидова норма вектора x1(t), характеризующая расстояние между возмущенным и невозмущенным движениями системы.

Геометрическая интерпретация определения устойчивости по Ляпунову в фазовом пространстве q1,…, qn , q1,…, qn

На рис.1.2 представлена геометрическая интерпретация определения устойчивости по Ляпунову на фазовой плоскости q, q для случая системы с одной степенью свободы (n=1). Начав при t=t0 движение из точки М0(q(t0), q(t0)), лежащей в Bδ-окрестности начальной точки ξ0, ξ0 невозмущенного движения ξ(t), ξ(t), где Bδ: |q(t0) — ξ0| < δ1, |q(t0) — ξ0| < δ2, изображающая точка М в случае устойчивости невозмущенного движения ξ(t), ξ(t) при t ≥ t0 не выйдет из Bε-окрестности данного невозмущенного движения, где Bε: |q(t) — ξ(t)|<ε1,
|q(t)-ξ(t)| < ε2.



Рис.1.2 Геометрическая интерпретация определения устойчивости по Ляпунову на фазовой плоскости


Комментарии запрещены.




Статистика