УДК 531.36+517.977
Вестник СПбГУ. Сер. 1. Т. 1 (59). 2014. Вып. 2
ОБ УПРАВЛЯЕМОСТИ СИСТЕМЫ А. Ю. ИШЛИНСКОГО*
А. С. Кулешов, В. В. Рыбин
Московский государственный университет им. М.В.Ломоносова, Российская Федерация, 119991, Москва, ГСП-1, Ленинские горы, 1
В 1965 году А. Ю. Ишлинский привёл пример неголономной механической системы малой размерности, которая не является системой Чаплыгина. Данная система состоит из трёх шероховатых цилиндров, два из которых, одинакового радиуса а, катаются без проскальзывания по неподвижной горизонтальной плоскости. Третий цилиндр, радиуса R, катается без проскальзывания по первым двум цилиндрам. В работе рассматриваются вопросы управляемости системы Ишлинского. При помощи теоремы Чжоу—Рашевского доказана управляемость данной системы. Библиогр. 3 назв. Ил. 1.
Ключевые слова: неголономная система, качение без проскальзывания, управляемость.
1. Введение. Рассматривается задача управления неголономной механической системой, предолженной А. Ю. Ишлинским [1]. При помощи теоремы Чжоу—Рашевского [2] доказана управляемость данной системы. Следует отметить, что ранее вопрос об управляемости системы Ишлинского рассматривался в работе [3]. К сожалению, не все результаты, полученные в [3], оказываются верными. В нашей работе мы уточняем и усиливаем результаты, полученные в работе [3].
Ниже сформулированы некоторые факты из теории управления, которыми мы будем пользоваться в дальнейшем.
2. Основные определения. Рассмотрим нелинейную управляемую систему с линейными управлениями
m
Uigi (x) , x e M, u =(ui,...,um) e Rm, (1)
i=1
где M — связное гладкое многообразие, gi — аналитические векторные поля на M, а ui (■) —измеримые локально ограниченные управления.
Для системы (1) можно сформулировать следующую задачу управления. Для заданных точек p, q e M и времени T > 0 требуется найти управление u(t), t e [0, T], переводящее систему (1) из точки p в точку q:
x (0)= p, x (T) = q.
Система называется вполне управляемой, если для любых p, q e M сформулированная задача разрешима для некоторого T > 0. Для линейной по управлениям системы (1) критерий полной управляемости даётся в терминах алгебры Ли векторных полей, порождённых полями в правой части системы:
Lie (gi,..., gm) = span (gi,..., gm, [gi, gj], [gi, [gj, gfc]],...).
Здесь span (...) —линейная оболочка векторных полей, а [gi, gj ] — скобка Ли векторных полей (коммутатор), которая в координатах имеет выражение
г т dgj dgi
__М = " -
* Работа выполнена при финансовой поддержке РФФИ (грант №13-01-00230).
Система (1) имеет полный ранг в точке x £ M, если
Lie (gi,..., gm) (x) = TxM.
Справедливо следующее утверждение:
Теорема (Чжоу—Рашевский [2]). Система (1) вполне управляема на M тогда и только тогда, когда она имеет полный ранг в любой точке x £ M.
Таким образом, вопрос управляемости для линейных по управлению систем (1) имеет полное решение. Ответ сводится к дифференцированию векторных полей и вычислению размерности их линейной оболочки.
3. Система А. Ю. Ишлинского. Постановка задачи. Рассмотрим механическую систему, состоящую из трёх шероховатых однородных цилиндров: два одинаковых цилиндра радиуса a катаются по горизонтальной плоскости, а третий цилиндр, радиуса R, катается по этим цилиндрам (рисунок).
Введём неподвижные системы координат Oxyz и Ox'y' z', оси Ox и Ox' которых направлены параллельно образующим нижних цилиндров, расположенных относительно друг друга под углом а. В силу того, что нижние цилиндры катаются по плоскости без проскальзывания, угол а является постоянным. Введём обобщённые координаты: x, y — координаты центра масс верхнего цилиндра (координата z = 2a + R остаётся постоянной), в — угол между образующей верхнего цилиндра и осью Ox, y — угол поворота верхнего цилиндра относительно своей оси, y i и y2 — соответствующие углы поворота нижних цилиндров. Обозначим через ri и Г2 радиусы-векторы, проведённые из центра масс верхнего цилиндра в точки соприкосновения цилиндров. В проекции на оси системы координат Oxyz эти радиусы-векторы имеют вид
r 1 = (ay 1 - y) ctg вех + (ayi - y) ey - Rez;
cos в sin в .
r2 = ~—77,-7 (a<ß2 + isma — y cosa) еж+——тт-г (atf 2 + xsma — у cosa) ey—Rez.
sin (в — a)
sin (в — a)
Система А. Ю. Ишлинского.
Скорости точек соприкосновения верхнего цилиндра с нижними в проекции на те же оси равны
vi = 2афiey, V2 = —2аф2 sin aex +2аф2 cos aey.
Скорость центра масс верхнего цилиндра и его угловая скорость равны соответственно
v = Xex + yey, ш = ф cos + ф sin + . Условия качения без проскальзывания верхнего цилиндра по нижним имеют вид
v + [ш X ri] = Vi, v + [ш X Г2] = V2. В проекциях на оси системы координат Oxyz эти условия имеют вид X — Дф sin в — (аф1 — y) в = 0,
y + Дф cos в + (аф1 — y) в ctg в — 2аф1 = 0,
X sin (в — а) — Дф sin в sin (в — а) + 2аф2 sin а sin (в — а) —
— (аф2 + X sin а — y cos а) в sin в = 0, (2)
y sin (в — а) + Дф cos в sin (в — а) — 2аф2 cos а sin (в — а) +
+ (аф2 + X sin а — y cos а) в cos в = 0.
Таким образом, рассматриваемая система имеет шесть обобщённых координат (x, y, в, ф1, ф2, ф) и две степени свободы.
4. Постановка задачи управления. Следуя работе [3], переобозначим обобщённые координаты, характеризующие положение системы, следующим образом:
X = 91, y = 92, в = 9з, ф1 = 94, ф2 = 95, ф = qe,
и, кроме того, введём обозначения
аф1 — y = aq4 — 92 = z,
аф2 + x sin а — аф1 cos а = aq5 + 91 sin а — aq4 cos а = w.
Разрешим уравнения связей (2) относительно обобщённых скоростей 91, 92, 9з, 95, а в качестве управляющих воздействий выберем две оставшиеся обобщённые скорости: 94 = U1 и 9e = U2. Тогда система уравнений (2) перепишется следующим образом:
2az sin 93 sin а
91 =-«i + Дм2 sm 93,
w sin 93 + z sin а cos 93
2aw sin 93
92 =-«i — Дм2 cos93,
w sin 93 + z sin а cos 93
2a sin а sin 93
® = -:-¡-:-/o\
w sin 93 + z sin а cos 93 (3)
<94 = U1,
sin 93
<75 = —,-
sin (93 — а)
9e = u2,
и будет иметь вид нелинейной системы с линейными управлениями (1). Векторные поля 01 и 02 имеют в данном случае вид
01
( 2аг эт д3 эт а \ ад эт д3 + г эт а соэ д3
2аад эт д3 ад эт д3 + г эт а соэ д3
2а эт а эт д3 ад эт д3 + г эт а соэ д3
1
вт дз
эт (дз - а)
\
0
02
I Д эт д3 \ —Д соэ д3 0 0
0 1
/
Таким образом, для системы Ишлинского можно сформулировать описанную в начале работы задачу управления. Докажем теперь, что система Ишлинского является полностью управляемой.
5. Управляемость системы Ишлинского. Для доказательства управляемости системы Ишлинского нам необходимо из векторных полей 01, 02 и их скобок Ли,
0з = [01, 02] —скобка Ли первого порядка;
04 = [02, [01,02]] = [02,03]
05 = [[01,02] ,01] = [03,01]
06 = ^ [02, [01,02]]] = [02,04]
07 = ^ [[01,02] ,01]] = [02, 05]
08 = [01, [02, [01,02]]] = [01,04]
09 = [0Ь [[01,02] ,01]] = [01, 05]
— скобки Ли второго порядка;
— скобки Ли третьего порядка,
выбрать такие, чтобы матрица, составленная из выбранных векторных полей, имела полный ранг (ранг 6) в каждой точке конфигурационного многообразия. Тогда по теореме Чжоу—Рашевского система Ишлинского будет вполне управляемой всюду на конфигурационном многообразии М.
Впервые вопрос управляемости системы Ишлинского рассматривался в работе [3]. Доказательство управляемости, предложенное в этой работе, также опиралось на теорему Чжоу—Рашевского. При этом отмечалось, что векторные поля, образующие матрицу полного ранга, можно выбрать, рассматривая поля 01, 02 и их скобки Ли до третьего порядка включительно.
Однако это утверждение не является справедливым. Непосредственной проверкой можно установить, что скобки Ли g4, ge, g7 и g8 выражаются через другие скобки Ли по следующим формулам:
2R sin а 6R2 sin2 а
9А =--:-;-:-9з, 96 = ---:-:-з9з,
w sin 93 + z sin а cos 93 (w sin q3 + z sin а cos q3)2
2R sin а 2aR (sin2 (q3 — а) — sin2 q3) sin а
97 =--:-;-:-95 + --:-:-,2 . ,-793, 98 = ~97-
w sin q3 + z sin а cos q3 (w sin q3 + z sin а cos q3)2 sin (q3 — а)
Следовательно, среди векторных полей gi, g2 и их скобок Ли до третьего порядка включительно независимыми являются лишь
gi, g2, g3, g5, g9.
Матрица, составленная из этих векторных полей, будет иметь ранг 5 и не будет матрицей полного ранга. Поэтому для доказательства управляемости системы Ишлинского необходимо рассмотрение скобок Ли четвёртого порядка:
g10 = [g2, ge] , gil = [g2,g7] , g12 = [g2, g8] , g13 = [g2, g9] , g14 = [gl, ge] , g15 = [gi,gr], gie = [gi,g8] , g17 = [gi,g9].
Для этих скобок непосредственной проверкой можно установить следующие соотношения:
24R3 sin3 а
9 ю = ---:-:-3g3,
(w sin q3 + z sin а cos q3)
6R2 sin2 а 12aR2 (sin2 (q3 — а) — sin2 q3) sin а
'w sin q3 + z sin а cos q3)2 (w sin q3 + z sin а cos q3)3 sin (q3 — а)
6Ra2 (sin3 q3 cos (q3 — а) — cos q3 sin3 (q3 — а) + sin а sin2 q3) sin2
а
913 = -:-'(—:—;—:-\3 . 2 /-Ñ-:-
(w sin q3 + z sin а cos q3) sin (q3 — а)
3R sin а
"g9,
w sin q3 + z sin а cos q3
4Ra2 (sin3 q3 cos (q3 — а) — cos q3 sin3 (q3 — а) + sin а sin2 q3) sin2 а
gis = --—:-:-,3 . o .-г-93-
(w sin q3 + z sin а cos q3) sin2 (q3 — а)
4aR (sin2 (q3 — а) — sin2 q3) sin а 2R sin а
95--:-;-:-g9-
(w sin q3 + z sin а cos q3)2 sin (q3 — а) w sin q3 + z sin а cos q3 g12 = —g11, g14 = —g11, gie = —g15.
Таким образом, среди скобок Ли четвёртого порядка независимой с векторными полями gi, g2, g3, g5, g9 является только скобка
gi7 = [gb g9] = [g1, [g1, [[g1, g2] , gi]]]
Окончательно можно сделать вывод, что система Ишлинского является управляемой на всём конфигурационном многообразии. Векторные поля
01, 02, 03 = [01,02], 05 = [[01,02] ,01],
09 = [0Ь [[01, 02] ,01]] , 017 = [01, [01, [[01, 02] ,01]]]
являются независимыми и составляют матрицу полного ранга — ранга 6. Тем самым, система Ишлинского является вполне управляемой по теореме Чжоу—Рашевского.
6. Обсуждение результатов. После того, как управляемость системы Ишлинского установлена, возникает вопрос о построении управления U1 (t), М2 (t), переводящего систему (3) за время T > 0 из начального состояния p в терминальное состояние q с любой наперёд заданой точностью. Однако задача построения управления для такого типа систем является весьма нетривиальной в силу анизотропии пространства состояний. Иными словами, для таких систем возможность перемещения неодинакова по различным направлениям. Величина смещения в направлении полей! 01, 02 за малое время t есть O (t), в направлении коммутатора 03 = [01,02] есть O (t2), в направлении 04 = [02, 0з] и 05 = [03,01] есть O (t3) и т. д. В дальнейшем мы надеемся полностью изучить задачу построения явного вида управления для системы Ишлин-ского.
Литература
1. Неймарк Ю.И., Фуфаев Н. А. Динамика неголономных систем. М.: Наука, 1967. 519 с.
2. Аграчёв А. А., Сачков Ю. Л. Геометрическая теория управления. М.: Физматлит, 2005. 392 с.
3. Jarzebowska E., McClamroch N. H. On nonlinear control of the Ishlinsky system as an example of a nonholonomic non-Chaplygin system // Proceedings of the American Control Conference. Chicago, Illinois, USA. June 2000. Vol. 5. P. 3249-3253.
Статья поступила в редакцию 26 декабря 2013 г.
Сведения об авторах
Кулешов Александр Сергеевич — кандидат физико-математических наук, доцент; kuleshov@mech.math.msu.su
Рыбин Вадим Владимирович — аспирант; RybinVV@gmail.com
CONTROLLABILITY OF THE ISHLINSKY SYSTEM
Alexander S. Kuleshov, Vadim V. Rybin
Lomonosov Moscow State University, GSP-1, Leninskie Gory, 1, Moscow, 119991, Russian Federation; kuleshov@mech.math.msu.su, RybinVV@gmail.com
In 1965 A. Yu. Ishlinsky proposed the example of a low dimensional nonholonomic non-Chaplygin mechanical system. The Ishlinsky system consists of three cylinders. One of them with the radius R rolls without sliding on top of two other identical cylinders of a radius a, which each roll whithout sliding on a fixed horizontal plane. In this paper we analyze controllability conditions for the system of cylinders. We prove that the Ishlinsky system is completely controllable by the Chow — Rashevsky theorem. This example of a non-Chaplygin nonholonomic system clearly illustrates the analytic and control coplexity of such systems. Refs 3. Figs 1.
Keywords: nonholonomic system, rolling without sliding, controllability.