УДК 531.38
БОТ: 10.18698/2308-6033-2018-7-1786
Конечное поперечное передвижение упругого стержня с массой на конце при гашении колебаний в момент остановки
© СВ. Русских1, 2, Ф.Н. Шклярчук1, 2
Московский авиационный институт (национальный исследовательский университет), Москва, 125993, Россия 2Институт прикладной механики Российской академии наук, Москва, 125040, Россия
Рассмотрены нестационарные поперечные колебания упругого консольного стержня с массой на правом конце при управляемом перемещении левого конца вдоль вертикальной направляющей. Изгибные колебания стержня отыскиваются в виде разложения по собственным формам колебаний. Поставлена следующая задача: найти закон управляемого передвижения левого конца стержня, при котором он перемещается на требуемое расстояние за заданное время при одновременном гашении упругих колебаний по нескольким низшим собственным формам в момент остановки. В качестве неизвестной функции рассмотрено ускорение левого конца стержня, которое считают пропорциональным некоторой финитной функции. Эту функцию в свою очередь, записывают в виде ряда по синусам или косинусам с неизвестными коэффициентами на заданном времени управляемого движения. Приведены примеры расчета с гашением упругих колебаний по первым трем собственным формам.
Ключевые слова: упругий консольный стержень, нормальные координаты, гашение упругих колебаний, финитная функция, конечные перемещения
Введение. Задачи управления составными упругими системами, которые могут совершать различные конечные передвижения (перемещения и повороты), являются актуальными для быстроходных ма-нипуляционных роботов [1, 2 и др.], виброударных механизмов, ракет с разделяющимися ступенями, различных космических конструкций и т. д. Общие вопросы теории управления большими механическими системами, задачи и методы управления колебаниями конструкций рассмотрены в [3, 4]. В большинстве работ упругие колебания, которые возникают при конечных перемещениях, поворотах механических систем или их составных частей, считаются малыми и описываются нормальными координатами, представляющими собой относительные движения по собственным формам колебаний. Задачи терминального управления с учетом начальных и конечных условий в такой постановке сводятся к системам интегральных уравнений для нормальных координат.
Для того чтобы решить получающиеся интегральные уравнения для собственных форм, которые следует погасить, управляющие функции отыскивают в виде кусочно-линейных функций или рядов
по заданным функциям с неизвестными коэффициентами. В работе
[2] для управления поперечными изгибными колебаниями балки с грузом на конце в качестве таких функций использованы синусы и косинусы с частотами собственных колебаний балки. В этом случае искомые управляющие воздействия имеют вид трудно реализуемых на практике рядов с высокочастотными гармониками по синусам и косинусам. Воспроизвести такие управляющие воздействия с необходимой точностью практически невозможно: даже незначительное отклонение реального времени управления от расчетного или весьма малые различия между собственными частотами реальной системы и ее математической модели могут привести к тому, что конечные условия для нормальных координат, которые должны быть нулевыми, не будут выполняться.
В работах [5-8] предложен альтернативный подход к определению силовых или кинематических управляющих воздействий для конечных передвижений упругой системы из одного состояния в другое за заданное время при гашении упругих колебаний по нескольким низшим собственным формам в конечный момент времени. Для этого используются «простые» управляющие финитные функции с «настройкой» на них нескольких низших собственных частот колебаний системы.
В работах [8-10] решение линейной задачи терминального управления строят в нормальных координатах с разложением управляющей функции в ряд по синусам или косинусам с неизвестными коэффициентами, для которых алгебраические уравнения получают путем удовлетворения начальных и конечных условий для нормальных координат.
В данной статье на основании такого подхода приведено решение задачи о конечном поперечном перемещении стержня с массой на конце на заданное расстояние за заданное время при гашении в момент остановки упругих колебаний по первым трем низшим собственным формам колебаний. Эта же задача в [2] была решена с использованием разложения в ряды Фурье по высокочастотным гармоникам.
Постановка задачи. Будем рассматривать поперечные колебания в вертикальном направлении упругого стержня, консольно прикрепленного к подвижному основанию (рис. 1). Система имеет следующие параметры: EI и m — изгибная жесткость и масса 1 м стержня; До и — масса подвижного основания и масса, присоединенная на конце стержня х = l. Перемещение левого конца стержня х = 0 и основания ц (t) отыскивают на интервале 0 < t < T, где T — заданное
время управляемого конечного перемещения, которое подчиняется начальным и конечным условиям:
U о (0) = Uo(0) = 0; и o(T) = 0; Uo(T) = ц.. (1)
Рис. 1. Общий вид рассматриваемой системы
Дополнительное поперечное перемещение и(^, /) в результате изгиба стержня представляют в виде разложения по собственным формам колебаний Уп (£), £ = х/1:
О = У УпЮШ п _ 1,2, ...,
(2)
где /п (/) — неизвестные нормальные координаты. Собственные формы колебаний неподвижно закрепленной консольной балки постоянного поперечного сечения с массой | на конце при нормировке
Уп (1) = 1 записывают в виде
Уп (£) = Сп (С08И(Кп£) - С08(Кп£)) - ^ (^ККп£) - Кп^ С _ 8шЬ(Кп ) + 8Ш(Кп ) „ _ С08И(Кп ) + С08(Кп ) С _ Ап ' "п _ Лп '
Ап _ 2(С08И(Кп) 8ш(Кп) - 8тИ(Кп)С08(Кп)).
Параметры кп для п _ 1, 2, ... находят из частотного уравнения
1 + С08И(к) С08(К) - — к(С08И(К) 8Ш(К) - 8шИ(К) С08(К)) _ 0. т1
После этого определяют собственные частоты колебаний юп:
2 Е1 4
Юп п _ 1,2, ••• .
т1
Дифференциальные уравнения движения системы в нормальных координатах записывают в виде
тп(!п +®2п/п) _-топо^Х п _1,2, • ••;
топ _¡ГтУпй^ + ^Уп(1); тп _¡ГтУ^^ + |^2(1).
(3)
о о
Реакция при конечном управляемом перемещении системы с массой Мо _|о +т1 +1 может быть найдена по выражению
п
Р0^) = М033 0 Ш0 п /. (4)
п
При конечном перемещении рассматриваемой системы из состояния покоя при г = 0 и гашении N низших собственных форм колебаний в момент окончания операции г = Т получаем следующие терминальные условия для /п:
/п(0) = /п (0) = 0, п = 1,2,...; /п(Т) = /п (Т) = 0, п = 1,2,..., N.
Вместо перемещения 30 (г) в качестве неизвестной функции будем рассматривать ускорение основания 30(г):
30 = Т^ ф(г), (6)
где ф(г) — безразмерная финитная функция при 0 < г < Т. С учетом (1) и (6) ф(г) должна удовлетворять условиям [9, 10]
Т I Т
|ф(г )йг = 0; зТ =--2 {ф(* )*Ж. (7)
0 Т 0
Из уравнений колебаний системы (3) с учетом (6) получаем:
/ (¡) = -Ш_ Шоп! Г ф(т)8;п[Юп (г-Х)¥х при 0 <г < Т;
Шп тп Т2 0 (8)
1 т 1 Т
/п (г) =----Т2 |ф(ф!п[Шп(г -х)]ёх при г >Т.
Шп тп Т 0
Решение в виде ряда по синусам. Представим функцию ф(г) в виде ряда по синусам [9, 10], которые удовлетворяют начальным условиям (5):
ф(г) = 2 ьк *п(кОг); О = —, к = 1,2, ... . (9)
к 1
Здесь Ък — неизвестные коэффициенты разложения, подлежащие определению. Подставляя (9) в выражения (8), получаем:
/(г) = -тп2Ък 2 ^[^(кОг)-—втКг)] при 0<г<Т;
тп Т2 Т ®п - к2О2 Шп
/(г) = -±010*1^^(г-Т))-8шКг)]2Ък 2 кО, 2 при г>Т.
Шп тп Т2 к Ш2п - к2О2
(10)
Для гашения N собственных форм колебаний в момент I = Т из конечных условий (5), а также с учетом второго выражения (7) получаем следующую систему N +1 линейных уравнений для нахождения коэффициентов Ък, к = 1, 2,..., N +1:
N+1 к О N+1 ъ 2-л-
2 К = 0; 2 к = Т»Т, п = 1,2,.... N. (11)
к=1 -п -к О к=1 к /
Таким образом, во всех суммах по к в выражениях (9) и (10) необходимо положить к = 1, 2,..., N +1.
Приведем формулы для коэффициентов Ък системы (11) при N = 1, 2,3.
Для N = 1
= 8гсиТ(-2 -О2) ; = 4гсиТ(-2 -4О2)
Ъ1 = 2 ' Ъ2 = 2 .
3/-1 3/-1
Для N = 2
= 3гсиТ(-2 -О2)(-2 -О2) =-12жуТ(-2 -4О2)(- -4О2) Ъ1 = , 2 2 ' Ъ2 = 2 2 ;
/-1 -2 5/-1 -2
= 3гсиТ (-2 - 9О2)(-2 - 9О2)
Ъ3 = ^^ .
Для N = 3
5/-1 -2
1 б^ит (-2 - о2)(-2 - о2)(-2 - о2)
Ъ3 =
Ъл =-
С 2 2 2 5/-1 -2-3 5
2 1бл:иТ (- - 4О2)(-2 - 4О2)(-2 - 4О 2);
С 2 2 2 5/-1 -2-3 5
2 48л:иТ (- - 9О2)(-2 - 9О2)(-2 - -9О2) ;
35/-2-2-2 5
2 8л:иТ(-1 - 1бо 2)(-2 - 16о2)(-2 - 16О2)
4 2 2 2
35/-1 -2-3
Решение в виде ряда по косинусам. Функция ф(^) может быть также представлена в виде ряда по косинусам [9, 10], которые удовлетворяют начальным условиям (5):
ф(0 = 2 % со{кО') О = у, к = 1,3,..., (12)
2
где а, — неизвестные коэффициенты. Подставляя (12) в выражения для нормальной координаты (8), будем иметь
/п (г) = 2 * 1
О Т2 ^ к шп - к2О2/4
, ч I ко
с°Б(шпг) - с°Б | г
при 0< г < Т;
о / 1
/п (г) ^^32[с08(шп (г - Т)) + с°8(шпг )]2 ак 2 ,2„2/л при г > Т.
тп Т2 Т шп - к2о74
(13)
Для гашения N собственных форм колебаний в момент окончания операции из конечных условий (5) и второго выражения (7) получаем линейную систему N +1 уравнений для нахождения коэффициентов Ък, к = 1,3,..., 2 N +1:
2 N+1 1 2 N+1 а 2
2 ак 2 ,2О2/4 = 2 -Ок = \^ п = 1,2, . Ж (14)
к=1,3 ... шп - к О /4 к=1,3... к 21
Таким образом, во всех суммах по к в выражениях (12) и (13) по аналогии с предыдущим случаем необходимо положить к = = 1, 3,..., 2N +1.
Приведем формулы для коэффициентов о, системы (14) при N = 1, 2,3.
Для N = 1
= 9 л2иТ(4ш2 -О2) ; = 9л2иТ(4ш2 -9О2)
а1 = 2 ; а3 = 2 .
64/ш2 64/ш2
Для N = 2
= 75л2иТ(4ш2 -О2)(4ш2 - О2)
а1 = Г" 9 2 ;
2048/ш2ш2
а =-
225л2зТ (4ш2 - 9О2)(4ш2 - 9О2) ;
а3 = ,99 ;
4096/ш2 ш2 = 75л2иТ (4ш2 - 25О2)(4ш2 - 25О2)
а5 =-2—2-.
4096/ш2ш2
Для N = 3
а=
1225л2зТ(4ш2 -О2)(4ш2 -О2)(4ш2 - О2)
1 131072/ш2ш2ш2
а =-
2205л2иТ (4-2 - 9О2)(4-2 - 9О2)(4-2 - 9О2);
3 131072/-2-2-2
а5 =
1225л2иТ(4-2 - 25О2)(4-2 - 25О2)(4-2 - 25О2) ;
131072/-2-2-2 '
245л2иТ(4-2 - 49О2)(4-2 - 49О2)(4-2 - 49О2)
а7 131072/-2-2-2
Примеры расчета. В качестве примера рассмотрим систему со следующими параметрами: / = 2 м; д 0 = 5 кг; ^ = 25 кг; т = 1,5 кг/м; Е1 = 4000 Н • м2; иТ = 3 м; Т = 2 с. Число гасимых собственных форм колебаний принимаем последовательно равным N = 1, 2, 3. При выбранных параметрах во всех расчетных случаях колебания консольного стержня с массой на конце остаются в рамках линейной теории, т. е. максимальный угол поворота поперечных сечений стержня | Щ, г)Ни'& г)| < 0,5 рад.
Собственные частоты колебаний рассматриваемой системы таковы: - = 7,639 рад/с; -2 = 200,490 рад/с; -3 = 646,582 рад/с; -4 = = 1347,372 рад/с.
В первом расчетном случае рассмотрим функцию ф(г) в виде ряда по синусам (9). Коэффициенты Ък имеют следующие значения: для N = 1 Ъ1 = 10,441; Ъ2 = -2,032; для N = 2 Ъ1 = 11,743; Ъ2 = -3,654; Ъ3 = -1,474; для N = 3 Ъ1 = 12,526; Ъ2 = -4,872; Ъ3 = -3,368; Ъ4 = 1,830.
На рис. 2 показан вид функций ф(г) в виде ряда (9) при N = 1, 2, 3. Здесь и далее на рисунках сплошная линия на графике соответствует значению N = 1, пунктирная линия — N = 2, штрихпунктирная линия — N = 3. На рис. 3 представлен график зависимости прогиба упругой балки на конце и(1, г), на рис. 4 — график изменения реакции Р0 (г) при конечном перемещении упругой балки с массами Д0 и Д1 на концах.
Следует отметить, что в данном расчетном случае даже при N = 1 остаточные колебания при г > Т конца балки и(1, г) по непогашенным собственным формам п = 2, 3,... составляют менее 10-6 м, остаточная реакция Р0(г) не превышает 0,05 Н.
Во втором расчетном случае функция ф(г) представлена в виде ряда по косинусам (12). Приведем по аналогии значения коэффициентов ак: для N = 1 а1 = 7,975; а3 = -5,158; для N = 2 а1 = 8,307; а3 = -8,055; а5 = -0,248; для N = 3 а1 = 8,480; а3 = -9,867; а5 = = -0,506; а7 = 1,892.
ф(0
Рис. 2. Финитная функция в первом расчетном случае 1) (1,0.м
Рис. 4. Реакция системы в первом расчетном случае
На рис. 5 показан вид функций ф(г) для N = 1, 2, 3, на рис. 6 — график зависимости прогиба упругой балки на конце и(1, г), на рис. 7 — реакция /0(г) при конечном перемещении системы. В дан-
ном расчетном случае при N = 1 остаточные колебания при I > Т конца балки и(1,1) по непогашенным собственным формам п = 2, 3,...
составляют менее 10"5 м, остаточная реакция ) не превышает 0,001 Н.
Рис. 5. Финитная функция во втором расчетном случае
-о (1,0, м
Р0 (О, Н
Рис. 7. Реакция системы во втором расчетном случае
В рассмотренном примере Т1 = 2л/- < Т, -2/-1 ~ 26, и поэтому низшая собственная форма оказывает доминирующее влияние на нестационарные колебания при передвижении балки за время Т. В этом случае при гашении колебаний только по одной низшей собственной форме п = N = 1 остаточные колебаний по непогашенным высшим формам п = 2, 3, . оказываются весьма малыми.
Выводы.
1. Решена задача о гашении колебаний в момент остановки при конечном передвижении упругого стержня с массой на конце. Получены решения данной задачи в рядах по синусам и в рядах по косинусам. При гашении колебаний по N низшим собственным формам задача сводится к системе N +1 линейных уравнений для коэффициентов этих рядов.
2. Рассмотрены примеры расчета остаточных колебаний системы после остановки ее передвижения при гашении различного числа ^ = 1, 2, 3) низших собственных форм колебаний, эти колебания оценены.
Работа выполнена в рамках государственного задания (номер государственной регистрации темы AAAA-A17-117032010144-8) и при частичной финансовой поддержке Российского фонда фундаментальных исследований (проект № 18-08-00778а).
ЛИТЕРАТУРА
[1] Черноусько Ф.Л., Болотник Н.Н., Градецкий В.Г. Манипуляционные роботы: динамика, управление, оптимизация. Москва, Наука, 1989, 363 с.
[2] Бербюк В.Б. Динамика и оптимизация робототехнических систем. Киев, Наукова Думка, 1989, 187 с.
[3] Черноусько Ф.Л., Акуленко Л.Д., Соколов Б.Н. Управление колебаниями. Москва, Наука, 1976, 383 с.
[4] Черноусько Ф.Л., Ананьевский И.М., Решмин С. А. Методы управления нелинейными механическими системами. Москва, Физматлит, 2006, 326 с.
[5] Гришанина Т.В. Управляемый поворот упругого стержня на конечный угол. Вестник МАИ, 2004, т. 11, № 1, с. 64-68.
[6] Гришанина Т.В. Устранение колебаний упругой системы после ее быстрого передвижения и поворота. Вестник МАИ, 2004, т. 11, № 2, с. 68-75.
[7] Гришанина Т.В. Динамика управляемого движения упругих систем при конечных перемещениях и поворотах. Изв. РАН. МТТ, 2004, № б, с. 171-186.
[8] Гришанина Т.В., Шклярчук Ф.Н. Динамика упругих управляемых конструкций. Москва, Изд-во МАИ, 2007, 328 с.
[9] Русских С.В. Управляемый поворот космического аппарата с упругими панелями солнечных батарей. Известия высших учебных заведений. Машиностроение, 2016, № 12, с. 97-105.
[10] Гришанина Т.В., Русских С.В., Шклярчук Ф.Н. Управление конечным поворотом упругой системы из одного состояния в другое с гашением коле-
баний в момент окончания операции. Ученые записки Казанского университета. Серия физико-математические науки, 2017, т. 159, кн. 4, с. 429-443.
Статья поступила в редакцию 20.04.2018
Ссылку на эту статью просим оформлять следующим образом: Русских С.В., Шклярчук Ф.Н. Конечное поперечное передвижение упругого стержня с массой на конце при гашении колебаний в момент остановки. Инженерный журнал: наука и инновации, 20187, вып. 7. http://dx.doi.org/10.18698/2308-6033-2018-7-1786
Русских Сергей Владимирович — канд. физ.-мат. наук, доцент кафедры «Проектирование и прочность авиационно-ракетных и космических изделий» МАИ (НИУ), старший научный сотрудник ИПРИМ РАН. Автор 35 опубликованных научных работ. Область научных интересов: механика и динамика деформируемого твердого тела, динамика упругих управляемых конструкций, динамика космических систем. e-mail: sergey.russkih@rambler.ru
Шклярчук Федор Николаевич — д-р техн. наук, профессор, заслуженный деятель науки Российской Федерации. Главный научный сотрудник ИПРИМ РАН, профессор кафедры 602 «Проектирование и прочность авиационно-ракетных и космических изделий» МАИ (НИУ). Автор 215 опубликованных научных работ. Область научных интересов: механика и динамика деформируемого твердого тела, динамика упругих управляемых конструкций, динамика космических систем. e-mail: shklyarchuk@list.ru
Final transverse movement of an elastic rod with mass at the end when the vibrations are damped at the moment of stopping
© S.V. Russkikh1,2, F.N. Shklyurchuk1, 2
1Moscow Aviation Institute (National Research University), Moscow, 125993, Russia 2Institute of Applied Mechanics of Russian Academy of Sciences, Moscow, 125040, Russia
The article considers nonstationary transverse vibrations of an elastic cantilever rod with a mass at the right end under controlled displacement of the left end along a vertical guide. The flexural vibrations of the rod are sought in the form of eigenmode vibration expansion. The following problem is posed: to find the law of controlled movement of the left end of the rod, governing its movement to the required distance in a given time, while simultaneously several lower eigenforms of the elastic vibrations are damped at the moment of stopping. The acceleration of the left end of the rod is considered to be an unknown function. The acceleration is supposed to be proportional to some finite function. This function, in turn, is written in the form of a series of sines or cosines with unknown coefficients at a given time of controlled motion. Examples of calculation with damping of elastic vibrations in the first three eigenforms are given.
Keywords: elastic cantilever rod, normal coordinates, damping of elastic vibrations, finite function, finite displacements
REFERENCES
[1] Chernousko F.L., Bolotnik N.N., Gradetsky V.G. Manipulyatsionnye raboty: dinamika, upravlenie, optimizatsiya [Manipulation robots: dynamics, management, and optimization]. Moscow, Nauka Publ., 1989, 363 p.
[2] Berbuk V.B. Dinamika i optimizatsiya robototekhnicheskikh system [Dynamics and optimization of robotic systems]. Kiev, Naukova dumka Publ., 1989, 187 p.
[3] Chernousko F.L., Akulenko L.D., Sokolov B.N. Upravlenie kolebaniyami [Control of vibrations]. Moscow, Nauka Publ., 1976, 383 p.
[4] Chernousko F.L., Ananyevsky I.M., Reshmin S.A. Metody upravleniya nelineynymi mekhanicheskimi sistemai [Methods of nonlinear mechanical system control]. Moscow, Fizmatlit Publ., 2006, 326 p.
[5] Grishanina TV. Vestnik Moskovskogo aviatsionnogo instituta (Scientificjournal of Moscow Aviation Institute), 2004, vol. 11, no. 1, pp. 64-68.
[6] Grishanina TV. Vestnik Moskovskogo aviatsionnogo instituta (Scientific journal of Moscow Aviation Institute), 2004, vol. 11, no. 2, pp. 68-75.
[7] Grishanina TV. Izvestiya RAN. Mekhanika tverdogo tela — Proceedings of the Russian Academy of Sciences. Mechanics of Solids, 2004, no. 6, pp. 171-186.
[8] Grishanina TV., Shklyarchuk F.N. Dinamika uprugikh upravlyaemykh kon-struktsiy [Dynamics of elastic controllable structures]. Moscow, MAI Publ., 2007, 328 p.
[9] Russkikh S.V. Izvestiya vysshikh uchebnykh zavedeniy. Mashinostroenie — Proceedings of Higher Educational Institutions. Маchine Building, 2016, no. 12, pp. 97-105.
[10] Grishanina TV., Russkikh S.V., Shklyarchuk F.N. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki — Proceeding of Kazan University. Physics and Mathematics Series, 2017, vol. 159, book 4, pp. 429-443.
Final transverse movement of an elastic rod with mass at the end when the vibrations...
Russkikh S.V., Cand. Sc. (Phys. & Math.), Assoc. Professor, Department of Design and Strength of Aeronautical Missile and Space Products, Moscow Aviation Institute (National Research University), Senior Research Fellow, Institute of Applied Mechanics of Russian Academy of Sciences. Author of 35 research publications. Research interests: mechanics and dynamics of a deformable rigid body, dynamics of elastic controllable structures, dynamics of space systems. e-mail: sergey.russkih@rambler.ru
Shklyarchuk F.N., Dr. Sc. (Eng.), Professor, Honored Scientist of the Russian Federation, Chief Research and Development Officer, Institute of Applied Mechanics of Russian Academy of Sciences, Professor of the Department of Design and Strength of Aeronautical Missile and Space Products, Moscow Aviation Institute (National Research University). Author of 215 research publications. Research interests: mechanics and dynamics of a deformable rigid body, dynamics of elastic controllable structures, dynamics of space systems. e-mail: shklyarchuk@list.ru