CC BY
43
ОЛЬШАНСКИИ Василий Павлович
(Харьковский национальный технический университет сельского хозяйства)
OL'SHANSKIY Vasiliy Pavlovich
(Kharkiv, Ukraine, Kharkiv Petro Vasylenko National Technical University of Agriculture)
ОЛЬШАНСКИЙ Станислав Васильевич
(Национальный технический университет «Харьковский политехнический институт»)
OL'SHANSKIY Stanislav Vasil'evich
(Kharkiv, Ukraine, National Technical University «Kharkiv Polytechnic Institute»)
УДК 534.1
О прохождении резонанса в механической системе переменной массы
В.П. Ольшанский, С.В. Ольшанский
При исследовании прохождения резонанса в механической системе переменной массы машин, работающих в зарезонансных режимах, основной задачей является определение коэффициента динамичности в осцилляторе линейно-переменной массы при его нестационарных резонансных колебаниях под действием внешней силы переменной частоты. В отличие от существующих публикаций впервые получено аналитическое решение задачи Коши в квадратурах для случая, когда частота возмущающей силы меняется по линейному закону, как и масса. В статье с помощью ВКБ-метода приближенно описаны нестационарные резонансные колебания осциллятора линейно-переменной массы под действием внешней силы, частота которой меняется по линейному закону. Рассмотрены случаи разгона и выбега системы при увеличении и уменьшении массы осциллятора. Расчеты показали, что и в осцилляторе переменной массы при переходе через резонанс максимальные амплитуды колебаний наблюдаются не в момент равенства меняющихся во времени собственной частоты и частоты внешнего воздействия, а позже, т. е. при некотором смещении от этого момента времени. Величина смещения (запаздывания) уменьшается при увеличении силы вязкого трения. Полученные формулы позволяют рассчитать амплитуды опасных резонансных колебаний.
Ключевые слова: осциллятор, система переменной массы, коэффициент динамичности, резонанс, ВКБ-метод.
On the passage through a resonance in a variable-mass mechanical system
V.P. Ol'shanskiy, S.V. Ol'shanskiy
The study ofpassage through resonance is an essential problem in many technical applications especially when analyzing the dynamics of above-resonance machines. In this case, it is important to determine the dynamic factor of transient resonant oscillations of a linear variable-mass oscillator under the action of an external force of variable frequency. In contrast to the existing literature, a new analytical solution to the Cauchy problem is found by quadratures for case when the disturbing force frequency and mass vary linearly. In this paper, the WKB approximation is used to describe transient resonance oscillations of a linearly varying mass under the action of an external force whose frequency varies linearly as well. The acceleration and coasting of the system with an increasing and decreasing oscillating mass are considered. The calculations showed that the resonance mode of a variable mass oscillator is not observed when the time-varying natural frequency coincides with the frequency of an external action but at a later time. The offset time decreases with increasing viscous forces. The obtained formulas make it possible to calculate the amplitudes ofdangerous resonant oscillations.
Keywords: oscillator, variable-mass system, dynamic factor, resonance, WKB method.
Т/Тсследованием резонанса системы с одной степенью свободы при действии возмущающей силы переменной частоты, а также переменной во времени амплитуды, занимались многие авторы [1—3]. Менее исследованным является вопрос прохождения резонанса осциллятором с переменной жестокостью или массой. Проблема такого резонанса возникает в устройствах, в которых систематически используют механизмы с переменной массой или жесткостью звеньев [4]: механизмы экскаваторов и грейферов; машины шахтного подъема; вибрационные конвейеры; емкости и бункеры для загрузки жидкостью или сыпучим веществом; подрессоренные платформы транспортных средств, загружаемые сыпучим грузом и пр. В отличие от публикации А.П. Бессонова в обзоре [5], где моделировалось прохождение резонанса на аналоговой машине, в данной работе строится приближенное аналитическое решение задачи Коши в квадратурах для случая, когда частота возмущающей силы меняется по линейному закону, как и масса.
Цель работы — определение коэффициента динамичности в осцилляторе линейно-переменной массы при его нестационарных резонансных колебаниях под действием внешней силы переменной частоты.
Решение задачи в квадратурах. С учетом дис-сипативной и реактивной сил, вынужденные резонансные колебания осциллятора описываются дифференциальным уравнением
d2 х m 0 ys + ц dx
dt2
- +
m
,(1+ Yt) dt
+ -
ш
Q
m 0 (1+ Yt)
sin
1 + Yt
,2 \
x
at +
btl 2
(1)
где т0 — начальная масса осциллятора; т0 у — скорость изменения массы во времени t; ш2 = с / т0; с — коэффициент жесткости пружины; ц — коэффициент вязкого трения; в — коэффициент реактивности переменной массы, 0 <в <1; 0 — амплитуда возмущающей силы; а, Ь — постоянные коэффициенты.
Решим уравнение (1) при нулевых начальных условиях
„ dx
х (0) = 0; — dt
= 0.
(2)
t=0
В соответствии с методом Лагранжа, решение задачи Коши, представленной выражениями (1) и (2), имеет вид
х (
(t ) = C1 (t )х 1 (t) + с 2 (t )х 2 (t).
Здесь
Q i-х21 (t)sin(at + bt2 /2)
^ 2 (t) = +—J —-ттт:-ч-dt;
m
1 х 2
00
A(t )(1+ Yt)
(3)
(4)
A(t ) = х 1 (t) % - х 2 (t) ^;
dt dt
х 1 (t) , х2 (t) удовлетворяют однородному урав-
нению
d2х m0ys + ц dx ш2 dt2 m0 (1 + Yt) dt 1+ Yt
-х = 0.
(5)
При ц = в = 0 в (5) функции х 1 (t) и х2 (0 выражаются через функции Бесселя индекса единица [6]. Далее, для упрощения задачи используем их ВКБ-приближения (приближения Вентцеля—Крамерса—Бриллюэна) [7]. Чтобы получить эти приближения подставим в (5) произведение
х (t) = y(t)
exp
-1J
2
m0 ys + ц
m
(1+ Yt)
dt
= y(t )(1+ Yt )-2( Е+ц/ m0Y).
Для определения функции y(t) используем уравнение
d2 y dt2
+
ш
1 + Yt
m 0 Ys + ц
m
(1+ Yt)
+
+
1 Y(m0 Ys + Ц)
2
m
(1+ Yt)2
y=0,
которое решаем приближенно ВКБ-методом [8], в предположении, что ш/|y|>>1 Пренеб-
регая слагаемыми, порядок малости которых выше (ш / |y|) , получаем
2
lit):
= exp
±i®f
dt
+1fdt
Vl+ Yt 4J 1+ yt,
/ = V-1
Вычислив эти интегралы, с учетом записанного выше выражения x(г) и известной формулы Эйлера, находим ВКБ-приближения фундаментальных решений уравнения (5):
х,
где
1 (t) = (п / По У cos(n); х2 (t) = (п / По =)ß sin(n),
Ц
2® г,- 1
по =—; п=nw1+Yt; ß=ö-£
y 2 ym
Подставляя (6) в (4), с учетом, что
/ \2P-1
п
f.Jill. д(() = m
dn УПо 2 VПо/
решение (3) сводим к виду
; 1+Yt =
п2
öl I
х(t) =--1 [01 (^cos^ + a2 (п) sin^J (7)
Здесь
п/по
h (п) = -по f 7- ¡^п(по 7)х
Xsin
1 (72 -1)la + b(72 -1)
2Y
d7;
Xsin
i (п) = по f 7 в cos(по 7)x
1
1 ( 72-1)fa + b ( 72-1)
2Y
dy. (8)
Используя формулы (7) и (8) для вычисления х(г), находим выражение
х (t ) = Q f(l±*
то® о U + Yv
sin[nо(л/1+ Yt -л/1+ YT)]
X-, -sin
X
1
ат + 2 Ьт 2jdx. (9)
^ ут
В пределе, когда у ^ о, ß ^ — а>, а
lim по W1+ Yt -л/1+ ут ) = ®(t - т),
lim
у^0
/л, лв/2 1 + Yt \
1+ y^
= exp
Ц
2то
(t - т)
выражение (9) сводится к более компактной формуле:
х (t ) = Q f
' Ч t-т )
sin (№
то® 'о
(t - т)
sin
1 2 \ ат + — Ьт Idт.
2 I
Это приближенное выражение, полученное предельным переходом для осциллятора постоянной массы m = m0 = const, совпадает с известным точным решением [2], после замены в нем
®2 -
\2
Ц
то I
. Следовательно, по-
та на ®1 = ^
грешность принятого ВКБ-приближения мала
с 1 при — >> -
то 4
\2 Ц\
то I
Выделяя в (9) амплитуду вынужденных колебаний атх(г), получаем формулу максимального коэффициента динамичности перемещений осциллятора
K д =
сат х1
(t)
Q
= ю^Ь12 (t) + Ь22 (t), (1Q)
где
Ь1 (t ) = f И Гх
о 11+ Yt71 sin[nо (л/1+ Yt -V1+ Y7)]
X
Xcos
л/1+ Yt
X
7
it ( 7-1) + ^ ( 72-1)
d7; (11)
Ь2 (t) = f \/2 x
2 i U + Yt71
sin[^ (V 1+ Yt -V1+ Yt7)]
X
Xsin
>/1+ Yt
7
it (7-1) + ^ ( 72-1)
X
d7.
Интегралы (11) не выражаются через элементарные или затабулированные специаль-
т
e
о
a
ные функции, поэтому их вычисляют на компьютере.
Численные результаты и их анализ. Рассмотрим колебания осциллятора возрастающей массы, для следующих исходных данных: т0 = 10 кг; с = 160 Н/м; в = 0,5; у = 0,025 с1; 0 = 1Н. Зависимости коэффициента динамичности при разгоне осциллятора для а = 0,266 с-1, Ь = 0,15 с-2 представлены на рис. 1, а; для а = 6,266 с-1, Ь = — 0,15 с-2 при выбеге осциллятора — на рис. 1, б.
Для принятых исходных данных совпадение собственной частоты осциллятора возрастающей массы с частотой внешнего возмущения происходит при t = 20 с. Однако, как при разгоне, так и при выбеге, максимум амплитуды колебаний достигается при t > 20 с. На рисунке видно смещение нестационарного резонанса вправо, т. е. запаздывание его во времени. При выбеге уровень амплитуд колебаний осциллятора возрастающей массы немного выше, чем при разгоне.
Сравним результаты, полученные приближенно с помощью интегралов (11) с результатами численного интегрирования уравнения (1). Зависимость вынужденных колебаний в резонансной области, к которым приводит численное решение задачи Коши для выражений (1) и (2) на компьютере, при разгоне и выбеге представлены на рис. 2. В расчетах использованы те же исходные (численные) данные, что и при получении кривых 2 на рис. 1.
Сравнивая амплитудные значения перемещений при разгоне (см. рис. 1, а и 2, а), а также
Рис. 1. Зависимости коэффциента динамичности Кд от времени t при разгоне (а) и при выбеге (б):
1 — ц = 3 кг/с; 2 — ц = 6 кг/с; 3 - ц = 9 кг/с
при выбеге (см. рис. 1, б и 2, б), можно отметить близкие значения коэффициентов динамичности. Следовательно, приближенные формулы (11) обладают достаточной точностью.
Рис. 2. Зависимость вынужденных колебаний х от времени t при разгоне (а) и при выбеге (б)
Рассмотрим колебания осциллятора убывающей массы для следующих исходных данных: т0 = 10 кг; с = 160 Н/м; в = 0,5; у = —0,05 с-1; 0 =1 Н. Зависимости коэффициента динамичности, соответствующие разгону осциллятора при а = 1,6569 с-1 и Ь = 0,4 с-2, представлены на рис. 3, а, а при выбеге осциллятора при а = 9,6569 с-1 и Ь = —0,4 с-2 — на рис. 3, б.
В данном примере совпадение собственной частоты осциллятора убывающей массы с частотой внешней силы происходит при t = 10 с. Но, как при разгоне, так и при выбеге, максимальные резонансные амплитуды достигаются при t > 10 с. На кривых нестационарный резонанс смещается вправо. Уровень амплитуд
Рис. 3. Зависимости коэффциента динамичности Кд от времени t при разгоне (а) и при выбеге (б):
1 — ц = 3 кг/с; 2 — ц = 5 кг/с; 3 — ц = 7 кг/с
колебаний при разгоне осциллятора убывающей массы оказался выше уровня амплитуд при выбеге.
Выводы
1. Применение ВКБ-метода позволяет построить в квадратурах приближенное решение задачи о нестационарном резонансе в осцилляторе переменной массы при линейном изменении частоты внешней силы.
2. Расчеты показали, что и в осцилляторе переменной массы при переходе через резонанс максимальные амплитуды колебаний наблюдаются не в момент равенства меняющихся во времени собственной частоты и частоты внешнего воздействия, а при некотором запаздыванием от этого момента времени.
3. Установлено, что величина смещения (запаздывания) уменьшается при увеличении силы вязкого трения. Сила трения также оказывает значительное влияние на уровень резонансных амплитуд колебаний.
Литература
[1] Katsuhiko Ogata. System Dynamics. University of Minnesota, 2005. 617 p.
[2] Ольшанский В.П., Ольшанский С.В. Метод ВКБ в расчетах нестационарных колебаний осцилляторов. Харьков, Мюькдрук, 2014. 264 с.
[3] Митропольский Ю.А. Избранные труды. В 2 т. Т. 2. Киев, Наукова думка, 2012. 504 с.
[4] Cveticanin L. Dynamics of Machines with Variable Mass. Taylor & Francis Ltd, 1998. 300 p.
[5] Cveticanin L. A review on dynamics of mass variable systems. Journal of the Serbian Society for Computational Mechanics, 2012, vol. 6, no. 1, pp. 56—74.
[6] Ольшанский В.П., Ольшанский С.В. Моделирование колебаний осциллятора линейно-переменной массы при импульсном нагружении. ВСник НТУ «ХП1»: Матема-тичне моделювання в технiцi та технологiях, 2013, № 37(1010), с. 125—130.
[7] Karnakov B.M., Krainov V.P. WKB Approximation in Atomic Physics. Springer, 2013. 176 p.
[8] Зайцев В.Ф., Полянин А.Д. Справочник по обыкновенным дифференциальным уравнениям. Москва, ФИЗМАТЛИТ, 2001. 576 с.
References
[1] Katsuhiko Ogata. System Dynamics. Fourth Edition. University of Minnesota, 2005. 617 p.
[2] Ol'shanskii V.P., Ol'shanskii S.V. Metod VKB raschetakh nestatsionarnykh kolebanii ostsilliatorov [WKB method for unsteady oscillator]. Khar'kov, Mis'kdruk publ., 2014. 264 p.
[3] Mitropol'skii Iu.A. Izbrannye trudy [Selected Works]. Kiev, Naukova dumka publ., vol. 2, 2012. 504 p.
[4] Cveticanin L. Dynamics of Machines with Variable Mass. Taylor & Francis Ltd, 1998. 300 p.
[5] Cveticanin L. A review on dynamics of mass variable systems. Journal of the Serbian Society for Computational Mechanics, 2012, vol. 6, no. 1, pp. 56—74.
[6] Ol'shanskii V.P., Ol'shanskii S.V. Modelirovanie kolebanii ostsilliatora lineino-peremennoi massy pri impul'snom nagruzhenii [Modeling oscillator linearly variable mass under impact loading]. VisnikNTU«KhPI»: Matematichne modeliuvannia v tekhnitsi ta tekhnologiiakh [Herald of NTU «KPI»: Mathematical modeling in engineering and technologies]. 2013, no. 37(1010), pp. 125—130.
[7] Karnakov B.M., Krainov V.P. WKB Approximation in Atomic Physics. Springer, 2013. 176 p.
[8] ZaitsevV.F, PolianinA.D. Spravochnik po obyknovennym differentsial'nym uravneniiam [Handbook of Differential Equations]. Moscow, FIZMATLIT publ., 2001. 576 p.
Статья поступила в редакцию 30.01.2014
Информация об авторах
ОЛЬШАНСКИЙ Василий Павлович (Харьков) — доктор физико-математических наук, профессор кафедры «Теоретической механики и деталей машин». Харьковский национальный технический университет сельского хозяйства (61002, Харьков, Украина, ул. Артема, 44).
ОЛЬШАНСКИЙ Станислав Васильевич (Харьков) — кандидат физико-математических наук, доцент кафедры «Системы и процессы управления». Национальный технический университет «Харьковский политехнический институт» (61002, Харьков, Украина, ул. Фрунзе, 21, e-mail: stasolsh77@gmail.com).
Information about the authors
OL'SHANSKIY Vasiliy Pavlovich (Kharkiv) — Dr. Sc. (Phys. Math.), Professor of «Theoretical Mechanics and Machine Parts» Department. Kharkiv Petro Vasylenko National Technical University of Agriculture (KhPVNTUA, Artema str., 44, 61002, Kharkiv, Ukraine).
OL'SHANSKIY Stanislav Vasil'evich (Kharkiv) — Cand. Sc. (Phys. Math.), Associate Professor of «Control Processes and Systems» Department. National Technical University «Kharkiv Polytechnic Institute» (NTU KhPI, Frunze str., 21, 61002, Kharkiv, Ukraine, e-mail: stasolsh77@gmail.com).
