Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

МАТЕМАТИЧЕСКИЕ МОДЕЛИ СИСТЕМ СВЯЗАННЫХ ТВЁРДЫХ ТЕЛ В ИМПУЛЬСАХ ПУАССОНА

Иванов В.Н. 1 Полосков И.Е. 1 Шимановский В.А. 1
1 ФГБОУ ВО «Пермский государственный национальный исследовательский университет»
Работа посвящена разработке методов компьютерного моделирования динамики сложных механических систем. В статье представлена новая матричная форма уравнений движения систем абсолютно твердых тел со структурой дерева. Рассмотрен случай голономных связей. В качестве независимых параметров, однозначно определяющих положение и распределение скоростей тел системы в пространстве, выбраны обобщенные координаты и переменные, имеющие размерность импульсов. Особенность системы уравнений состоит в том, что она разрешена относительно производных от обобщённых импульсов и не содержит реакций связей. Приведен вывод предлагаемой новой формы уравнений движения системы связанных твердых тел из уравнений динамики свободного твердого тела в квазикоординатах Эйлера – Лагранжа с использованием матрично-геометрического подхода. Предложен метод разрешения уравнений движения относительно старших производных, ориентированный на использование ЭВМ. Получены рекуррентные формулы для определения всех кинематических и динамических переменных, входящих в уравнения. Временна?я сложность решения системы уравнений с помощью данного алгоритма растёт по линейному закону в зависимости от числа тел в механической системе, что говорит о его эффективности. Рассматриваемый алгоритм подобен методу, предложенному А.Ф. Верещагиным для разрешения общих уравнений динамики цепочки твёрдых тел относительно ускорений. На примере механической системы с четырьмя степенями свободы продемонстрированы все этапы подготовки первичной информации и составления уравнений движения в предлагаемой форме.
система абсолютно твердых тел
уравнения движения
динамика
математическое моделирование
обобщенные координаты
импульсы Пуассона
матрично-геометрический метод
1. Величенко В.В. Матрично-геометрические методы в механике с приложениями к задачам робототехники. – М.: Наука, 1988. – 280 с.
2. Верещагин A.Ф. Метод моделирования на ЦВМ динамики сложных механизмов роботов-манипуляторов // Изв. АН СССР. Техническая кибернетика. – 1974. – № 6. – С. 89–94.
3. Верещагин A.Ф. Принцип наименьшего принуждения Гаусса для моделирования на ЭВМ динамики роботов-манипуляторов // Докл. АН СССР. – 1975. – Т. 220, № 1. – С. 51–53.
4. Иванов В.Н. Математическое моделирование динамики механических систем со структурой дерева.– Пермь: Перм. ун-т. – 1983. – 11 с. – Деп. в ВИНИТИ: 28.10.83, № 5876.
5. Иванов В.Н., Домбровский И.В., Набоков Ф.В. и др. Классификация моделей систем твердых тел, используемых в численных расчетах динамического поведения машиностроительных конструкций // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. – 2012. – № 2. – С. 139–155.
6. Иванов В.Н., Суслонов В.М. Уравнения движения механических систем со структурой дерева // Проблемы механики управляемого движения. Нелинейные динамические системы. – Пермь, 1984. – С. 154–159.
7. Лилов Л.К. Моделирование систем связанных тел. – М.: Наука, 1993. – 272 с.
8. Лурье А.И. Аналитическая механика. – М.: Гос. изд-во физ-мат. лит, 1961. – 824 с.
9. Шимановский В.А., Иванов В.Н. Формирование уравнений движения механических систем в обобщенных координатах // Проблемы механики и управления. Нелинейные динамические системы. – Пермь, 2005. – Вып. 37. – С. 188–201.
10. Шимановский В.А., Иванов В.Н. Методы составления уравнений движения систем связанных твердых тел в декартовых координатах // Проблемы механики и управления: Нелинейные динамические системы. – 2007. – Вып. 39. – С. 248–262.
11. Wittenburg J. Dynamics of multibody systems. – Berlin: Springer-Verlag, 2008. – 223 p.
12. Shabana A.A. Computational dynamics. – New York: Wiley, 2009. – 542 p.

