Научная статья на тему 'ДИНАМИКА ПЛОСКОГО ДВИЖЕНИЯ КОСМИЧЕСКОГО КРАНА-МАНИПУЛЯТОРА ТИПА РУКИ С УЧЕТОМ ИЗГИБА ЗВЕНЬЕВ'

ДИНАМИКА ПЛОСКОГО ДВИЖЕНИЯ КОСМИЧЕСКОГО КРАНА-МАНИПУЛЯТОРА ТИПА РУКИ С УЧЕТОМ ИЗГИБА ЗВЕНЬЕВ Текст научной статьи по специальности «Физика»

CC BY
45
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КРАН-МАНИПУЛЯТОР ТИПА РУКИ / ДИНАМИКА ПЛОСКОГО ДВИЖЕНИЯ / УПРУГИЕ КОЛЕБАНИЯ / УРАВНЕНИЯ В ОБОБЩЕННЫХ КООРДИНАТАХ

Аннотация научной статьи по физике, автор научной работы — Русских С. В., Шклярчук Ф. Н.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

PLANE MOTION DYNAMICS OF A SPACE-BASED CRANE-MANIPULATOR OF THE ARM TYPE TAKING INTO ACCOUNT THE LINKS’ BENDING

The paper considers the plane motion dynamics of a three-link crane-manipulator consisting of two rods elastic in bending and a massive solid body attached to them at their end (load capture). All links are interconnected by hinge joints with the given (controlled) relative rotation angles. The first link is also connected to the mobile base. The system mathematical model was developed to calculate its non-stationary oscillations and dynamic loading of links at arbitrary in time and large relative angles of the links rotation as the solid bodies. The rods bending deformations (elastic transverse displacements and rotation angles of the axes) were considered to be insignificant (linear). Each rod ending was represented by the Ritz method using two given shape functions with the unknown coefficients, which were taken as the generalized coordinates. In this case, the "cantilever" bending forms of each rod with the given transverse displacement and angle of rotation at their ends were used. A simplified version was considered: each rod bending was represented by only one given form of its bending by a single moment applied at the end. System oscillation equations were obtained based on the principle of probable displacements in the generalized coordinates (four equations or two in a simplified form). Since rotation angles at the ends of the elastic rods were added to the given rotation angles of the non-deformable rods, total angles were left under the signs of sines and cosines without linearization for convenience of calculations and increasing accuracy. As a result, such equations of non-stationary oscillations remained non-linear. They were integrated by numerical methods using standard programs in one of the computer algebra systems. Examples are provided with the comparative results estimates obtained in the two-term and one-term approximations of the function in regard to the bending form of the system bending rod joints.

Текст научной работы на тему «ДИНАМИКА ПЛОСКОГО ДВИЖЕНИЯ КОСМИЧЕСКОГО КРАНА-МАНИПУЛЯТОРА ТИПА РУКИ С УЧЕТОМ ИЗГИБА ЗВЕНЬЕВ»

УДК 534.12 doi: 10.18698/0536-1044-2023-5-112-122

Динамика плоского движения космического крана-манипулятора типа руки с учетом изгиба звеньев*

С.В. Русских1,2, Ф.Н. Шклярчук2,1

1 ФГБОУ ВО «Московский авиационный институт (национальный исследовательский университет)»

2 ФГБУН «Институт прикладной механики Российской академии наук»

Plane motion dynamics of a space-based crane-manipulator of the arm type taking into account the links' bending

S.V. Russkikh12, F.N. Shklyarchuk21

1 Moscow Aviation Institute (National Research University)

2 Institute of Applied Mechanics of the Russian Academy of Sciences

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

Ключевые слова: кран-манипулятор типа руки, динамика плоского движения, упругие колебания, уравнения в обобщенных координатах

* Работа выполнена при финансовой поддержке РНФ (проект № 22-29-01206).

