УДК 534.015.1 Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2 МБС 70J99
Динамика двойного маятника с вязким трением в шарнирах. I. Математическая модель движения и построение диаграммы режимов
А. С. Смирнов1'2, И. А. Кравчинский2'3
1 Санкт-Петербургский политехнический университет Петра Великого, Российская Федерация, 195251, Санкт-Петербург, ул. Политехническая, 29
2 Институт проблем машиноведения Российской академии наук, Российская Федерация, 199178, Санкт-Петербург, Большой пр. В. О., 61
3 ООО «Центр технологического консалтинга»,
Российская Федерация, 195277, Санкт-Петербург, Большой Сампсониевский пр., 28/2
Для цитирования: Смирнов А. С., Кравчинский И. А. Динамика двойного маятника с вязким трением в шарнирах. I. Математическая модель движения и построение диаграммы режимов // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2. С. 371-384. https://doi.org/10.21638/spbu01.2024.210
В работе обсуждаются вопросы динамического поведения двойного математического маятника с идентичными параметрами звеньев и концевых грузов, который испытывает действие вязкого трения в обоих своих шарнирах с различными в общем случае диссипативными коэффициентами. Дается линейная математическая модель движения системы при малых отклонениях и выводится характеристическое уравнение, содержащее два безразмерных диссипативных параметра. Для случая малого демпфирования находятся приближенные аналитические выражения, позволяющие оценить и сопоставить друг с другом факторы затухания при движении системы по каждой из форм колебаний. Строится диаграмма диссипативных режимов движения, которая возникает при разделении плоскости безразмерных параметров дискриминантными кривыми на области с качественно различным характером движения системы. Отмечается, что в рассматриваемой системе может иметь место диссипативный внутренний резонанс и устанавливаются условия его существования в аналитической форме, а также дается их графическая иллюстрация. Статья является первой частью исследования динамики диссипативного двойного маятника, продолжение которого будет представлено в виде отдельной работы «Динамика двойного маятника с вязким трением в шарнирах. II. Диссипативные формы колебаний и оптимизация параметров демпфирования».
Ключевые слова: двойной маятник, вязкое демпфирование, малое трение, фактор затухания, дискриминантная кривая, диаграмма диссипативных режимов, диссипатив-ный внутренний резонанс.
1. Введение. Динамике и управлению движением многозвенных маятниковых систем посвящено весьма большое количество научных работ, что связано с их широким использованием в прикладных задачах робототехники и биомеханики как элементов разнообразных манипуляторов и многочисленных биоморфных устройств [1-5]. При этом важную роль играет вопрос об учете сил сопротивления
© Санкт-Петербургский государственный университет, 2024
в шарнирных сочленениях подобных конструкций, которые специально вводятся для подавления их колебаний, и влиянии этих сил на динамическое поведение и характер движения рассматриваемых систем [6-8]. На этом пути естественным образом возникает необходимость определения, построения и детального исследования форм колебаний в диссипативных системах, которые обобщали бы известное понятие форм колебаний в консервативных системах. Указанным аспектам уделяется явно недостаточное внимание в существующей литературе, а в классических трудах и справочниках по теории механических колебаний данная тема освещается крайне скупо лишь в общих чертах и без каких-либо конкретных числовых примеров, которые прояснили бы ключевые моменты и тонкости, возникающие при построении диссипативных форм колебаний [9-13].
Ввиду того, что наиболее простой многозвенной конструкцией является двойной математический маятник с идентичными звеньями и концевыми грузами, вполне логично именно на этом весьма распространенном примере продемонстрировать последовательно все шаги, необходимые для построения диссипативных форм его колебаний при наличии вязкого трения в обоих шарнирах. Сложность предлагаемого исследования состоит в том, что имеется два независимых коэффициента вязкого демпфирования, комбинации которых могут быть совершенно различными, в отличие от работ [14-16], где диссипативные параметры в шарнирах считались для упрощения одинаковыми, что позволило получить конечные результаты аналитическим путем без особых затруднений. Поэтому представляет несомненный интерес установить зависимость всех основных величин, описывающих диссипативные формы колебаний, от коэффициентов демпфирования, обсудить их качественные и количественные характеристики, а также выявить многие нетривиальные особенности, которые требуют серьезного внимания.
Кроме того, отдельного разговора заслуживают и попутно возникающие задачи оптимизации процессов затухания движений двойного маятника, в которых также можно распоряжаться выбором наилучшего по некоторому критерию качества значения не одного лишь диссипативного параметра, когда трение в обоих шарнирах двойного маятника полагается идентичным, а сразу двух коэффициентов вязкого демпфирования, принимая их оба в качестве оптимизируемых параметров. Это обстоятельство вновь существенно усложняет исследование, однако посредством надлежащего подбора указанных коэффициентов оно позволяет обеспечить наиболее эффективное подавление колебаний системы и выявить тем самым предельные возможности в отношении пассивного гашения колебаний [17, 18].
Предлагаемая работа, публикуемая в двух частях, охватывает все указанные вопросы динамического поведения двойного маятника с учетом вязкого трения в обоих его шарнирных сочленениях. Настоящая первая часть работы посвящена изучению математической модели движения диссипативного двойного маятника с рассмотрением сначала варианта малого трения, а затем построением диаграммы диссипа-тивных режимов при произвольных коэффициентах демпфирования и, кроме того, анализом интересного случая диссипативного внутреннего резонанса. Вторая часть работы, которая будет опубликована в виде отдельной статьи, будет посвящена продолжению данного исследования, а именно построению и анализу диссипативных форм колебаний двойного маятника с выявлением условия их независимости от демпфирования и детальным обсуждением данной ситуации, а также вопросам оптимизации параметров вязкого демпфирования.
2. Математическая модель движения диссипативного двойного маятника. Рассмотрим расчетную схему двойного математического маятника, полагая для простоты массы его концевых грузов идентичными и равными т и считая также его длины одинаковыми и равными I. Кроме того, предположим, что этот маятник испытывает действие вязкого трения в неподвижном и подвижном шарнирах с различными в общем случае коэффициентами демпфирования 61 и 62 соответственно (рис. 1).
Ч
Рис. 1. Расчетная схема двойного маятника с вязким трением в шарнирах.
т
Построим математическую модель движения этой системы, полагая, что рассматриваются ее малые отклонения. Принимая углы у>1 и у>2, образуемые звеньями маятника с вертикалью, в качестве обобщенных координат, запишем кинетическую и потенциальную энергии двойного маятника в квадратичной аппроксимации, которые определяются хорошо известными выражениями [14]:
Т = ^тп12 (2ф\ + 2фхф2 + ф1) = ф,
П
^тд1 (2^+^) = ^<ртС<р,
(1)
где ф = [у>1,у>2 ]Т — столбец обобщенных координат, а А и С — матрицы инерционных и квазиупругих коэффициентов линейной модели соответственно:
А = тР
2 1 1 1
С = тд1
2 0 0 1
(2)
Для учета влияния вязкого трения в шарнирах двойного маятника запишем дисси-пативную функцию Рэлея в виде [19]
к = ¿(М? + М22) = \ [(Ь1 + Ъ2)ф\ - 2Ъ2ф1ф2 + Ъ2ф\} = ^фТВф,
(3)
где 01 = , 62 = ^2 — — шарнирные углы поворота, а матрица диссипативных коэффициентов В имеет вид
В
61 + 62 -62
-62
62
(4)
Следовательно, уравнения движения диссипативного двойного маятника можно записать в известной матричной форме [9]:
Аф + Вф + Сф = 0.
(5)
Разыскивая решение этого матричного уравнения в виде ф = Фвхг, приходим к алгебраическому матричному уравнению относительно столбца Ф:
(ЛЛ2 + БЛ + С)Ф =0, (6)
для которого условие существования нетривиального решения имеет вид
ёе^ЛЛ2 + БЛ + С)} = 0. (7)
Из этого характеристического уравнения можно определить собственные значения диссипативного двойного маятника. Подставляя в (7) выражения для матриц Л, Б и С согласно (2) и (4), после раскрытия определителя и ряда преобразований придем к уравнению четвертой степени относительно собственных значений Л:
аоЛ4 + ахЛ3 + а2Л2 + азЛ + а4 = 0, (8)
коэффициенты которого имеют следующую структуру:
а0 = т214, а1 = т12(Ъ1 + 5Ъ2), а2 = 4т2д13 + Ъ1Ъ2, аз = тд1(Ъ1 + 3Ъ2), а4 = 2т2д212.
Видно, что все коэффициенты (9) положительны, а диссипативные параметры Ъ1 и Ъ2 входят линейно только в те коэффициенты уравнения (8), которые стоят перед нечетными степенями Л и которые обращаются в нуль при Ъ1 = Ъ2 = 0 [18]. Для упрощения дальнейшего исследования преобразуем характеристическое уравнение к более удобной форме:
Я(р) = Р4 + 2(^1 + 5^2 )р3 + 4(1 + ^)р2 + 2(^1 + 3v2)p + 2 = 0, (10) где введены следующие обозначения:
/ Гв Л Ь2
Поясним эти обозначения: к — частота малых колебаний обычного маятника длиной I, которая необходима для обезразмеривания; р — безразмерный корень характеристического уравнения системы; а параметры V! и ^2 можно назвать безразмерными диссипативными коэффициентами, причем обезразмеривание величин Ъ1 и Ъ2 осуществлено однотипным образом. Решая уравнение (10) при помощи численных процедур, можно найти его четыре корня р8, в = 1,4, при любых значениях параметров V1 и V2.
3. Исследование случая малого трения. Начнем с рассмотрения наиболее простого случая малого трения, когда Vl и V2 являются малыми величинами. Отметим сначала, что в рамках консервативной модели, т.е. при Vl = V2 = 0, частоты колебаний двойного маятника кю и к2о имеют следующие значения [14]:
к10 = \J2-V2k, к20 = ^2 + \F~lk. (12)
Нетрудно понять, что при малых значениях диссипативных параметров Vl и V2 частоты колебаний двойного маятника с учетом вязкого трения практически не изменятся по сравнению с консервативным вариантом (12) [20]. Тем не менее сколь
угодно малое трение вызывает затухание колебаний, наблюдаемое на практике, поэтому основными характеристиками такого процесса будут являться факторы затухания каждой из форм колебаний, определяемые модулем вещественной части каждой из пар корней характеристического уравнения (10), которые при малых VI и ^2, очевидно, являются комплексно-сопряженными. Обращаясь к этому уравнению, определим его решения с точностью до первого порядка малости по величинам VI и V2, используя асимптотические методы решения алгебраических уравнений [21]. Ясно, что решения невозмущенного характеристического уравнения (при отсутствии диссипации, т.е. при VI = V2 = 0) имеют вид
Р1,2 = ±«^10, Р3,4 = ±¿^20, (13)
где сою = кю/к = у2 — \/2 и и>2о = к2о/к = + — безразмерные частоты колебаний двойного маятника, вычисленные согласно (12). Будем разыскивать решения уравнения (10) при ненулевых значениях VI и V2 в виде
Р1,2 = -П1 ± Шю, Р3,4 = -П2 ± Ш20, (14)
где величины П1 и П2 представим линейными комбинациями параметров V! и V2:
П1 = + П2 = a2Vl + ^2, (15)
причем коэффициенты «1, в1, «2 и в подлежат определению после подстановки (14) в уравнение (10) и выполнения его с требуемой точностью. В результате указанных действий можно получить следующие значения:
2-^2 2 + а/2 10-7^2 10 + 7^2
«1 = —1—, «2 = —^—, Р\ = -^-, Р 2 =-^-• (16)
Видно, что эти величины оказываются вещественными, т. е. поправка первого порядка малости действительно входит только в вещественные части решений (14), которая нас и интересует, тогда как поправка в мнимую часть будет иметь, очевидно, второй порядок малости. Эти выводы подтверждают, что частоты колебаний при малых значениях Vl и V2 могут быть приняты такими же, как и в консервативной модели. Таким образом, окончательные приближенные выражения для величин П1 и П2, представляющих собой искомые факторы затухания, получаются в результате подстановки (16) в (15):
2~ У2 , Ю-7У2 2 + У2 , 10 + 7У2
п\ = ---VI Л----V21 п2 = ---VI Н----V2. (17)
4 4 4 4
Чтобы сопоставить эти значения, образуем отношение величин П2 и П1, вводя при этом в рассмотрение безразмерный параметр п = v2/v1:
, = п2 = 2 + а/2 + (10 + 7л/2)г] п\ 2 — а/2 + (10 — 7А/2)?У '
Подчеркнем, что хотя значения Vl и V2 предполагаются малыми, параметр п может принимать любые значения в диапазоне 0 < п < Видно, что на левой границе этого интервала, когда г/ = 0 {и2 = 0, т. е. в подвижном шарнире трение отсутствует), имеем / = 3 + 2 а/2 « 5.828, тогда как на правой границе этого интервала,
150 / 100 50 0
0 10 20 30 40 50 60 70 80 90 100
1
Рис.2. График зависимости величины f = пч/пх от параметра п = ич/их.
когда г/ = оо (г/1 = 0, т.е. в неподвижном шарнире трение отсутствует), имеем / = 99 + 70\/2 « 197.99. Наконец, в случае 1] = 1 (г/1 = г/2, т.е. в обоих шарнирах трение одинаковое) будем иметь / = 17 + 12а/2 « 33.971. Таким образом, видно, что в зависимости от параметра п = V2/vl движения системы по первой и второй формам колебаний при малом трении могут гаситься совершенно различно в количественном плане по отношению друг к другу. Тем не менее ясно, что П2 > «1, а потому вторая форма всегда будет подавляться быстрее, чем первая, т. е. в общем случае задания начальных условий движения будут происходить достаточно быстрое исчезновение колебаний по второй форме и их длительное затухание по первой форме (конечно, формально полное затухание происходит лишь за бесконечное время). График дробно-рациональной функции /(п) согласно (18) для наглядности приведен на рис. 2.
4. Построение диаграммы диссипативных режимов движения.
В предыдущем пункте было показано, что при малых величинах Vl и V2 корни уравнения (10) являются комплексно-сопряженными. Ясно, что такая ситуация будет сохраняться и при достаточно небольших значениях этих параметров, однако при их постепенном увеличении два корня этого уравнения уже станут вещественными и отрицательными, тогда как два других корня по-прежнему останутся комплексно-сопряженными, а если Vl и V2 вовсе будут являться весьма значительными, то все четыре корня окажутся вещественными и отрицательными. Чтобы разобраться, при каких значениях данных параметров будет реализовываться тот или иной вариант, перейдем к построению и обсуждению диаграммы диссипативных режимов движения двойного маятника на плоскости безразмерных параметров Vl и V2. Поскольку оба эти параметра положительны, то нас будет интересовать лишь первый квадрант указанной плоскости. Разделим его на области с качественно разным характером движения, чему отвечает наличие различного числа вещественных корней характеристического уравнения (10). Ясно, что границам этих областей будут отвечать те кривые на плоскости VlV2, для которых уравнение (10) имеет кратные корни. Чтобы установить эти кривые, которые мы будем называть дискриминант-ными, необходимо решить систему уравнений ^(р) = 0, Я'(р) = 0, которая в явном виде будет:
Q(p) = Р4 + 2^1 +5v2)p3 +4(1 + VlV2)p2 +2^1 +3v2)p + 2 = 0,
(19)
Q'(p) = 4р3 + 6^1 + 5v2 )р2 + 8(1 + Vl V2)p + 2^1 + 3V2) = 0.
Выразим из второго уравнения (19), например, VI:
2р(р2 + 2) + 3(5р2 + 1)г/2 (2(])
3р2 + 1 + 4^2
и подставим это соотношение в первое уравнение (19), в результате чего после упрощений получим:
4р2(5р2 - 3^2 + 4р(р4 + 2р2 - 2)^ + р6 - р4 - 2р2 - 2 = 0. (21)
Решения этого квадратного уравнения относительно V2 имеют вид
р4 + 2р2 - 2 ± (2р2 + 1)^-р4 + 4р2-2 -2р(3-5р2)-' (22)
Подставляя эти выражения в формулу (20), получим после серии несложных преобразований:
р4 - 2р2 - 2 Т (2р2 + 1)У-р4 + 4р2-2 1/1 = -2р(1-р2)-' (23)
причем в соотношениях (22) и (23) верхние и нижние знаки соответствуют друг другу. Уравнения (22) и (23) представляют собой параметрические уравнения искомых дискриминантных кривых, отвечающих случаю наличия кратных корней характеристического уравнения, причем величина р, которая, конечно, является вещественной и отрицательной, играет роль параметра. Для удобства последующего построения этих кривых запишем их далее по отдельности. Так, беря нижние знаки в формулах (22) и (23), получим:
. . р4 -2р2 -2 + (2р2 + 1 )у/-р4 + Ар2 -2 =--'
Р4 + 2Р2 - 2 -
р4 + 2р2 - 2- (2р2 + 1 )л/-р4 + Ар2 -2
(24)
2р(3 - 5р2)
а беря в этих формулах верхние знаки, будем иметь:
. . р4 -2р2 - 2- (2р2 + 1)л/-р4 + 4р2 -2
=-ы^р1)-'
. , р4 + 2р2 -2 + (2р2 + 1)л/-р4 + 4р2 -2 ^ =-2р(3-5р2)-'
(25)
Очевидно, что для вещественности величин VI и V2 в обоих случаях следует потребовать положительность подкоренного выражения, т. е. выполнение неравенства -р4 + 4р2 - 2 > 0, откуда ввиду отрицательности р вытекает, что варьировать этот параметр необходимо в следующих пределах:
-1.848 « — у 2 + а/2 < р < —у 2 — \/2 к, -0.765. (26)
Далее рассмотрим по отдельности каждую из зависимостей (24) и (25) и разберемся с тем, как они будут себя вести при уменьшении р в диапазоне (26) от верхней
01234567
vi
Рис. 3. Диаграмма диссипативных режимов движения на
плоскости Ь>1Ь>2-
границы до нижней. При р = —\/2 — а/2 « —0.765 зависимости (24) и (25) дают одну и ту же точку ь>\ = щ = л/10 + 7 а/2 « 4 .461, лежащую на биссектрисе первого квадранта. Далее при уменьшении р от — \/2 — а/2 « —0.765 до — у^З/5 к, —0.775, т. е. в весьма узких пределах, кривые (24) и (25) расходятся от указанной выше точки, причем при р —> —\/3/5 + 0 для кривой (24) будем иметь v\ —> a/Ï5 « 3.873, а v2 —> +оо, тогда как для кривой (25) имеем в этих условиях ь>\ —> 41/(2a/Ï5) « 5.293, а для щ после раскрытия неопределенности вида 0/0 имеем щ —> 19/(2a/Ï5) « 2.453. Дальнейшее уменьшение р приводит к тому, что кривая (24) оказывается вне первого квадранта и попадает в него лишь при достижении значения рФ, подлежащего определению. Что же касается кривой (25), то она в указанных условиях продолжается от точки v\ = 41/(2 а/15), vo = 19/(2a/Ï5), причем при р —> —1 + 0 будем иметь —>■ +оо, а щ —^1, а при дальнейшем понижении р она также оказывается вне первого квадранта вплоть до того, как р примет значение р**, также пока неизвестное. Нетрудно понять, что кривая (24) попадет в первый квадрант в результате пересечения с осью vo = 0, а кривая (25) окажется в этом квадранте вследствие пересечения с осью v\ = 0. Приравнивая нулю выражение для щ из (24) и выражение для из (25), можно установить соответственно, что р* = —1.507, при этом v\ = 1.647, а р** = —1.837, при этом щ = 0.368. Наконец, будем уменьшать р от р, = —1.507 до + а/2 « -1.848 в формулах (24), а также уменьшать р от = —1.837 до + а/2 « -1.848 в формулах (25), вследствие чего получим оставшиеся участки кривых, расположенные в первом квадранте. Отметим, что эти участки при р =
-V2+V2
стыкуются в точке v\ = vo = VlO - 7V2 « 0.317,
которая лежит на биссектрисе первого квадранта, как и та точка, из которой мы начали построение кривых прир = —л/2 — л/2. На рис. 3 участки кривых, отвечающих формулам (24), построены синим цветом, а участки кривых, отвечающих формулам (25) — красным цветом. Объединяя попарно участки, стыкующиеся друг с другом на биссектрисе первого квадранта (которая представлена пунктирной прямой), мы получаем две результирующие кривые, разделяющие первый квадрант плоскости VIV2 на три области, классификация которых такова: в области I имеются две пары комплексно-сопряженных корней, в области II — пара комплексно-сопряженных корней и два отрицательных корня, и, наконец, в области III — четыре отрицательных корня. Это означает, что обе диссипативные формы колебаний существуют только в области I, в области II будет иметь место только одна диссипативная форма колебаний, тогда как в области III диссипативных форм колебаний уже не существует. Разумеется, данные заключения носят лишь качественный характер, тогда как количественный анализ диссипативных форм колебаний представляет собой отдельную и весьма трудоемкую задачу.
5. Условие внутреннего резонанса. Обратимся к более детальному анализу области I на диаграмме диссипативных режимов движения, в которой имеются две пары комплексно-сопряженных корней. Выясним, возможно ли для рассматриваемой диссипативной системы существование случая диссипативного внутреннего резонанса, означающего совпадение частот ее колебаний при движении по каждой из форм колебаний. Этот случай важен не только с теоретической точки зрения, но также может представлять и определенный практический интерес, так как движение двойного маятника в общем случае будет являться одночастотным. С этой целью положим, что корни уравнения (10) есть
р1,2 = -П1 ± '1ш, рз,4 = -П2 ± гш, (27)
т. е. все они имеют одинаковую по абсолютной величине мнимую часть, причем ш и есть общая частота колебаний. Тогда характеристический полином Я(р) можно записать в следующем виде:
Я(р) = (р + п1 - гш)(р + П1 + гш)(р + П2 - гш)(р + П2 + гш) =
= р4 + 2(п1 + п2)р3 + (п2 + п2 + 2ш2 + 4п1п2)р2+ (28)
+2 [щ(п2 + ш2) + п2(п1 + ш2)] р + (п2 + ш2)(п2 + ш2).
Приравнивая коэффициенты при одинаковых степенях в выражениях (10) и (28), получим следующую систему уравнений:
П1 + П2 = V! + 5V2,
п2 + п2 + 2ш2 +4П1П2 =4(1 + V1V2),
(29)
п1(п2 + ш2) + п2(п1 + ш2) = V1 +3V2, (п1 + ш2)(п2 + ш2) = 2. Переписывая третье из этих уравнений в виде
(п1 + п2)(п1п2 + ш2) = V1 + 3V2 (30)
и исключая отсюда и\ + посредством первого уравнения (29), можно выразить комбинацию П1П2:
VI + 3V2 2
П1П2
VI + 5v2
Обращаясь теперь ко второму уравнению (29) и записывая его в виде (п1 + П2)2 + 2П1П2 + 2ш2 = 4(1 + V1V2),
(31)
(32)
вновь исключим из него П1 + П2 посредством первого уравнения (29) и П1П2 посредством уравнения (31), в результате чего будем иметь:
V + 5V2)3 + 2^1 + 3V2) = 4(1 + V1V2)(V1 + 5V2).
(33)
Отметим, что при этом слагаемые с ш сокращаются, так что уравнение (33) и есть неявное уравнение искомой кривой на плоскости VlV2, отвечающей случаю внутреннего резонанса, которая, очевидно, должна располагаться внутри области I. Уравнение (33) также может быть переписано в альтернативной форме:
V? + 125 VI + + 5V2) = 2^1 + 7V2).
(34)
Отсюда вытекает, что при V! =0 будем иметь V2 = л/14/125 « 0.335, а при V2 =0 — VI = а/2 « 1.414. Подчеркнем также, что последнее из уравнений системы (29) нами не было использовано при получении уравнения (33). Перепишем это последнее уравнение (29) в виде
(П1П2)2 + [(П1 + П2)2 - 2П1П^ ш2 + (
(35)
и подставим в него П1 + П2 согласно первому уравнению (29) и П1П2 согласно (31), и в результате несложных преобразований получим биквадратное уравнение относительно ш:
4ш4 +
V + 5V2)2 - 4
VI + Зг/2
Vl + 5V2
ш2 +
VI + Зг/9 Vl + 5V2
-2 = 0.
(36)
Рассматривая его совместно с уравнением (33), можно при необходимости установить, как изменяется частота ш при движении по кривой, отвечающей случаю внутреннего резонанса. Остается лишь построить область I в увеличенном масштабе и провести внутри нее кривую (33) (она изображена зеленым цветом), которая разделит эту область на две подобласти I, а и I, Ь (рис. 4). Очевидно, что в подобласти I, а,
0.4 0.3
V 0.2 0.1 0
/ /.......-ч /
/ / I, а
/ / I, ь
0.2 0.4 0.6 0.8
1.2 1.4 1.6 1.8
Рис. 4. Кривая, отвечающая случаю внутреннего резонанса, в области I на плоскости
4
2
V
которая примыкает к точке vi = V2 =0 (отвечающей консервативному варианту) и содержит внутри себя зону малых значений vi и V2, подробно изученную ранее, частота колебаний при движении по первой форме колебаний диссипативной системы будет меньше частоты колебаний при движении по второй форме колебаний, и эта ситуация является типичной. Что же касается подобласти I, b, то в ней будет реализовываться противоположная ситуация, которая является весьма нетривиальной. Эти заключения лишний раз подчеркивают, что установление характера диссипативных форм колебаний требует весьма скрупулезного анализа, чему и будет посвящено дальнейшее исследование.
6. Заключение. В настоящей работе была рассмотрена динамика двойного математического маятника, обладающего идентичными звеньями и концевыми грузами, с вязким трением в обоих его шарнирных сочленениях. Для случая малых значений диссипативных коэффициентов были получены приближенные аналитические выражения для факторов затухания каждой из форм колебаний и показано, что вторая форма при любых соотношениях между параметрами демпфирования подавляется силами сопротивления намного быстрее, чем первая форма. При рассмотрении произвольных значений диссипативных коэффициентов была построена диаграмма режимов движения двойного маятника на плоскости безразмерных параметров демпфирования, позволяющая выявить поведение корней характеристического уравнения при изменении этих параметров. Помимо этого, было установлено, что в диссипативной системе возможен внутренний резонанс, отвечающий равенству частот колебаний при движении по каждой из форм колебаний и заслуживающий отдельного внимания ввиду его возможного практического использования. Полученные результаты были приведены в удобной графической форме, позволяющей наиболее детально проследить все обнаруженные качественные и количественные особенности. В следующей статье «Динамика двойного маятника с вязким трением в шарнирах. II. Диссипативные формы колебаний и оптимизация параметров демпфирования» будет рассмотрен непосредственно вопрос о построении и анализе диссипативных форм колебаний двойного маятника на основе разумного сочетания аналитических и численных методов исследования, а также будет поставлено и решено несколько сопутствующих задач об оптимальном пассивном гашении колебаний рассматриваемой системы.
Литература
1. Формальский А. М. Управление движением неустойчивых объектов. Москва, Физматлит (2014).
2. Перегудова О. А., Макаров Д. С. Синтез управления трехзвенным манипулятором. Автоматизация процессов управления 2 (40), 109—113 (2015).
3. Анохин Н. В. Приведение многозвенного маятника в положение равновесия с помощью одного управляющего момента. Известия РАН. Теория и системы управления 5, 44—53 (2013).
4. Смольников Б. А., Юревич Е. И. К проблеме биоморфного управления движениям роботов. Робототехника и техническая кибернетика 1 (6), 17—20 (2015).
5. Тяжелов А. А., Кизилова Н. Н., Фищенко В. А., Яремин С. Ю., Карпинский М. Ю., Карпинская Е. Д. Анализ стабилограмм на основе математической модели тела человека как многозвенной системы. Травма 13 (4), 17-25 (2012).
6. Sawant K. R., Shrikanth V. Energy dissipation and behavioral regimes in an autonomous double pendulum subjected to viscous and dry friction damping. European Journal of Physics 42 (5), 055008 (2021).
7. Bendersky S., Sandler B. Investigation of a spatial double pendulum: An engineering approach. Discrete Dynamics in Nature and Society 2006, 1-22 (2006).
8. Smirnov A. S., Smolnikov B.A. Dissipative Model of Double Mathematical Pendulum with Noncollinear Joints. Lecture Notes in Mechanical Engineering. ¡Selected Contributions from the Conference "Modern Engineering: Science and Education". St. Petersburg, Russia, June 2021, 38—47 (2022).
9. Бидерман В. Л. Теория механических колебаний. Москва, Высшая школа (1980).
10. Тимошенко С. П. Колебания в инженерном деле. Москва, Наука (1967).
11. Магнус К. Колебания: введение в исследование колебательных систем. Москва, Мир (1982).
12. Болотин В. В. (ред.). Вибрации в технике. Справочник. Т. 1. Колебания линейных систем. Москва, Машиностроение (1978).
13. Биргер И. А., Пановко Я. Г. (ред.). Прочность, устойчивость, колебания. Т. 3. Москва, Машиностроение (1968).
14. Леонтьев В. А., Смирнов А. С., Смольников Б. А. Оптимальное демпфирование колебаний двухзвенного манипулятора. Робототехника и техническая кибернетика 2 (19), 52—59 (2018).
15. Смирнов А. С., Смольников Б. А. Оптимизация режимов гашения колебаний пространственного двойного маятника. I. Постановка задачи. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 9 (67), вып. 2, 357—365 (2022). https://doi.org/10.21638/spbu01.2022.215
16. Смирнов А. С., Смольников Б. А. Оптимизация режимов гашения колебаний пространственного двойного маятника. II. Решение задачи и анализ результатов. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 10 (68), вып. 1, 121—138 (2023). https://doi.org/10.21638/spbu01.2023.111
17. Болотник Н.Н. Оптимизация амортизационных систем. Москва, Наука (1983).
18. Смольников Б. А. Проблемы механики и оптимизации роботов. Москва, Наука (1991).
19. Карман Т., Био М. Математические методы в инженерном деле. Москва; Ленинград, ГИТТЛ (1946).
20. Блехман И. И., Мышкис А. Д., Пановко Я. Г. Механика и прикладная математика. Логика и особенности приложений математики. Москва, Наука (1983).
21. Найфэ А. Введение в методы возмущений. Москва, Мир (1984).
Статья поступила в редакцию 7 августа 2023 г.;
доработана 28 октября 2023 г.; рекомендована к печати 9 ноября 2023 г.
Контактная информация:
Смирнов Алексей Сергеевич — ассистент, мл. науч. сотр.; smirnov.alexey.1994@gmail.com Кравчинский Игорь Андреевич — аспирант; igor.krav2000@gmail.com
Dynamics of a double pendulum with viscous friction in the joints. I. Mathematical model of motion and construction of the regime diagram
A. S. Smirnov1'2, I. A. Kravchinskiy2'3
1 Peter the Great St. Petersburg Polytechnic University,
29, ul. Polytechnicheskaya, St. Petersburg, 195251, Russian Federation
2 Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences, 61, Bolshoy pr. V. O., St. Petersburg, 199178, Russian Federation
3 LLC "Change Mobility Together",
28/2, Bolshoy Sampsonievsky pr., St. Petersburg, 195277, Russian Federation
For citation: Smirnov A. S., Kravchinskiy I. A. Dynamics of a double pendulum with viscous friction in the joints. I. Mathematical model of motion and construction of the regime diagram. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11 (69), issue 2, pp. 371-384. https://doi.org/10.21638/spbu01.2024.210 (In Russian)
The paper discusses the dynamic behavior of a double mathematical pendulum with identical parameters of links and end loads, which is under the influence of viscous friction in both of its joints with generally different dissipative coefficients. A linear mathematical model of the system motion for small deviations is given, and a characteristic equation
containing two dimensionless dissipative parameters is derived. For the case of low damping, approximate analytical expressions are found that make it possible to evaluate and compare with each other the damping factors during the motion of the system on each of the oscillation modes. A diagram of dissipative motion regimes is constructed, which arises when the plane of dimensionless parameters is divided by discriminant curves into regions with a qualitatively different character of the system motion. It is noted that a dissipative internal resonance can take place in the system under consideration, and the conditions for its existence in an analytical form are established, as well as their graphic illustration is also given. This article is the first part of the study of the dynamics of a dissipative double pendulum, the continuation of which will be presented as a separate work "Dynamics of a double pendulum with viscous friction in the joints. II. Dissipative oscillation modes and optimization of damping parameters".
Keywords: double pendulum, viscous damping, low friction, damping factor, discriminant curve, dissipative regime diagram, dissipative internal resonance.
References
1. Formalskii A. M. Motion control of unstable objects. Moscow, Fizmatlit Publ. (2014). (In Russian)
2. Peregudova O.A., Makarov D.S. Control synthesis for three-link manipulator. Automation of Control Processes 2 (40), 109-113 (2015). (In Russian)
3. Anokhin N. V. Bringing a multilink pendulum to the equilibrium position using a single control torque. Izvestiia RAN. Teoriia i sistemy upravleniia 5, 44-53 (2013). (In Russian) [Eng. transl.: Journal of Computer and Systems Sciences International 52 (5), 717-725 (2013)].
4. Smolnikov B. A., Yurevich E. I. About the problem of biomorphic robot motion control. Robotics and Technical Cybernetics 1 (6), 17-20 (2015). (In Russian)
5. Tyazhelov A.A., Kizilova N.N., Fischenko V.A., Yaremin S.Yu., Karpinsky M.Yu., Karpinskaya Ye. D. Analysis of posturography based on mathematical model of human body as multilink system. Travma 13 (4), 17-25 (2012). (In Russian)
6. Sawant K. R., Shrikanth V. Energy dissipation and behavioral regimes in an autonomous double pendulum subjected to viscous and dry friction damping. European Journal of Physics 42 (5), 055008 (2021).
7. Bendersky S., Sandler B. Investigation of a spatial double pendulum: An engineering approach. Discrete Dynamics in Nature and Society 2006, 1-22 (2006).
8. Smirnov A. S., Smolnikov B. A. Dissipative Model of Double Mathematical Pendulum with Noncollinear Joints. Lecture Notes in Mechanical Engineering. Selected Contributions from the Conference "Modern Engineering: Science and Education". St. Petersburg, Russia, June 2021, 38-47 (2022).
9. Biderman V.L. Theory of mechanical oscillations. Moscow, Vysshaya shkola Publ. (1980). (In Russian)
10. Timoshenko S. P. Vibration problems in engineering. Toronto; New York; London, D. Van Nostrand Company (1955) [Rus. ed.: Timoshenko S.P. Kolebaniya v inzhenernom dele. Moscow, Nauka Publ. (1967)].
11. Magnus K. Schwingungen. Eine Einführung in die theoretische Behandlung von Schwingungsproblemen. Stuttgart, Teubner (1961) [Rus. ed.: Magnus K. Kolebaniia: vvedenie v issledovanie kolebatel'nykh sistem. Moscow, Mir Publ. (1982)].
12. Bolotin V.V. (ed). Vibrations in technology. Directory. Vol. 1. Oscillations of linear systems. Moscow, Mashinostroenie Publ. (1978). (In Russian)
13. Birger I.A., Panovko Ya. G. (eds). Strength, stability, oscillations. Vol.3. Moscow, Mashinostroenie Publ. (1968). (In Russian)
14. Leontev V. A., Smirnov A. S., Smolnikov B. A. Optimal damping of two-link manipulator oscillations. Robotics and Technical Cybernetics 2 (19), 52-59 (2018). (In Russian)
15. Smirnov A. S., Smolnikov B. A. Optimization of oscillation damping modes of spatial double pendulum. I. Formulation of the problem. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 9 (67), iss. 2, 357-365 (2022). https://doi.org/10.21638/spbu01.2022.215 (In Russian) [Eng. transl.: Vestnik St. Petersburg University, Mathematics 55, iss. 2, 243-248 (2022). https://doi.org/10.21638/spbu01.2022.215].
16. Smirnov A. S., Smolnikov B. A. Optimization of oscillation damping modes of spatial double pendulum. II. Solving the problem and analyzing the results. Vestnik of Saint
Petersburg University. Mathematics. Mechanics. Astronomy 10 (68), iss. 1, 121—138 (2023). https://doi.org/10.21638/spbu01.2023.111 (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mathematics 56, iss. 1, 93-106 (2023). https://doi.org/10.21638/spbu01.2023.111].
17. BolotnikN.N. Optimization of amortization systems. Moscow, NaukaPubl. (1983). (In Russian)
18. Smolnikov B.A. Problems of mechanics and robotoptimization. Moscow, Nauka Publ. (1991). (In Russian)
19. Karman von T., Biot M.A. Mathematical methods in engineering. New York, McGraw-Hill (1940) [Rus. ed.: Karman T., Bio M. Matematicheskie metody v inzhenernom dele. Moscow; Leningrad, GITTL Publ. (1946)].
20. Blekhman 1.1., Myshkis A.D., Panovko Ya. G. Mechanics and applied mathematics. Logic and features of applications of mathematics. Moscow, Nauka Publ. (1983). (In Russian)
21. Nayfeh A. H. Introduction to perturbation techniques. New York, Chichester, Brisbane, Toronto, John Wiley & Sons (1981) [Rus. ed.: Nayfeh A. Vvedenie v metody vozmushchenii. Moscow, Mir Publ. (1984)].
Received: August 7, 2023 Revised: October 28, 2023 Accepted: November 9, 2023
Authors' information:
Alexey S. Smirnov — smirnov.alexey.1994@gmail.com Igor A. Kravchinskiy — igor.krav2000@gmail.com
ХРОНИКА
22 ноября 2023 г. на заседании секции теоретической механики им. проф. Н. Н. Поляхова в Доме ученых им. М. Горького (Санкт-Петербург) выступил кандидат техн. наук, доцент А. А. Суханов (Санкт-Петербургский политехнический университет Петра Великого) с докладом на тему «Об одной модели трения скольжения».
Краткое содержание доклада:
В докладе рассмотрен оригинальный взгляд на природу возникновения трения скольжения между твердыми телами, называемого сухим трением. В отличие от общепризнанной модели сухого трения, заключающейся в контактном взаимодействии и деформировании шероховатостей соприкасающихся поверхностей тел, предложена альтернативная модель, в которой сила трения не является деформационной силой, а представляет собой касательную проекцию силы нормального давления со стороны шероховатостей, выступающих в роли усредненных наклонных плоскостей, вдоль которых «поднимается» одно тело над другим. Предложенная модель сухого трения дает полное теоретическое обоснование эмпирическому закону Амонтона — Кулона, не противоречит результатам сложных относительных движений с ограничениями и, в частности, легко объясняет так называемый парадокс Пенлеве. Проведенные эксперименты наглядно демонстрируют и подтверждают предложенную модель сухого трения.