Компьютерное моделирование динамики механических систем широко используется в современной инженерной практике. Оно позволяет существенно уменьшить объём натурных испытаний и в конечном счёте сократить время и стоимость новых разработок. Требование точности компьютерного моделирования ведет к возрастанию сложности математических моделей технических систем. Но параллельно с ростом размерности математической модели увеличивается и трудоёмкость моделирования. Это отражается в серьезном росте времени, которое затрачивается на стадии оптимизации параметров проектируемой конструкции и исследовании ее поведения при различных условиях эксплуатации. Поэтому разработка методов, позволяющих ускорить процесс математического моделирования сложных технических систем, является актуальной задачей.

В настоящей работе выводится новая форма уравнений движения систем связанных твёрдых тел с использованием квазискоростей, обобщённых координат и обобщённых импульсов (импульсов Пуассона).

Описание механической системы. Матрица кинематической структуры

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

Пусть N – число тел в системе (не считая тела «0», движение которого во времени относительно инерциальной системы координат (СК) Oxyz задано). Тогда количество шарниров всегда будет равно N, если учитывать шарнир между системой и телом «0».

Пронумеруем тела и шарниры системы, причем телам присвоим номера так, чтобы для любого тела в графе системы номер предшествующего ему тела был меньше, а шарниру, связывающему i-е тело с предшествующим, присвоим номер i. В этом случае для полного описания структуры взаимосвязей в такой системе достаточно одного целочисленного массива k = {k1, …, kN}, на i-м месте которого расположен индекс тела, предшествующего i-му. C каждым телом системы свяжем следующие множества: Pi – упорядоченное множество индексов шарниров, составляющих путь между нулевым и i-м телами; Ui – множество индексов шарниров, для которых i-е тело является предшествующим.

Положение и ориентацию i-го тела системы относительно инерциальной СК будем определять с помощью радиуса-вектора Ivanov01.wmf фиксированной на нём точки Oi и системы ортонормированных векторов Ivanov02.wmf Ivanov03.wmf Ivanov04.wmf которые вместе с точкой Oi определяют СК Oixiyizi, неизменно связанную с i-м телом.

Введем следующие обозначения: Ivanov05.wmf – матрица-столбец координат точки Oi в ki-й СК; ri = OOi – матрица-столбец координат точки Oi в инерциальной СК; Ivanov06.wmf – матрица направляющих косинусов между базисными векторами Ivanov07.wmf Ivanov08.wmf Ivanov09.wmf и Ivanov10.wmf Ivanov11.wmf Ivanov12.wmf (матрица преобразования координат из j-й СК в i-ю).

Предположим, что физическая связь в i-м шарнире моделируется системой li ≤ 6 идеальных голономных линейно-независимых связей. Тогда множество относительных положений i-го тела в ki-й СК образует ni-мерное конфигурационное многообразие [1], где ni = 6 – li. Уравнения этого многообразия можно записать в параметрической форме, введя матрицу-столбец Ivanov13.wmf криволинейных (обобщенных) координат его текущей точки. Тогда матрицы ρi и Ivanov14.wmf являются функциями обобщенных координат:

ρi = ρi(qi, t); Gi = Gi(qi, t).

Задание матрицы-столбца q = (q1, …, qN)T как функции времени позволяет однозначно определить взаимное положение всех тел и их движение относительно абсолютной системы координат в любой момент t.

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

Ivanov15.wmf Ivanov16.wmf

Принимая движение i-го тела за относительное, а предшествующего ему ki-го тела за переносное, в соответствии с правилом сложения скоростей можно записать рекуррентные формулы для вычисления проекций линейной νi и угловой ωi скоростей тел механической системы на оси i-й СК [8]:

Ivanov17.wmf (1)

где

Ivanov18.wmf Ivanov19.wmf

Ivanov20.wmf Ivanov21.wmf

Ivanov22.wmf Ivanov23.wmf

Ivanov24.wmf Ivanov25.wmf

Здесь и далее символ «~» используется для обозначения кососимметричной матрицы (это соответствует матричной записи векторного произведения [8]).