The paper considers the plane motion dynamics of a three-link crane-manipulator consisting of two rods elastic in bending and a massive solid body attached to them at their end (load capture). All links are interconnected by hinge joints with the given (controlled) relative rotation angles. The first link is also connected to the mobile base. The system mathematical model was developed to calculate its non-stationary oscillations and dynamic loading of links at arbitrary in time and large relative angles of the links rotation as the solid bodies. The rods bending deformations (elastic transverse displacements and rotation angles of the axes) were considered to be insignificant (linear). Each rod ending was represented by the Ritz method using two given shape functions with the unknown coefficients, which were taken as the generalized coordinates. In this case, the "cantilever" bending forms of each rod with the given transverse displacement and angle of rotation at their ends were used. A simplified version was considered: each rod bending was represented by only one given form of its bending by a single moment applied at the end. System oscillation equations were obtained based on the principle of probable displacements in the generalized coordinates (four equations or two in a simplified form). Since rotation angles at the ends of the elastic rods were added to the given rotation angles of the non-deformable rods, total angles were left under the signs of sines and cosines without linearization for convenience of calculations and increasing accuracy. As a result, such equations of non-stationary oscillations remained non-linear. They were integrated by numerical methods using standard programs in one of the computer algebra systems. Examples are provided with the comparative results estimates obtained in the two-term and one-term approximations of the function in regard to the bending form of the system bending rod joints.

Keywords: crane-manipulator of the arm type, plane motion dynamics, elastic oscillations, equations in generalized coordinates

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

Космические манипуляторы, предназначенные для работы в условиях невесомости, способны перемещать при плавном, медленном ходе массивные грузы на большие расстояния. Например, манипулятор ERA на Международной космической станции, имеющий два звена в виде тонкостенных стержней общей массой 630 кг, может выдвигаться до 10 м и переносить 8-тонный груз со скоростью до 10 см/с при точности позиционирования до 5 мм [1-3].

Для сборки крупногабаритных космических конструкций, силовым каркасом которых являются регулярные фермы, образованные большим количеством однотипных элементов, могут потребоваться космические роботы-манипуляторы, чтобы выполнять однотипные быстрые операции с высокой точностью позиционирования и устранением упругих колебаний после каждой операции [4-6].

Подходы и методы моделирования динамики манипуляторов с упругими звеньями рассмотрены в работах [7-11]. Многие из этих методов отражены и реализованы на примерах расчета в монографии [12]. Там же подробно

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

В статье [13] для моделирования двухзвен-ного манипулятора со звеньями в виде стержней постоянного поперечного сечения использованы разложения по их собственным формам колебаний с условиями сопряжения в подвижном шарнире.

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

Многие задачи динамики упругих составных систем можно решать с помощью общих подходов и методов, изложенных в монографиях [18-20].

Постановка задачи. Рассмотрим плоское движение космического крана-манипулятора типа руки, состоящего из трех шарнирно связанных между собой звеньев (рис. 1). Относительные углы поворота соединяемых в шарнирных узлах 0, 1, 2 звеньев 00(t), 0i(t), 02(t) и перемещения основания x0(t), y0(t) будем считать заданными (управляемыми).

Первое и второе звенья принимаем упругими на изгиб стержнями с постоянными по длине жесткостями EIi и переменными в общем случае погонными массами mi (£) (i = 1, 2 ), поперечные перемещения, обусловленные изгибом стержней, — малыми, а их продольные перемещения — пренебрежимо малыми.

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

Внутри первого и второго полых стержней расположены механизмы и приводы, осуществляющие заданное жесткое кинематическое управление изменением углов между соединяемыми звеньями 0i (t) (i = 0,1,2) с целью передвижения и позиционирования твердого тела 3.

Для вывода уравнений движения рассматриваемой системы использован принцип возможных перемещений (принцип Даламбера — Лагранжа)

5П-5Аин = 0, (1)

где 8П — вариация потенциальной энергии изгиба первого и второго звеньев; 8Аин — вариация работы инерционных сил рассматрива-

(2)

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

Рассмотрим слагаемые, входящие в состав формулы (1).

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

х(1) = х0 (t) + £ cos ф1 - U(1) sin ф1;

У(1) = у0 (t) + £ sin ф1 + u(1) cos ф1, где £ — продольная координата местной, связанной с левым концом рассматриваемого стержня подвижной системы координат (0<£</1); ф1 = 00(t) — угол поворота первого звена как твердого тела; ü(1)(£, t) — поперечное перемещение первого стержня за счет изгиба.

