Устойчивость установившихся движений

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

Уравнения движения системы запишем в форме уравнений Лагранжа II рода

d/dt (∂T/∂q)T — (∂T/∂q)T = — (∂П/∂q)T + QН. (4.1)

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

Устойчивость относительного равновесия. Измененная потенциальная энергия W1

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

Предположим, что механическая система, имеющая n степеней свободы, имеет r позиционных (r < n) и (n−r) циклических координат. Вектор позиционных координат системы обозначим через q = (q1,…,qr)T, а вектор циклических координат − через φ = (φ1,…,φn−r)T, где для удобства принято qr+11, qr+22,…, qnn−r.

Тогда кинетическую энергию можно записать в виде

T = ½ ∑ aijqiqj = ½q•TA(q)q + q•TB(q)φ + ½φ•TC(q)φ (4.20)

Элементы матриц A=A(q), B=B(q) и C=C(q), а также потенциальная энергия П = П(q) зависят только от позиционных координат q.

Пусть обобщенных неконсервативных сил Qн = Qн(q) нет, а по циклическим координатам действуют обобщенные управляющие силы U.

Тогда уравнения движения данной механической системы в форме уравнений Лагранжа II рода имеют следующий вид

d/dt (∂T/∂q)T — (∂T/∂q)T = — (∂П/∂q)T (4.21)

d/dt (∂T/∂φ) = U

Или в скалярном виде:

d/dt ∂T/∂qi — ∂T/∂qi = — ∂П/∂qi (i = 1,…,r)

d/dt (∂T/∂φj) = Uj (j = 1,…, n − r)

Вторая группа уравнений (4.21) после подстановки (4.20) имеет вид

•• + dС/dt φ + d/dt [BTq] = U . (4.22)

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

U0 = d/dt [BTq] dС/dt φ (4.23)

При таком управлении U0 уравнение (4.22) принимает наиболее простой вид

••= 0,

из которого следует первый (векторный) интеграл движения

φ = ω = const, (4.24)

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

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

Пример: Если в системе (4.21) все циклические координаты соответствуют углам поворота роторов, входящих в систему, то условия

φ = ω = const

означают, что все роторы, входящие в систему, управляются двигателями неограниченной мощности и вращаются равномерно с постоянными угловыми скоростями ωr+1, ωr+2, … , ωn независимо от величин сил, действующих в системе.

При значении U0 управления (4.23) число степеней свободы механической системы (4.21) уменьшается, поскольку движение по циклическим координатам задано

φ = ωt + φ0 = const. (4.24)

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

Кинетическая энергия (4.20) в данном случае зависит только от позиционных координат и скоростей q, q и параметров ω и ее можно представить в виде суммы трех однородных относительно позиционных скоростей q функций

T (q, q) = T2 (q, q) + T1 (q, q) + T0 (q), (4.25)

где индекс означает степень однородности функции,

T2 (q, q) = ½qTA(q)q , T1 (q, q) = qTB(q)ω , T0 (q) = ½ ωTC(q)ω.

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

d/dt (∂T2/∂q)T — (∂T2/∂q)T = — (∂W1/∂q)T) + Г1 (4.26)

Или в скалярном виде так:

d/dt ∂T2/∂qi — ∂T2/∂qi =
— ∂W1/∂qi + Г1 (i = 1,…,r) (4.26′)

Определение 4.3:
Функция W1(q)=П(q) − T0(q) называется измененной потенциальной энергией системы.

Вектор

Г1(q, q) = d/dt (∂T1/∂q)T — (∂T1/∂q)T

с элементами

Г1i(q, q) = d/dt ∂T1/∂qi — ∂T1/∂qi = ∑gtjqj

представляет собой обобщенные гироскопические силы, коэффициенты которых образуют кососимметрическую матрицу G = (gij) = — GT с элементами

gij (q, q) = — gji (q, q) = ∑ [∂bkj/∂qi — ∂bki/∂qj]ωk,

где bkj − элементы матрицы B.

Гироскопические силы могут отсутствовать в уравнениях (4.26), если все коэффициенты gij=0. В этом случае система (4.26) называется гироскопически несвязанной, в
противном случае − гироскопически связанной.

Умножим обе части уравнений (4.26) на q•T (или уравнения (4.26′) на q j и сложим их).

Согласно основному свойству гироскопических сил сумма работ гироскопических сил на любом действительном перемещении системы равна нулю, тогда уравнение баланса энергии примет вид

dH1/dt = 0, (4.27)

где H1 (q, q) = ∑ ∂T2/∂qi — T2 + W1 = T2 + W1. При вычислении функции H1 использована теорема Эйлера об однородных функциях.

Из уравнений (4.27) следует первый интеграл − обобщенный интеграл энергии (интеграл Якоби)

H1 = T2 + W1= const.

Уравнения (4.26) имеют частное решение

qp(t) = qp = const, qp(t) ≡ 0 , (4.28)

описывающее положение относительного равновесия системы.

Подставим (4.28) в (4.26) и найдем уравнения относительного равновесия

0 = — ∂W1/∂q . (4.29)

Именно решением этих уравнений (4.29) является решение (4.28)

Определение 4.4. Равновесие (4.28) системы (4.26) с регулируемыми циклическими координатами (4.24), при котором постоянны все позиционные координаты qp =const и
циклические скорости φ=ω=const называется положением относительного равновесия системы.

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

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

Для этого согласно теореме Ляпунова необходимо найти знакоопределенную функцию Ляпунова, производная которой в силу уравнений возмущенного движения будет функцией знакопостоянной противоположного знака или тождественно равной нулю.

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

U = H1 =T2 + W1. (4.30)

Поскольку согласно (4.25) выражение квадратичной формы кинетической энергии T2 представляет собой определенно положительную функцию позиционных скоростей q, то mдостаточные условия устойчивости относительного равновесия (4.28) могут быть найдены с помощью следующей модифицированной теоремы Лагранжа.

Теорема 4.2 (модифицированная теорема Лагранжа):
Если в положении относительного равновесия (4.28) измененная потенциальная энергия системы W1(q) имеет изолированный минимум, то такое положение относительного равновесия устойчиво относительно позиционных координат q и скоростей q.

Доказательство аналогично доказательству теоремы 4.1.

Замечание: В данном случае измененная потенциальная энергия системы

W1(q) = П(q) − T0 = П(q) − ½ ωTC(q)ω = W1(q, ω)

зависит от параметров ω=const и при исследовании устойчивости относительных равновесий qp=qp(ω) (4.28) также имеет место закон смены устойчивости в точках бифуркации.

Пример: Устойчивость относительного равновесия тяжелого колечка A на вращающейся с постоянной скоростью окружности (см.рисунок).

Задача:
Тяжелое колечко A массой m может скользить без трения по окружности радиуса OA=R, которая вращается вокруг своего вертикального диаметра. Угловая скорость вращения окружности постоянна φ = ω, φ − угол поворота окружности. Найти относительные равновесия колечка на окружности и исследовать их устойчивость. На плоскости параметров системы θ, ω построить бифуркационную диаграмму θ = θ (ω), где θ − угол, определяющий положение колечка на окружности.




Оставить комментарий





Статистика

Рейтинг@Mail.ru