Введём блочную 6N×6N-матрицу S с квадратными подматрицами порядка 6 по следующей формуле:

Ivanov26.wmf Ivanov27.wmf (2)

Заметим, что для любой кинематической структуры эта матрица содержит в каждой строке только два ненулевых блока E6 и –Ci. Поскольку матрица S содержит информацию как о топологической структуре системы, так и об относительном положении тел в системе, ее называют матрицей кинематической структуры. С использованием матрицы S уравнения кинематики системы тел (1) можно записать следующим образом:

Ivanov28.wmf (3)

где v = (v1, ω1, …, vN, ωN)T, A = diag(A1, …, AN).

Рекуррентные формулы (3) можно записать в виде явных выражений

Ivanov29.wmf (4)

где обратная к S матрица T =S–1 является блочной 6N×6N-матрицей, подматрицы которой могут быть вычислены по рекуррентным формулам:

Ivanov30.wmf Ivanov31.wmf (5)

Для дальнейшего изложения нам потребуется соотношение для производной матрицы T. Для вывода этого соотношения получим формулу для производной матрицы Ci:

Ivanov32.wmf

где введено обозначение

Ivanov33.wmf

Из этого равенства и формулы (2) следует, что

Ivanov34.wmf (6)

где Ω = diag(Ω1, …, ΩN).

Продифференцируем равенство ST = E и из полученного выражения выразим производную Ivanov35.wmf. Тогда с учётом формулы (6) получим

Ivanov36.wmf (7)

Расширенная система уравнений движения в импульсах Пуассона

Для составления уравнений движения систем связанных твёрдых тел со структурой дерева будем использовать принцип освобождения от связей. В этом случае мысленно уничтожаем связи между телами, а движение системы, осуществляемое при наличии связей, обеспечим введением дополнительных сил – реакциями связей. Введение этих реакций позволяет записать дифференциальные уравнения движения любого тела в форме уравнений Эйлера – Лагранжа для свободного твёрдого тела [8]:

Ivanov37.wmf (8)

где

Ivanov38.wmf Ivanov39.wmf

Ivanov40.wmf

mi – масса i-го тела; Ji – тензор инерции i-го тела; Ivanov41.wmf – матрица-столбец проекций радиуса-вектора центра масс i-го тела на оси связанной с ним СК; Ivanov42.wmf Ivanov43.wmf – проекции главного вектора и главного момента активных сил, действующих на i-е тело в i-й СК; Ivanov44.wmf Ivanov45.wmf – проекции главного вектора и главного момента реакций связей в i-м шарнире в i-й СК.

С использованием матрицы кинематической структуры S уравнение (8) запишем в краткой форме:

Ivanov46.wmf (9)

где

M = diag(M1, …, MN); F = (F1, …, FN)T;

R = (R1, …, RN)T.

Декартовым импульсом системы назовём 6N-мерную матрицу-столбец

pa = Mv. (10)

Тогда n-мерная матрица-столбец Ivanov47.wmf

p = ATTTpa = BTpa

является обобщённым импульсом, или импульсом Пуассона [1].

Построим выражение для производной обобщённого импульса Ivanov48.wmf. Для этого умножим слева на BT левую и правую части уравнения (9). Получим, что

Ivanov49.wmf

или

Ivanov50.wmf (11)

Здесь было учтено равенство

BTSTR = ATTTSTR = ATR = 0,

отражающее принцип идеальности связей. В уравнении (11) Q = BTF – матрица-столбец обобщенных сил.

Преобразуем выражение Ivanov51.wmf с учётом формулы (7):

Ivanov52.wmf

где введено обозначение μ = TT pa.

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

Ivanov53.wmf (12)

Особенность уравнений (12) состоит в том, что они разрешены относительно производных обобщённых импульсов Ivanov54.wmf. Первые три из этих уравнений образуют линейную систему с симметричной, блочной трёхдиагональной разреженной матрицей коэффициентов относительно скоростей v, Ivanov55.wmf и переменных μ, которые являются множителями Лагранжа [4, 6]. Уравнения (12) для цепочки твердых тел впервые были выведены В.Н. Ивановым из принципа Гамильтона – Остроградского в работе [4]. В работе [6] эти уравнения были получены в координатной форме для систем твердых тел со структурой дерева.