Малый угол поворота оси стержня за счет его изгиба определяется как du(1)/d £. Прогиб v1(t) = u(1)(/b t) и угол поворота «1(t) на конце стержня далее будем рассматривать как обобщенные координаты.

Вариация работы сил инерции первого звена, взятая с обратным знаком, имеет вид /1

5Аин = J Ш1 [х(1) 5х(1) + У(1) 5у(1) ]d£, (3)

0

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

5у(1) =[(-х0 sinф1 + у0 cosф1) + А1 ]5t>(1);

A1 = ^ф1 (1) -и(1) ср2.

(4)

(5)

Второе звено. Координаты произвольной точки второго стержня х(2)(£, t) и у(2)(£, t) в неподвижной системе координат ху (см. рис. 1) определяются как

х(2) = х1 + £ cos ф2 - U(2) sin ф2;

у(2) = у1 + £ sin ф2 + u(2) cos ф2,

где ф2 — угол поворота второго звена; u(2)(£, t) — поперечное перемещение второго стержня за счет изгиба, 0 <£< /2; х1^) = = х(1>(/ь t); у^) = у(1)(/1, t).

Согласно рис. 1, угол поворота второго звена как твердого тела

ф2 =ф1 + «1 +0x(t),

причем

Рис. 1. Модель трехзвенного космического крана-манипулятора типа руки

5ф2 =5^1.

Прогиб -2(t) = o(2)(¡2, t) и угол поворота V2(t) на правом конце второго стержня за счет его изгиба далее будем рассматривать в качестве обобщенных координат.

Вариация работы сил инерции второго звена, взятая с обратным знаком, имеет вид ¡2

-ЗАин = J Ш2[х (2)Sx(2) + y (2)Sy (2>]d S, (6)

0

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

x (2)Sx(2) + y (2)Sy(2) = = S— {(-X0 sin ф! + y0 cos ф! ) + Ai(¡i) +

+ A2 cos^2 -Ф1) - B2 sin (ф2 -ф! )} +

