ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2013 Математика и механика № 1(21)
УДК 534.014
И.П. Попов КОЛЕБАТЕЛЬНЫЕ СИСТЕМЫ, СОСТОЯЩИЕ ТОЛЬКО ИЗ ИНЕРТНЫХ ИЛИ ТОЛЬКО УПРУГИХ ЭЛЕМЕНТОВ, И ВОЗНИКНОВЕНИЕ В НИХ СВОБОДНЫХ ГАРМОНИЧЕСКИХ КОЛЕБАНИЙ
Рассматриваются механические колебательные системы, состоящие только из инертных (гага-, гагага-, ти-системы) или только упругих элементов (кк-, ккк-, к”-системы). Показана возможность возникновения в них свободных гармонических колебаний. В гага-, гагага-, ти-системах происходит взаимный обмен кинетической энергией между инертными элементами, в кк-, ккк-, к1-системах - взаимный обмен потенциальной энергией между упругими элементами.
Ключевые слова: колебательные, инертные, упругие, гармонические, частота.
Свободные гармонические колебания основаны на обмене энергией между элементами колебательной системы.
В механическом линейном гармоническом осцилляторе происходит обмен энергией между разнородными элементами - грузом массой т (инертным элементом) и пружиной с коэффициентом упругости к (упругим элементом). При этом кинетическая энергия груза преобразуется в потенциальную энергию пружины и наоборот.
Существуют электромеханические колебательные системы [1], в которых свободные гармонические колебания осуществляются за счет взаимного преобразования потенциальной энергии пружины в энергию электрического поля конденсатора или кинетической энергии груза в энергию магнитного поля катушки индуктивности.
Таким образом, свободные гармонические колебания сопровождаются самыми разнообразными вариантами преобразования энергии. В этой связи представляет интерес возможность возникновения свободных гармонических колебаний, осуществляемых за счет преобразования кинетической энергии в кинетическую или потенциальной энергии в потенциальную. Реализующие такие колебания системы должны состоять либо только из инертных (тт-, ттт-, т1-системы), либо только из упругих элементов (кк-, ккк-, к”-системы). Механизм обмена энергией между однородными элементами в таких системах позволит, в частности, расширить возможности нейтрализации реакции этих элементов на внешние периодические воздействия.
Как таковая теория колебаний была в основном разработана в 30-х годах 20-го века [2]. Современные источники в части свободных гармонических колебаний преимущественно сохраняют преемственность [3, 4]. В указанных и подобных им источниках колебательные системы с однородными элементами не рассматриваются. В последнее время появляются работы в направлении темы настоящего исследования [5], которое является их развитием.
Синтез mm-системы
Синтез системы осуществляется на основе двух исходных условий.
Первое исходное условие. Система содержит два инертных элемента - два груза массой m каждый. Элементы совершают гармонические колебания x1 = Asin(Z + Cl), x2 = Asin(Z + Z2), где x1, x2 - текущие координаты грузов, A - амплитуда колебаний, Z - фаза, Zb Z2 -начальные фазы.
Второе исходное условие. Энергия системы при колебаниях не меняется
W1 + W2 = const.
Одновременный учет обоих исходных условий дает представление о характере связи между инертными элементами. Действительно,
m ( dx V m ( dx2 V
— I +—\ — I = const,
2 ^ dt ) 2 ^ dt )
cos2(Z + Z1) + cos2(Z + Z2) = const2.
Последнее справедливо при условиях
Zi - Z2 = ±п/2.
Полученное соотношение позволяет определить связующее звено между инертными элементами. Таким звеном является устройство, изображенное на рис. 1.
Рис. 1. mm-система
Анализ тт-системы
Внешние усилия к грузам не приложены. Масса промежуточного стержня и трение не учитываются. Координаты грузов соответственно
х1 = I соБф, х2 = /соє(л/2 - ф). (1)
В качестве обобщенной координаты удобно использовать ф. Система имеет одну степень свободы и уравнение Лагранжа второго рода для нее записывается в виде
а (дТЛ дТ
, , ■ = Q.
ё ^ дф ) дф
Обобщенная сила Q = 0, поскольку активные силы отсутствуют. Кинетическая энергия
^ т (йху )2 т (ёх2 )2 т12 . 2 . 2 т12 2 . 2 т12 . 2
=~2 [) + "2 [) = ~81п фСф + ~008 фСф = ~ф .
К = 0, К = т/2ф, 4“) = т/2ф = 0.
дф дф ё ^ дф)
Решение последнего уравнения
ёф/Ж = С1, ф = + С2. (2)
Пусть начальные условия
ф(0) = ф0, (0) = ®0. (3)
ш
Тогда С2 = ф0, С1 = ©0. (4)
При этом (1) принимает вид
х1 = /со8(ю0/ + ф0), х2 = /со8(л/2 - Ю0/ - ф0). (5)
Пусть
1) х1(0) = х10, со8ф0 = х10/1, ф0 = агссо8(х10/1) = агс81п(х20//);
ёх2
2) ——(0) = у20 , /©0со8(©00 + ф0) = ^20, ©0 = v2o/Xlo = -Vlo/x2o ,
т
х1 = I со8[(у20/х10)/ + агссо8(х10//)], х2 = I со8[п/2 - (-у10/х20)/ - агс81п(х20//)], (6)
VI = /(Ую/х20)81п[(-У10/х20)^ + агс81и(хю//)], у2 = /(v20/x10)co8[(v20/x10)t + агссо8(х20//)].
Таким образом, грузы массой т совершают свободные гармонические колебания (внешние усилия к грузам не приложены).
В рассмотренной колебательной системе происходит взаимный обмен кинетической энергией между инертными элементами. При ф = 0 кинетическая энергия первого груза равна нулю, а второго - максимальна. После этого первый груз начинает ускоряться за счет энергии второго груза, который приобретает отрицательное ускорение.
Синтез mmm-системы
Пусть три координатные оси Ox1, Ox2, Ox3 лежат в одной плоскости Z, последовательно повернуты относительно друг друга на п/3 и пересекаются в одной точке. Точка пересечения O является началом произвольно направленного вектора R, принадлежащего Z.
Теорема 1. Координаты x1, x2, x3 проекций конца вектора R на оси Ox1, Ox2, Ox3 являются вершинами равностороннего треугольника, размер которого не зависит от направления R.
Доказательство. Координаты проекций
x1 = R cos ф, x2 = R cos(n/3 - ф), x3 = R cos(2n/3 - ф). (7)
Здесь ф - угол между R и Ox1. По теореме косинусов
(x1x2)2 = R2[cos^ + cos2(n/3 - ф) - 2cosфCOs(л/3 - ф)^(п/3)] =
R2
(
COS ф +
1
л/3 .
— COS ф +--------Sin ф
V2 2 J
-2cos ф
л/3 .
1
— COS ф +----------Sin ф
V2 2 J
2 1 2 3 . 2 V3 . 12л/з
COS ф + 4COS ф + 4Sin ф +——COSфSinф-2COS ф-— COSфSinф
(x1x3)2 = R2[cos2ф + cos2(2n/3 - ф) - 2cosфCOs(2л/3 - ф)cos(2л/3)] =
= R
2 Г 1 л/з . ї 2 Г і V3. ї Г 1 її
COS ф + COS ф+ Sin ф - 2cos ф COS ф+ Sin ф I -d
2 2 V J 2 2 V J V 2 J
2 1 2 3.2 л/3 . 1 2
cos ф+^cos ф+4sin ф——cosфSinф-—cos ф+——COSфSinф
л/3
(x2x3)2 = R2[cos2(n/3 - ф) + cos2(2n/3 - ф) -2cos(n/3 - ф)cos(2л/3 - ф)cos(л/3)] =
=R2
=R/
1 у/3. ї Г і л/з. Г1 л/з. Y і л/з. ї г
-cosф^-sinфі +I— cosфН-sinфі -2I-cosфН---------------sinф II— а^фН----sinфі-
1 2 3 2 л/3 1 2 3 2 л/з 1232
-cos ф+-sin ф+-C0SфSinф+— cos ф+-sin ф-cosфsmф+—cos ф—sin ф
442 442 44
==R2 4
Теорема доказана.
Теорема 2. Центр треугольника х1х2х3 совпадает с серединой вектора Я.
Доказательство. Пусть г - середина вектора Я. По теореме косинусов (х^г)2 = Л2[со52ф + 1/4 - 2со5ф(1/2)соєф] = К2/4,
(Х2Г)2 = Л2[со82(л/3 - ф) + 1/4 - 2соє(л/3 - ф)(1/2)со8(п/3 - ф)] = Л2/4,
(хзг)2 = Л2[со82(2л/3 - ф) + 1/4 - 2соє(2л/3 - ф)(1/2)со8(2л/3 - ф)] = Л2/4.
Точка г равноотстоит от точек хх и х2, следовательно, она расположена на прямой, перпендикулярной отрезку хіх2 и проходящей через его середину. Это же справедливо в отношении отрезка хіх3. Таким образом, точка г принадлежит двум высотам треугольника ххх3, следовательно, она лежит на их пересечении, которое для равностороннего треугольника является центром. Теорема доказана.
Теоремы 1 и 2 позволяют определить конфигурацию дадада-системы, упрощённая схема которой показана на рис. 2.
Рис. 2. гагага-система
Анализ mmm-системы
Внешние усилия к грузам не приложены. Скорости грузов с учетом (7) dx2/dt = -Rsincp dcp/dt, dx2/dt = Rsin(n/3 - cp)dcp/dt, dx3/dt = Rsin(2n/3 - cp)dcp/dt.
Условием возникновения свободных гармонических колебаний является постоянство полной, в рассматриваемом случае кинетической энергии системы T = 0,5mR2[sin2cp + sin2(n/3 - cp) + sin2(2n/3 - cp)](dcp/dt)2 = 0,75mR2(dcp/dt)2 = const,
dcp/dt = const.
С учетом (2)-(4) (7) принимает вид
xi = Rcos(ro0t + p0), x2 = Rcos(n/3 - ra0t - p0), x3 = Rcos(2n/3 - ra0t - p0). (8)
Пусть x1(0) = x10, (dx1/dt)(0) = v10. Тогда
x I x2-
cosp0 = x10/R, cpn = arccos—10 = arcsin.. Д —20 ; (9)
0 R \ R
vln
-Rro0sin(ro00 + cp0) = v10, ю0 =—. . (10)
VR2 - xi20
Таким образом, все три груза совершают свободные гармонические колебания, обмениваясь между собою кинетической энергией.
Упругая kk-система
Система отличается от дада-системы тем, что вместо массивных грузов к вертикальному и горизонтальному стержням прикреплены пружины с коэффициентом упругости к, а середина промежуточного стержня связана с началом координат посредством кривошипа. При ф = 0 первая пружина максимально сжата, вторая - не деформирована. При ф = п/2 первая пружина не деформирована, вторая -максимально сжата. При ф = п первая пружина максимально растянута, вторая -не деформирована. Поскольку все действующие на систему силы потенциальные, то уравнение Лагранжа второго рода может быть записано в виде:
Таким образом, уравнение Лагранжа для системы является тождеством. Суммарный момент М = іїШйф = 0 при любой скорости вращения кривошипа. При начальных условиях (3) и в соответствии с аналогом первого закона Ньютона для вращательного движения
При этом (1) приводится к (5) и (6).
Таким образом, пружины с коэффициентом упругости к совершают свободные гармонические колебания.
В рассмотренной колебательной системе происходит взаимный обмен потенциальной энергией между упругими элементами. При ср = 0 потенциальная энергия первой пружины максимальна, а второй - равна нулю. После этого вторая пружина начинает сжиматься за счет энергии первой пружины, которая начинает разжиматься.
Система отличается от дадада-системы тем, что грузы заменены пружинами с коэффициентом упругости к. Кроме того, она оснащена кривошипом, связывающим точку пересечения координатных осей с центром треугольника x2x2x3. При этом центр треугольника жестко связан с его сторонами. С учетом (7) потенциальная (она же полная) энергия системы
П = 0,5kR2[cos2cp + cos2(n/3 - cp) + cos2(2n/3 - cp)] = 0,75kR2 = const.
Суммарный момент M = dn/dcp = 0 при любой скорости вращения кривошипа. При начальных условиях (3) и в соответствии с аналогом первого закона Ньютона для вращательного движения выполняется (11). При этом (7) приводится к (8)-
Таким образом, все три пружины совершают свободные гармонические колебания, обмениваясь между собою потенциальной энергией.
Потенциальная энергия
П = 0,5kxr2 + 0,5kx2 = 0,5kl2 (cos2 ф + sin2 ф) = 0,5kl2.
(11)
Упругая kkk-система
(10).
Обобщение на и-мерный случай
Пусть n координатных осей Ox1, Ox2, ... Oxn, лежат в одной плоскости Z и пересекаются в одной точке. Соседние оси ориентированы под углами п/n относительно друг друга. Точка пересечения O является началом произвольно направленного вектора R, принадлежащего Z.
Теорема 3. Координаты x1, x2, ... xn проекций конца вектора R на оси Ox2, Ox2, ..., Oxn являются вершинами правильного n-угольника, размер которого не зависит от направления R, а центр совпадает с серединой R.
Доказательство. Координаты проекций
x1 = Rcoscp, x2 = Rcos(n/n - cp ), x3 = Rcos(2n/n - p),
xi = Rcos[(i - 1)n/n - cp], xn = Rcos[(n - 1)n/n - p]. (12)
По теореме косинусов (xxi+1)2 = R2{cos2[(i-1)n/n-p] + cos2(in/n-p) - 2cos[(i-1)n/n-p]cos(in/n-p )cos(n/n)} = = R2{cos2(in/n - p - n/n) + (cos(in/n)cosp + sin(in/n)sinp)2 -- 2[cos(n/n)(cos(in/n)cosp + sin(in/n)sinp) + sin(n/n)(sin(in/n)cosp -- cos(in/n)sinp)]x(cos(in/n)cosp + sin(in/n)sinp)cos(n/n)} =
ООО OOOOO 0
= R {cos (n/n)cos (in/n)cos cp + cos (n/n)sin (in/n)sin p + sin (n/n)sin2(in/n)cos p +
+ sin2(n/n)cos2(in/n)sin2p + 2cos2(n/n)cos(in/n)cospsin(in/n)sinp +
+ 2cos(n/n)cos(in/n)cos2p sin(n/n)sin(in/n) - 2cos(n/n)cos2(in/n)cosp sin(n/n)sinp +
+ 2cos(n/n)sin2(in/n)sinpsin(n/n)cosp - 2cos(n/n)sin(in/n)sin2psin(n/n)cos(in/n) -- 2sin2(n/n)sin(in/n)cospcos(in/n)sinp + cos2(in/n)cos2p + sin2(in/n)sin2p +
+ 2cos(in/n)cospsin(in/n)sinp - 2cos2(n/n)cos2(in/n)cos2p -
- 2cos2(n/n)sin(in/n)sinp cos(in/n)cosp - 2sin(n/n)sin(in/n)cos2pcos(in/n)cos(n/n) +
+ 2sin(n/n)cos2(in/n)sinpcospcos(n/n) - 2cos2(n/n)cos(in/n)cospsin(in/n)sinp -
- 2cos2(n/n) sin2(in/n)sin2p - 2sin(n/n)sin2(in/n)cospsinpcos(n/n) +
+ 2sin(n/n)cos(in/n)sin2psin(in/n)cos(n/n)} = R2{cos2p(sin2(n/n)sin2(in/n) -
00 О ООО 00
- cos (n/n)cos (in/n) + cos (in/n)) + sin p(sin (n/n)cos (in/n) - cos (n/n)sin (in/n) +
0 0 0 0 0 о 0
+ sin (in/n))}= R {cos p(sin (n/n)sin (in/n) + cos (in/n)sin (n/n)) +
ООО о 0
+ sin p(sin (n/n)cos (in/n) + sin (in/n)sin (n/n))} =
ООО 0 0 0 0
= R {cos p sin (n/n) + sin psin (n/n)} = R sin (n/n).
Таким образом, все стороны многоугольника равны и не зависят от p. Пусть г
- середина вектора R. По теореме косинусов (xr)2 = R2{cos2[(i-1)n/n - p] + 1/4 -2cos[(i-1)n/n - p](1/2)cos[(i-1)n/n - p]} = R2/4. Все лучи, соединяющие г с вершинами многоугольника, равны между собой. Таким образом, многоугольник правильный. Теорема доказана.
Для инертной системы скорости грузов с учетом (12)
dx1/dt = -Rsinpdp/dt, dx2/dt = Rsin(n/n - p )dp/dt, dxjdt = Rsin[(i - 1)n/n - p]dp/dt, dxn/dt = Rsin[(n - 1)n/n - p]dp/dt. Условием возникновения свободных гармонических колебаний является постоянство полной, в рассматриваемом случае кинетической энергии системы T = 0,5mR2{sin2p + sin2(n/n - p) + ... + sin2[(i - 1)n/n - p] + ...
+ sin2[(n - 1)n/n - p]}(dp/dt)2 (13)
Теорема 4. При n > 2 справедлива формула
" (І - 1)п
X Sin2
i=1
-±ф
Доказательство.
n
X Sin2
І=1
(І - 1)п
n
-±ф
V /1 1
= Х (2 - ^cos
n
2.
(І - 1)2п
n
±2ф
n 1 „ n (i - 1)2п 1 . „ n . (І - 1)2п
-----cos2ф> cos--------------+ — Sin2ф> sin-------------
i=i n 2 i=i
2 2
n
Пусть п единичных векторов образуют центрально симметричную звезду таким образом, чтобы один из векторов совпал с осью абсцисс. В силу симметрии сумма векторов равна нулю. Сумма проекций векторов на любую ось тоже равна нулю. Таким образом, две последние суммы равны нулю. Теорема доказана. Аналогично доказывается Теорема 5. При п > 2 справедлива формула
' (і - 1)п
Xcos2
i=1
-±ф
n
2
Кинетическая энергия инертной системы (13)
T = 0,25nmR2(dp/dt)2 = const, dp/dt = const С учетом (2) - (4) (12) принимает вид
xr = Rcos(ra0t + ф0), x2 = Rcos(n/n -xt = Rcos[(i - 1)n/n - ro0t - ф0], При этом в соответствии с (9), (10)
„ = Vx10
ro0t - ф0), x3 = Rcos(2n/n - ro0t - ф0), xn = Rcos[(n - 1)n/n - ro0t - ф0]. (14)
Vr2 - x2
10
ф0 = arccos—^. 0R
(15)
Таким образом, все n грузов совершают свободные гармонические колебания, обмениваясь между собою кинетической энергией.
Многоугольник x1x2.x„ совершает двойное вращение - вокруг своего центра и вокруг точки O.
В отличие от инертной для упругой системы условием возникновения свободных гармонических колебаний являются постоянство потенциальной (полной) энергии и наличие кривошипа, связывающего точку пересечения координатных осей с центром многоугольника x1x2.x„, который в свою очередь жестко связан с его сторонами.
ООО о о
П = 0,5kR {cos p + cos (п/n - p) + ... + cos [(i-1)n/n - p] + ... +cos [(n-1)n/n - p]}.
В соответствии с теоремой 5
П = 0,25nkR2 = const.
Суммарный момент M = dn/dp = 0 при любой скорости вращения кривошипа. При начальных условиях (3), в соответствии с аналогом первого закона Ньютона для вращательного движения и (11), выполняются (14), (15).
Таким образом, все n пружин совершают свободные гармонические колебания, обмениваясь между собою потенциальной энергией.
Заключение
Установлена возможность возникновения свободных гармонических колебаний в системах, состоящих только из инертных (mm-, mmm-, даи-системах) или только упругих элементов (kk-, kkk-, ^-системах), которая реализуется при обеспечении сдвига по фазе между колебаниями элементов.
В отличие от традиционных [2-4] или смешанных [1] колебательных систем при энергообмене между однородными элементами представленных систем вид энергии не меняется. В mm-, mmm-, mи-системах происходит взаимный обмен кинетической энергией между инертными элементами. В kk-, kkk-, ^-системах - потенциальной энергией между упругими элементами. При этом суммарная энергия систем при колебаниях не изменяется.
Другим отличием является то, что частоты свободных колебаний систем с однородными элементами не зависят от параметров элементов и определяются исключительно начальными условиями. Иначе говоря, рассмотренные системы могут совершать свободные гармонические колебания с любой изначально заданной частотой.
Колебательные свойства mm-, mmm-, mи-систем могут учитываться при проектировании двигателей внутреннего сгорания, поршневых пневмосистем и прочих преобразователей возвратно-поступательного движения во вращательное в плане взаимной компенсации реактивного характера масс движущихся частей - поршней, штоков и пр.
Принцип обмена энергией между элементами в kk-, kkk-, ^-системах может использоваться для самонейтрализации квазиупругих воздействий в пневмосистемах.
ЛИТЕРАТУРА
1. Попов И.П. Реактивные элементы электрических цепей с «неэлектрическими» параметрами // Вестн. Самарского государственного технического университета. Технические науки. 2010. № 4(27). С. 166-173.
2. Баркгаузен Г. Введение в учение о колебаниях. М.: Госэнергоиздат, 1934. 116 с.
3. Tongue B. Principles of Vibration. Oxford University Press, 2001. 367 p.
4. Thompson W.T. Theory of Vibrations. Nelson Thornes Ltd., 1996. 295 p.
5. Попов И.П. Свободные гармонические колебания в системах с однородными элементами // ПММ. 2012. Т. 76. Вып. 4. С. 546-549.
Статья поступила 20.03.2012 г.
Popov I.P. VIBRATION SYSTEMS CONSISTING ONLY OF INERT OR ONLY ELASTIC ELEMENTS AND THE APPEARANCE OF FREE HARMONIC VIBRATIONS IN THEM. We consider mechanical vibration systems consisting only of inert (mm-, mmm-, mn-systems) or only elastic (kk-, kkk-, k”-systems) components. The possibility of the of free harmonic vibrations appearance in such systems is shown. In the mm-, mmm-, mn-systems, the mutual exchange of kinetic energy between the inert elements occurs; in kk-, kkk-, k”-systems, the mutual exchange of the potential energy between elastic elements.
Keywords: vibration, inert, elastic, harmonic, frequency.
POPOV Igor Pavlovich (Kurgan State University)
E-mail: [email protected], [email protected]