Структура первых трех кинематических уравнений системы (12) аналогична структуре расширенной системы уравнений динамики системы связанных твердых тел [5, 9, 10]:

Ivanov56.wmf (13)

где w – матрица-столбец абсолютных линейных и угловых декартовых ускорений всех тел системы в связанных с ними СК. Правые части уравнений F* и w*, помимо внешних сил, содержат гироскопические, центробежные и явно зависящие от времени ускорения.

Отличие уравнений (12) и (13) состоит в том, что первые три уравнения системы (12) служат для получения значений обобщенных скоростей Ivanov57.wmf, а не ускорений Ivanov58.wmf.

Совпадение структуры уравнений (12) и (13) означает, что все разработанные ранее алгоритмы решения системы уравнений (13), такие как метод «прогонки» (отдельных тел) А.Ф. Верещагина [2, 3], методы проекций уравнений в подпространства обобщенных координат или подпространства линейно-независимых компонент реакций связей (схемы явного формирования уравнений движения в форме Лагранжа 2 или 1 рода) [7, 8, 11, 12], применимы и для системы уравнений (12).

Решение расширенной системы уравнений движения в импульсах Пуассона методом «прогонки»

Учитывая формулы (1) и (2), систему уравнений (12) запишем в следующем виде:

Ivanov59.wmf (14)

Ivanov60.wmf (15)

Ivanov61.wmf (16)

Ivanov62.wmf (17)

где

Ivanov63.wmf Ivanov64.wmf (18)

Построим на основе этой системы рекуррентные формулы для вычисления обобщённых скоростей Ivanov65.wmf и множителей Лагранжа μi. Для этого сначала докажем, что для каждого тела системы множитель μi может быть представлен в виде

Ivanov66.wmf (19)

Нетрудно видеть, что для концевых тел системы такое представление наблюдается:

Ivanov67.wmf Ivanov68.wmf

Предположим, что формула (19) справедлива для всех номеров j ∈ Si. Покажем, что она справедлива и для номера i. Для этого с помощью формулы (15) исключим из (19) матрицу-столбец vi:

Ivanov69.wmf (20)

Подставим полученное равенство в уравнение (16) и выразим из него матрицу-столбец Ivanov70.wmf:

Ivanov71.wmf (21)

С помощью формул (20), (21) из уравнения (14) исключим множители μj, j ∈ Si, несомых тел:

Ivanov72.wmf

Собирая в полученном равенстве коэффициент перед vi, преобразуем его к виду (19), где

Ivanov73.wmf (22)

Ivanov74.wmf (23)

Таким образом, процесс вычисления обобщённых скоростей и производных обобщённых импульсов состоит из двух этапов. На первом из них по формулам (22), (23) и (18), начиная с концевых тел и заканчивая корневым, вычисляются матрицы Ivanov75.wmf Ivanov76.wmf Ivanov77.wmf На втором этапе, начиная с корневого тела, с помощью формул (21), (15), (19) и (17) определяются обобщённые скорости Ivanov78.wmf и производные обобщённых импульсов Ivanov79.wmf.

Данный метод является одним из самых эффективных методов численного моделирования систем с длинными кинематическими цепочками. Трудоёмкость решения системы уравнений (12) с помощью данного алгоритма растёт по линейному закону в зависимости от числа тел в механической систем [10]. При реализации этого алгоритма требуется обращение только матриц Ivanov80.wmf, порядок которых равен числу степеней свободы в i-м шарнире, причем эти матрицы симметричны и положительно определены, а их порядок всегда мал (не превышает шести). Именно этим и обусловлена эффективность представляемого метода.

Описанный алгоритм аналогичен методу отдельных тел А.Ф. Верещагина, предложенному в работах [2, 3] для решения общих уравнений динамики цепочки твёрдых тел относительно ускорений.

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

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

pic_14.tif