+ (So(2) + )x{(-X0 sinф2 + y0 cos ф2) + + Ai (¡i ) COs (ф2 -ф! ) + Bi (¡i ) sin (ф2 -ф! ) + A2 } + + (u(2^){(-x¿0 cosф2 - y0 sinф2 )-- Ai (¡i) sin(ф2 -ф!)+Bi (¡i) cos(ф2 -ф!) + B2}. (7) Здесь

A2 = Sf>2 + 0(2) - -0(2)<ф2;

B2 = £ф 2 + 20 (2)ф 2 +^(2)(p2;

Bi = £ф 2 +2-0 (1)ф i + о(1)фь

Третье звено. Координаты центра масс твердого тела в неподвижной системе координат xy (см. рис. i) определяются следующим образом:

Хз = Х2 + r?cos фз; Уз = У2 + r3sin фз, (8)

где ф3 — угол поворота твердого тела. С учетом формулы (5)

X2(t) = x(2)(¡2, t); У2(t) = У(2)(¡2, t);

фз =ф2 + V +02(t);

Sф3 = S#i + SV2. Вариация работы сил инерции с обратным знаком для третьего звена (как для абсолютно твердого тела) определяется выражением

-8А£> = (X3SX3 + y3Sy3) М3 + J 0ф з (( + S#2),

которое с учетом соотношения (8) можно записать как

-8А£> = [(X2SX2 + y2Sy2) + С3Г3 ] M3 + + [M3r3 (-X2 sinф3 + y2 cosф3) + J3ф3 ]x

x( +S#2), (9)

где

C3 = S-i [фз cos ( фз - ф! ) - ф3 sin ( фз - ф! ) ] + + (S-2 +¡2S#i) [фз cos (фз -ф2) - ф^ sin (фз -ф2 )] +

+ -S^i [фз sin (фз -ф2 ) + ф22 cos (фз -ф2)]; (i0)

J 3 = J 30 + М3Г32.

Дифференциальные уравнения движения.

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

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

-(k) =-kVk + Vk¥k, k = i, 2;

S2 S3 S3 S2

Vk(S) = 3^-2s-; ¥k(S) = ; (ii)

¡2 ¡3

¡2 ¡k

Vk(¡k) = i; ¥k(¡k) = 0; v'k(¡k) = 0; y'k(¡k) = i.

С учетом выражения (ii) потенциальная энергия изгиба k-го стержня

nk =

i \l "\2 2EIkJ(k) ) dS--

if11 (- ¡k -k Vk + 1¡2 V

(i2)

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

П = П1 + П2.

С учетом соотношений (3), (6), (9), (4), (7) и (10) выражение для вариации работы инерционных сил рассматриваемой системы в обобщенных координатах , , о2, принимает вид

-8А = -8А(1) - 8А(2) - 8А(3) = = !18^1 + !28#1 + 138о2 +148#2, (13)

где Ъ — нелинейные дифференциальные операторы заданных функций х0^), у0(^), 90(£), 01 (^), 92(£) и неизвестных функций о^), ), 02^), ), г = 1,2,3,4; ввиду громоздкости выражения для них не приведены.

С учетом выражений (1), (11)-(13) система дифференциальных уравнений плоского движения и упругих колебаний трехзвенного крана-манипулятора имеет вид

r ^ l2Eh 6EI1 q L1 +---= 0;

r 6EI1 ^

L2--+

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2 ¡2

т ^ 12EI2 L3 + ——-2

г 6EI2 ^ L4--—-2 +

12

¡12

4 EI11 ¡1

_ 6EI2

" ¡22 4EI2

¡2

►1 = 0; ^2 = 0; $2 = 0.

(14)

При одночленной аппроксимации прогиба к-го стержня как консоли в местной связанной системе координат Ел можно использовать выражение

Е2

и(к) =-0к «к; а к (Е) = к = 1,2. (15) 12

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

h = -к -р ¡k

(16)

Здесь обобщенной координатой для к-го стержня будет только щ (^).

С учетом выражений (15) и (16) потенциальная энергия изгиба к-го стержня

Пк =

1 \ I "\2 1 4Eh

2 EI к J ( ) dE = -к. (17)

Тогда вместо выражения (13) запишем -8Аин = ^ Ъ + 2 Ъ ^ 8^1 + ^ ¿з + ^ й ^ 8^2. (18)

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

~ 2 ~ 4 EI1

L1 + — L2 +--01 = 0;

1 ¡1 2 ¡13 1

~ 2 ~ 4 EI2 L3 + — L4 +о2 = 0.

¡2

(19)

¡ 3

2

Численные решения начальных задач с анализом точности и сходимости для системы, состоящей из четырех (14) и двух (19) дифференциальных уравнений, в обобщенных координатах получены с использованием стандартных программ, реализующих методы Адамса, Рун-ге — Кутты, Булирша — Штера, в программном комплексе компьютерной алгебры PTC Mathcad.

Пример расчета. Проведем численное сравнение результатов интегрирования двух нелинейных систем дифференциальных уравнений второго порядка (14) и (19) на двух модельных примерах расчета.

Во всех примерах расчета программное изменение углов 90(^), 81 (^), 92(£) взято идентичным и показано на рис. 2:

60(f) =

61(t) =

о 1 1 f t 00T---cos I п—

L 22 ^ T

00т при t > T; 0 при 0 < t < T;

0 Г11 f t _ T o1T---cos I П-

L 2 2 ^ T

01T при t > 2T; 0 при 0 < t < 2T;

0 Г11 f t_2T

02T---cos I п-

L 2 2 ^ T

02T при t > 3T,

при 0 < t < T;

при Г < t< 2T;

при 2T < t< 3T;

02 (t) =

где

8ог = —; 81т = —; 82т =--; Т = 15 с.

8 3 4

Узел 0 считаем неподвижным, поэтому

%о(г) = о, уок) = о.

Общие для обоих примеров расчета исходные данные: ¡\ = 12 = 5 м; Ш\ = т2 = 5 кг/м; г3 = 2 м; М3 = боо кг; Ю = 1оо кг-м2. Численное интегрирование выполнялось по методу Адамса в интервале времени от ^ = о до tк = 4Т = бо с, шаг интегрирования — 1о-4 с.

0о; 0ь 02> рад

Рис. 2. Изменение во времени t управляемых углов поворота стержней 8о, 9! и 82

В начальном положении стержни и схват находятся в покое и вытянуты вдоль оси x (00 = 01 = 02 = 0). В первом примере изгибная жесткость обоих упругих звеньев составляла 105 Па-м4, во втором примере — EI1 = EI2 = = 2 -106 Па-м4.

Расчетные зависимости прогибов о1, "2 и углов поворота , Ф2 стержней от времени £ при изгибной жесткости Е11 = Е12 = 105 Па-м4 приведены на рис. 3, а при Е11 = Е12 =

Зь рад

= 2-106 Па-м4 — на рис. 4. Здесь сплошными линиями показаны решения для одночленной аппроксимации прогиба, а штриховыми — для двучленной.

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

рад

-0,2 -

-0,4

0,125 -

0,05 -

-0,025 -

-0,100

Рис. 3. Расчетные зависимости прогибов о1 (я), о2 (б) и углов поворота (в), (г) стержней от времени I при изгибной жесткости Е11 = Е12 = 105 Па-м4

и1;м

0,02

-0,02

-0,04

и2, м

0,01

45 и с

-0,01

-0,02

45 и с

Рис. 4 (начало). Расчетные зависимости прогибов о1 (я), о2 (б) и углов поворота (в), (г) стержней от времени I при изгибной жесткости Е11 = Е12 = 2 • 106 Па-м4

Эьрад

0,01

рад

-0,01

-0,02

0,005

-0,005

-0,010

45 и с

Рис. 4 (окончание). Расчетные зависимости прогибов о1 (я), о2 (б) и углов поворота (в), Ф2 (г) стержней от времени £ при изгибной жесткости Е11 = Е12 = 2 • 106 Па-м4

у, ж

у, м

\

J_I_I I_I_1_

-12 -8 -4 0 4 8 х,м -12 -8 -4 -4

-4

\ )

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

/

I-

0 4 8 х,м

-12 —!

х, м -12 —!

-12 I-

Рис. 5. Формы крана-манипулятора при изгибной жесткости Е11 = Е12 = 105 Па-м4

в разные моменты времени:

я — £ = Т = 15 с; б — £ = 2Т = 30 с; в — £ = 3Т = 45 с; г — £ = 4Т = 60 с

Формы крана-манипулятора в различные Интегрирование по времени системы нели-

моменты времени для первого и второго рас- нейных дифференциальных уравнений (19)

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

рис. 5 и 6. стержней идет в среднем на 40...70 % быстрее,

Рис. 6. Формы крана-манипулятора при изгибной жесткости Ы1 = Е12 = 2 -Ю6 Па-м4

в разные моменты времени: а — t = Т = 15 с; б — t = 2Т = 3о с; в — t = 3Т = 45 с; г — t = 4Т = бо с

чем системы (14) для двучленной аппроксимации.

Выводы

1. Разработана математическая модель плоского движения трехзвенного крана-манипулятора типа руки с двумя упругими на изгиб стержневыми звеньями при кинематическом воздействии трех управляемых относительных углов поворота в соединительных шарнирных узлах. Изгиб каждого стержня описан заданными одно- или двучленными функциями

Литература

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

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

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

[1] Афанасьев И. Начало новой «Эры». Русский космос, 2021, № 10, с. 54-57.

[2] Яскевич А.В., Мирхайдаров В.М., Андреев О.С. Моделирование в реальном времени

сборочных операций, выполняемых на МКС манипулятором Era. Экстремальная робототехника, 2021, т. 1, № 1, с. 19-29.

[3] Голубев Ю.Ф., Яскевич А.В., Мирхайдаров В.М. Оценка динамики причаливания на

компьютерном стенде моделирования в реальном времени. Препринты ИПМ им. М.В. Келдыша, 2020, № 84, doi: https://doi.org/10.20948/prepr-2020-84

[4] Сомов Е.И., Бутырин С.А. Наведение и управление свободнолетающим роботом при