Введём неподвижную Oxyz и связанные с телами Oξ1η1ζ1, Oξ2η2ζ2 системы координат, как показано на рисунке. Система состоит из двух тел с массами m1 и m2 и осевыми моментами инерции J1x, J1y, J1z и J2x, J2y, J2z. Предположим, что:

а) в точке O цилиндрический шарнир обеспечивает вращение OA вокруг оси Oz;

б) точка A может перемещаться вдоль оси Oη1;

в) шарнир в этой точке обеспечивает вращение второго тела сначала вокруг оси Oη1, а затем вокруг оси Oξ2;

г) на тела системы действует сила тяжести;

д) точки O и A связаны пружиной.

Рассматриваемая система имеет четыре степени свободы. В качестве обобщённых координат выберем:

а) угол γ1 – угол поворота первого тела вокруг оси Oz;

б) изменяющуюся длину OA = y2;

в) углы β2 и α2 поворота второго тела вокруг осей Oη1 и Oξ2 соответственно, т.е. q1 = γ1; q2 = (y2, β2, α2)T, q = (y1, y2, β2, α2)T.

Подготовим элементы уравнений движения (12). Очевидно, что координаты точек O и A в системах координат Oxyz и Oξ1η1ζ1 соответственно равны

ρ1 = (0, 0, 0)T; ρ1 = (0, y2, 0)T.

Далее, матрицы направляющих косинусов между системами координат имеют вид

Ivanov81.wmf

Ivanov82.wmf

где

cγ1 = cos γ1; sγ1 = sin γ1; cα2 = cos α2;

sα2 = sin α2; cβ2 = cos β2; sβ2 = sin β2.

Масс-инерционные характеристики тел системы задаются матрицами

M = diag(M1, M2); rc1 = (0, yc1, 0)T;

rc2 = (0, 0, zc2)T,

где

Ivanov83.wmf

Ivanov84.wmf

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

Ivanov85.wmf A = diag(A1, A2),

где

Ivanov86.wmf Ivanov87.wmf

Ivanov88.wmf

Отсюда производная матрицы A равна

Ivanov89.wmf

где

Ivanov90.wmf

Вследствие того, что тело 0 системы стационарно, то вектор v* = 0. При этом вектор обобщённых сил имеет следующий вид:

Ivanov91.wmf

где ω0 = ω0(t) – заданная угловая скорость вращения первого тела; k – коэффициент усиления; c – коэффициент жёсткости пружины; y20 – длина ненапряжённого состояния пружины; g – ускорение свободного падения.

Теперь для получения уравнений движения в импульсах Пуассона достаточно подставить выписанные матрицы в уравнения (12).

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

Заключение

В работе получена новая расширенная матричная форма уравнений движения систем связанных твердых тел со структурой дерева. Уравнения содержат импульсы Пуассона, обобщенные координаты, квазискорости и множители Лагранжа. Они разрешены относительно производных импульсов Пуассона. Поэтому в отличие от классических уравнений в форме Лагранжа 1 или 2 рода представленные уравнения не требуется разрешать относительно старших производных (обобщенных ускорений) в процессе их численного интегрирования.

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

Для группы кинематических уравнений, входящих в расширенную систему уравнений, разработан матричный алгоритм разрешения относительно обобщенных скоростей, аналогичный известному методу «прогонки». Число арифметических операций в построенном алгоритме, как и в методе А.Ф. Верещагина [2, 3], растет линейно в зависимости от количества тел в механической системе. Поэтому для систем твердых тел с длинными кинематическими цепями разработанный алгоритм является более эффективным [10] в сравнении с процедурами моделирования, основанными на использовании уравнений в форме Лагранжа 2 рода, для которых аналогичная зависимость является кубической.


Библиографическая ссылка

Иванов В.Н., Полосков И.Е., Шимановский В.А. МАТЕМАТИЧЕСКИЕ МОДЕЛИ СИСТЕМ СВЯЗАННЫХ ТВЁРДЫХ ТЕЛ В ИМПУЛЬСАХ ПУАССОНА // Фундаментальные исследования. – 2016. – № 10-3. – С. 493-499;
URL: http://www.fundamental-research.ru/ru/article/view?id=40883 (дата обращения: 18.09.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074