завершении сближения с пассивным объектом в Дальнем космосе. Известия Самарского научного центра РАН, 2017, т. 19, № 4, с. 81-90.

[5] Тебуева Ф.Б., Петренко В.И., Антонов В.О. и др. Методика определения взаимополо-

жения суставов руки оператора для управления антропоморфным космическим манипулятором. Экстремальная робототехника, 2018, т. 1, № 1, с. 69-81.

[6] Даляев И.Ю., Кузнецова Е.М., Шардыко И.В. Перспектива создания роботизированных

сервисных спутников для технического обслуживания и продления сроков активного существования космических аппаратов. Робототехника и техническая кибернетика, 2015, № 3, с. 27-31.

[7] Акуленко Л.Д., Михайлов С.А., Черноусько Ф.Л. Моделирование динамики манипуля-

тора с упругими звеньями. Известия АН СССР. МТТ, 1981, № 3, с. 118-124.

[8] Михайлов С.А., Черноусько Ф.Л. Исследование динамики манипулятора с упругими

звеньями. Известия АН СССР. МТТ, 1984, № 2, с. 51-58.

[9] Черноусько Ф.Л. Динамика управляемых движений упругого манипулятора. Известия

АН СССР. Техн. кибернетика, 1981, № 5, с. 142-152.

[10] Рахманов Е.В., Стрелков А.Н., Шведов В.Н. Разработка математической модели упругого манипулятора на подвижном основании. Известия АН СССР. Техн. кибернетика, 1981, № 4, с. 109-114.

[11] Гукасян А. А. Исследование управляемых движений упругого манипулятора с тремя степенями подвижности. Известия АН АрмССР. Мехакника, 1983, т. 36, № 3, с. 12-20.

[12] Черноусько Ф.Л., Болотник Н.Н., Градецкий В.Г. Манипуляционные роботы. Москва, Наука, 1989. 363 с.

[13] Morris A.S., Madani A. Static and dynamic modelling of a two-flexible-link robot manipulator. Robotica, 1996, vol. 14, no. 3, pp. 289-300, doi: https://doi.org/10.1017/ S0263574700019603

[14] Гуляев В.И., Завражина Т.В. Динамическое управление плоскими движениями упругого двузвенного космического робота-манипулятора. Проблемы управления и информатики, 1998, № 1, с. 140-154.

[15] Гуляев В.И., Завражина Т.В. Динамика управляемых движений упругого робота-манипулятора. Известия АН. МТТ, 1998, № 5, с. 19-28.

[16] Гуляев В.И., Завражина Т.В. Динамика робота-манипулятора с упругоподатливыми звеньями и приводными механизмами. Известия АН. МТТ, 2003, № 6, с. 18-30.

[17] Завражина Т.В. Влияние упругой податливости звеньев на динамику и точность позиционирования робота-манипулятора с вращательными и поступательными сочленениями. Известия АН. МТТ, 2008, № 6, с. 17-32.

[18] Лурье А.И. Аналитическая механика. Москва, Физматлит, 1961. 824 с.

[19] Докучаев Л.В. Нелинейная динамика летательных аппаратов с деформируемыми элементами. Москва, Машиностроение, 1987. 232 с.

[20] Гришанина Т.В., Шклярчук Ф.Н. Динамика упругих управляемых конструкций. Москва, Изд-во МАИ, 2007. 328 с.

References

[1] Afanasyev I. Beginning of new "Era". Russkiy kosmos, 2021, no. 10, pp. 54-57. (In Russ.).

[2] Yaskevich A.V., Mirkhaydarov V.M., Andreev O.S. Real time simulation of assembly opera-

tions executing on ISS by Era robotic arm. Ekstremalnaya robototekhnika [Extreme Robotics], 2021, vol. 1, no. 1, pp. 19-29. (In Russ.).

[3] Golubev Yu.F., Yaskevich A.V., Mirkhaydarov V.M. Berthing dynamic estimation by using

real time computer simulator. Preprinty IPM im. M.V. Keldysha [KIAM Preprint], 2020, no. 84, doi: https://doi.org/10.20948/prepr-2020-84 (in Russ.).

[4] Somov E.I., Butyrin S.A. Guidance and control of free-flying robot during completion of the

rendezvous with a passive object in deep space. Izvestiya Samarskogo nauchnogo tsentra RAN [Izvestia RAS SamSC], 2017, vol. 19, no. 4, pp. 81-90. (In Russ.).

[5] Tebueva F.B., Petrenko V.I., Antonov V.O. et al. A method of determining the mutual posi-

tion of operator's arm joints for anthropomorphic space manipulator control. Ekstremalnaya robototekhnika [Extreme Robotics], 2018, vol. 1, no. 1, pp. 69-81. (In Russ.).

[6] Dalyaev I.Yu., Kuznetsova E.M., Shardyko I.V. Development prospects of robotic service sat-

ellites for maintenance purposes and active lifetime extension of spacecraft. Robototekhnika i tekhnicheskaya kibernetika [Robotics and Technical Cybernetics], 2015, no. 3, pp. 27-31. (In Russ.).

[7] Akulenko L.D., Mikhaylov S.A., Chernousko F.L. Modeling of manipulator dynamics with

elastic linkages. Izvestiya AN SSSR. MTT, 1981, no. 3, pp. 118-124. (In Russ.).

[8] Mikhaylov S.A., Chernousko F.L. Study on manipulator dynamics with elastic linkages.

Izvestiya AN SSSR. MTT, 1984, no. 2, pp. 51-58. (In Russ.).

[9] Chernousko F.L. Dinamika upravlyaemykh dvizheniy uprugogo manipulyatora. Izvestiya AN

SSSR. Tekhn. kibernetika, 1981, no. 5, pp. 142-152. (In Russ.).

[10] Rakhmanov E.V., Strelkov A.N., Shvedov V.N. Dynamics of controlled motions of an elastic manipulator. Izvestiya AN SSSR. Tekhn. kibernetika, 1981, no. 4, pp. 109-114. (In Russ.).

[11] Gukasyan A.A. IStudy on controlled motion of an elastic manipulator with three degrees of freedom. Izvestiya AN ArmSSR. Mekhaknika, 1983, vol. 36, no. 3, pp. 12-20. (In Russ.).

[12] Chernousko F.L., Bolotnik N.N., Gradetskiy V.G. Manipulyatsionnye roboty [Manipulation robors]. Moscow, Nauka Publ., 1989. 363 p. (In Russ.).

[13] Morris A.S., Madani A. Static and dynamic modelling of a two-flexible-link robot manipulator. Robotica, 1996, vol. 14, no. 3, pp. 289-300, doi: https://doi.org/10.1017/ S0263574700019603

[14] Gulyaev V.I., Zavrazhina T.V. Dynamic control of planar motions of an elastic two-link space robot-manipulator. Probtemy upravteniya i informatiki, 1998, no. 1, pp. 140-154. (In Russ.).

[15] Gulyaev V.I., Zavrazhina T.V. Controlled motion dynamics of an elastic manipulator robot. Izvestiya AN. MTT, 1998, no. 5, pp. 19-28. (In Russ.).

[16] Gulyaev V.I., Zavrazhina T.V. Dynamics of a robot manipulator with elastic-fit links and drive mechanisms. Izvestiya AN. MTT, 2003, no. 6, pp. 18-30. (In Russ.).

[17] Zavrazhina T.V. Influence of elastic compliance of links on the dynamics and accuracy of a manipulating robot with rotational and translational joints. Izvestiya AN. MTT, 2008, no. 6, pp. 17-32. (In Russ.). (Eng. version: Mech. SoHds, 2008, vol. 43, no. 6, pp. 850-862, doi: https://doi.org/10.3103/S0025654408060034)

[18] Lurye A.I. AnaUticheskaya mekhanika [Analytic mechanics]. Moscow, Fizmatlit Publ., 1961. 824 p. (In Russ.).

[19] Dokuchaev L.V. NeUneynaya dinamika ¡etatehykh apparatov s deformiruemymi ekmentami

[Nonlinear dynamics of aircraft with deformable elements]. Moscow, Mashinostroenie Publ., 1987. 232 p. (In Russ.).

[20] Grishanina T.V., Shklyarchuk F.N. Dinamika uprugikh upravfyaemykh konstruktsiy [Dynamics of elastic controlled constructions]. Moscow, Izd-vo MAI Publ., 2007. 328 p. (In Russ.).

Статья поступила в редакцию 14.11.2022

Информация об авторах

РУССКИХ Сергей Владимирович — доктор физико-математических наук, доцент, профессор кафедры 602 «Проектирование и прочность авиационно-ракетных и космических изделий». ФГБОУ ВО «Московский авиационный институт (национальный исследовательский университет)»; старший научный сотрудник. ФГБУН «Институт прикладной механики Российской академии наук» (125040, Москва, Российская Федерация, Ленинградский проспект, д. 7, e-mail: [email protected]).

Information about the authors

RUSSIKIKH Sergey Vladimirovich — Doctor of Science (Physics and Math), Associate Professor, Professor of Department 602 — Design and Durability of Aircraft, Rocket and Space Products. Moscow Aviation Institute (National Research University); Senior Researcher. Institute of Applied Mechanics of the Russian Academy of Sciences (125040, Moscow, Russian Federation, Leningradskiy Ave., Bldg. 7, e-mail: [email protected]).

ШКЛЯРЧУК Федор Николаевич — доктор технических наук, профессор, главный научный сотрудник. ФГБУН «Институт прикладной механики Российской академии наук»; профессор кафедры 602 «Проектирование и прочность авиационно-ракетных и космических изделий». ФГБОУ ВО «Московский авиационный институт (национальный исследовательский университет)» (125993, Москва, Российская Федерация, Волоколамское шоссе, д. 4, e-mail: [email protected]).

SHKLYARCHUK Fyodor Nikolaevich — Doctor of Science (Eng.), Professor, Chief Researcher. Institute of Applied Mechanics of the Russian Academy of Sciences; Professor of Department 602 — Design and Durability of Aviation-Rocket and Space Products. Moscow Aviation Institute (National Research University) (125993, Moscow, Russian Federation, Volokolamskoye Shosse, Bldg. 4, e-mail: [email protected]).

Просьба ссылаться на эту статью следующим образом:

Русских С.В., Шклярчук Ф.Н. Динамика плоского движения космического крана-манипулятора типа руки с учетом изгиба звеньев. Известия высших учебных заведений. Машиностроение, 2023, № 5, с. 112-122, doi: 10.18698/0536-1044-2023-5-112-122

Please cite this article in English as: Russkikh S.V., Shklyarchuk F.N. Plane motion dynamics of a space-based crane-manipulator of the arm type taking into account the links' bending. BMSTU Journal of Mechanical Engineering, 2023, no. 5, pp. 112-122, doi: 10.18698/0536-1044-2023-5-112-122

Издательство МГТУ им. Н.Э. Баумана предлагает читателям монографию

«Алгоритмы повышения точности предпосадочной подготовки беспилотных летательных аппаратов авианосного базирования»

Авторы: Л. Чжан, К.А. Неусыпин, М.С. Селезнева, А.В. Пролетарский

Исследованы алгоритмы обработки информации современных навигационных систем и комплексов беспилотных летательных аппаратов. Представлены наиболее популярные численные критерии степени наблюдаемости, управляемости и идентифицируемости, использующиеся для определения качественных характеристик моделей. Разработан оригинальный численный критерий степени наблюдаемости переменных состояния одного класса нелинейных моделей с применением SDC-метода. Приведен способ оценки степени управляемости переменных состояния моделей критерия степени наблюдаемости и принципа дуальности. Разработана методика выбора используемого в НК информационного поля и навигационных систем НК для проведения коррекции ИНС с учетом условий функционирования БЛА. Предложены критерии формирования динамической рабочей структуры НК в процессе полета БЛА в меняющейся помеховой обстановке, базирующиеся на анализе оценок дисперсий приведенных измеряемых шумов.

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

По вопросам приобретения обращайтесь:

105005, Москва, 2-я Бауманская ул., д. 5, стр. 1. Тел.: +7 499 263-60-45, факс: +7 499 261-45-97; [email protected]; https://bmstu.press

i Надоели баннеры? Вы всегда можете отключить рекламу.