2024 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 11 (69). Вып. 1
МАТЕМАТИКА. МЕХАНИКА. АСТРОНОМИЯ
МЕХАНИКА
УДК 533:534.13 МБС 76Б99
Затухающие колебания цилиндра с соосным диском и со стабилизатором
Д. В. Кауфман, А. Н. Рябинин
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Кауфман Д. В., Рябинин А. Н. Затухающие колебания цилиндра с соосным диском и со стабилизатором // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 1. С. 161-170. https://doi.org/10.21638/spbu01.2024.111
Изучаются затухающие вращательные колебания цилиндра, который в головной части снабжен соосным диском, а в хвостовой части имеет стабилизатор. Удлинение цилиндра (отношение длины к диаметру) равно девяти. Цилиндр крепится в рабочей части аэродинамической трубы малых скоростей на проволочной подвеске, содержащей стальные пружины. В положении равновесия ось цилиндра горизонтальна и параллельна вектору скорости набегающего потока. К одной из пружин подвески присоединен полупроводниковый тензопреобразователь, измеряющий во время колебаний зависимость натяжения пружин от времени. Напряжение на выходе тензопреобразо-вателя поступает на РС-осциллограф. Цифровой сигнал осциллографа передается на компьютер. После калибровки прибора определялась частота и амплитуда затухающих вращательных колебаний вокруг горизонтальной оси, проходящей через центр цилиндра и перпендикулярной вектору скорости набегающего потока. Под действием воздушного потока увеличивается скорость затухания вращательных колебаний цилиндра. Влияние воздушного потока описывается аналогами вращательных производных, которые в случае плохо обтекаемых тел зависят от амплитуды колебаний угла наклона тела и от амплитуды угловой скорости. Предложена простая модель влияния стабилизатора на вращательные производные.
Ключевые слова: аэродинамическая труба, вращательные колебания, тензометрия, цилиндр.
© Санкт-Петербургский государственный университет, 2024 https://doi.org/10.21638/spbu01.2024.111
1. Введение. В работе [1] изучались вращательные затухающие колебания цилиндра, в головной части которого на тонкой ножке находится коаксиально закрепленный диск. При переноске груза под вертолетом установка такого диска обеспечивает значительное снижение лобового сопротивления [2, 3]. Для ориентации груза в направлении, обеспечивающем минимальное лобовое сопротивление, используется стабилизатор, расположенный в хвостовой части груза [3, 4]. В настоящей работе экспериментально исследуется затухание вращательных колебаний цилиндра, имеющего в передней части коаксиальный диск, а в задней — стабилизатор. Отношение длины цилиндра Ь к его диаметру Б равно 9. Стабилизатор состоял из двух прямоугольных трапеций. Высота каждой трапеции равна Б/2. Площадь стабилизатора в два раза больше площади поперечного сечения цилиндра. Схема расположения стабилизатора приведена на рис. 1, б.
2. Экспериментальный метод. Эксперименты выполнялись на аэродинамической трубе малых скоростей АТ-12 Санкт-Петербургского университета [5]. Труба имеет открытую рабочую часть и сопло круглого сечения, выходной диаметр сопла 1.5 м. Цилиндр подвешивался в рабочей части аэродинамической трубы на проволочной подвеске, содержащей две пружины (рис. 1). В состоянии равновесия цилиндр располагался горизонтально. Его ось была параллельна вектору средней скорости набегающего потока. Подвеска позволяла цилиндру поворачиваться вокруг горизонтальной оси, перпендикулярной вектору скорости набегающего потока. Упругие силы пружин подвески при этом стремились возвратить цилиндр в исходное состояние равновесия.
Рис. 1. Схема эксперимента (а) и расположение стабилизатора на цилиндре (б): 1 — цилиндр; 2 — диск; 3 — стабилизатор; 4 — контргруз; 5 — пружины; 6 — полупроводниковый тензопреобразователь; 7 — РС-осциллограф; 8 — компьютер.
Нижняя пружина подвески прикреплена к полупроводниковому тензопреоб-разователю С-50, измеряющему натяжение пружины. РС-осциллограф Velleman-PCS500A преобразовывал аналоговый сигнал тензопреобразователя в цифровой и передавал его в компьютер. Во время экспериментов цилиндр отклоняли от положения равновесия на некоторый угол, затем отпускали. Начинались затухающие вращательные колебания. Натяжение пружины считывалось в течение 17 с при частоте 100 Гц. В эксперименте менялся диаметр с! соосного диска, принимая три значения:
0.625В, 0.75В и 0.875В, расстояние д от диска до переднего торца цилиндра также принимало три значения: 0.5В, 0.65В и 0.8В.
Были проведены калибровочные эксперименты, которые позволили связать амплитуду натяжения пружины с амплитудой вращательных колебаний.
Все эксперименты повторялись не менее 5 раз. Это позволило на основе разброса экспериментальных данных оценить случайную погрешность измеряемых величин. Методика обработки экспериментальных данных, основанная на методе наименьших квадратов, изложена в статье [1].
3. Описание затуханий колебаний цилиндра. Скорость затухания колебаний тел в воздушном потоке принято характеризовать вращательными производными [6]. Для удобообтекаемых тел коэффициент аэродинамического момента можно считать линейной функцией малых углов наклона и их производных:
. ачт2
1хв + гв + кв = (т% + твх)--в + т^Ьв, (1)
V
где I, — компонента момента инерции цилиндра, соответствующая оси вращения; в — угол наклона тела; г — коэффициент трения в подвеске; к — приведенная упругость пружин; а = (м2/2 — скоростной напор; р — плотность воздуха; в = пВ2/4 —
< ) й й
характерная площадь тела; V — скорость воздуха; т", таг, таг — вращательные производные коэффициента момента; ш — угловая скорость. В правой части уравнения (1) находятся члены, соответствующие аэродинамическому моменту. С использованием обозначений
аТ
ии2 = к/О, = ¡1 = (то" + тв)—---г/1г
VI,.
уравнение (1) преобразуется к виду
в + ю2в = цв + Пв. (2)
Однако для плохо обтекаемых тел отклонения от линейности могут быть существенными. На рис. 2 представлена зависимость коэффициента момента т, от угла наклона цилиндра в, имеющего удлинение Ь/В = 9. Момент измерен аэродинамическими весами с проволочной подвеской. Цилиндр во время измерения неподвижен.
На графике видно отклонение зависимости от линейности в области малых углов наклона. Можно ожидать, что отклонения от линейности присутствуют и в зависимости коэффициента момента от производной угла наклона по времени.
В общем случае, когда зависимость момента от в и в нелинейная, вместо уравнения (1) нужно использовать уравнение движения:
1гв + г в + кв = авЬ! (в, в). (3)
Введя обозначения
п = авЬ/1,, = г/1,, получим уравнение движения в виде
в + т2в = п1 (в, в) - П1 в. (4)
Рис. 2. Зависимость коэффициента момента от угла наклона неподвижного цилиндра.
Считая малыми параметры п и ni, уравнение (4) можно решать приближенным методом Крылова — Боголюбова [7]. Введем новые переменные: амплитуду A и ф фазу колебаний,
в = A cos ф, ф = wt +
Амплитуда A и сдвиг фазы являются медленно меняющимися переменными. В методе Крылова — Боголюбова в первом приближении для производных A и ф по времени используются выражения
г2п
г2п
А=--— I /(A cos ф7 — wAsin ф) sin фс1ф — —— / sin2 фс1ф,
2nw J0 2п J0
ф = w —
п
г2п
2Anw
f (A cos ф, —wA sin ф) cos фdф —
2тг
r2n
sin ф cos фdф.
Взяв интегралы во вторых слагаемых, получим
A = —
n f2n niA
f(Acos ф, — wAsin ф) sin фс1ф--,
ф = w —
2nw J о п
2
r2n
2Anw
f (A cos ф, —wA sin ф) cos фdф.
(5)
(6)
Применим процедуру Крылова — Боголюбова к уравнению (2), считая ц и Q
малыми параметрами:
A = A
2'
ф = w —
П
2ли
(7)
Приравняв правые части уравнения (5) и первого уравнения (7), правые части уравнения (6) и второго уравнения (7), получим
г2п
m. + mr = —
AwnL
f (A cos ф, —wA sin ф) sin фdф,
0
0
0
0
V
1 г2п
то® = —— / /(Асовф,—^лАвтф) совфг1ф. Ап .] о
>о
Отклонения от линейности могут приводить к зависимостям параметров т" + тйг и т, от амплитуд колебаний угла наклона А и производной угла наклона Аи. Однако эксперименты показали, что в некоторых диапазонах амплитуд колебаний зависимость логарифма амплитуды колебаний угла наклона от времени близка к линейной, и зависимость квадрата частоты колебаний линейно зависит от скоростного напора. Поэтому затухание можно описывать аналогами вращательных производных, которые, в точном смысле термина, вращательными производными не являются. При этом нужно указывать, в каком диапазоне амплитуд колебаний углов наклона и угловых скоростей осуществляется затухание колебаний.
В соответствии с алгоритмом определения вращательных производных, изложенном в статье [1], по линейной зависимости коэффициента затухания колебаний от скорости воздушного потока п = a+bv сумма вращательных производных т"+т, выражается через коэффициент Ь, который находится по экспериментальным данным методом наименьших квадратов:
в 4 Ых 4Ък
то_ + то_ — —
рвЬ2 рви!2 Ь2'
где частота колебаний и определяется в экспериментах с затухающими колебаниями с выключенным вентилятором аэродинамической трубы. Частота колебаний ф линейно зависит от скоростного напора а:
ф = и — т,авЬ/(2и1,) = и — т,авЬи/(2к).
Определив наклон экспериментальной зависимости ф = — Ь\а, можно вычислить вращательную производную тйг. В экспериментах коэффициент лишь приближенно совпадал с частотой колебаний, измеренной при затухающих колебаниях с выключенным вентилятором аэродинамической трубы. По-видимому, сказывалось торможение колебаний, вызванное аэродинамическим сопротивлением движущегося в спокойном воздухе стабилизатора.
4. Результаты экспериментальных измерений. В наших экспериментах амплитуда затухающих колебаний А уменьшалась с 0.1 до 0.03 рад, а амплитуда угловой скорости Аи уменьшалась с 16 до 5 с-1. Безразмерная амплитуда угловой скорости АиЬ/(2^и) менялась в пределах от 0.9 до 0.1.
Числа подобия в явлении колебания тел в потоке можно получить обезразмери-ванием уравнений движения (1) или (3). В дополнение к числам подобия Рейнольдса и Струхаля числами подобия являются логарифмический декремент колебаний 6 и рвЬ3/I,. ^ Скрутон в своей работе [8] впервые предложил в качестве критерия подобия для вращательных колебаний моста комбинацию из чисел подобия 61\ /(рЬ4), где — момент инерции единицы длины моста; Ь — характерный поперечный размер моста (ширина). Позднее в работе [9] было предложено использовать аналогичное число подобия для поступательных колебаний моста. Число подобия для поступательных колебаний впоследствии получило название числа Скрутона. Как для вращательных, так и поступательных колебаний число Скрутона определяет критическую скорость воздушного потока, при превышении которой под действием потока возникают колебания моста.
-и.4и ^—I—|—I—|—.—|—I—|—I—|—I—|—I—,
6 8 10 12 14 16 18 20
V, м/с
Рис. 3. Зависимость коэффициента затухания п от скорости набегающего потока V: 1, 3 — 4 = 0.625Д, д = 0.65Д; 2, 4 — 4 = 0.875Д, д = 0.65Д; 1, 2 — эксперимент; 3,4 — линейная аппроксимация.
В нашей работе изучались тела, для которых по образцу числа Скрутона можно сформулировать и применять критерий подобия при моделировании вращательных колебаний [10]:
^ = (8)
Число Рейнольдса, определенное по диаметру цилиндра, находилось в диапазоне от 0.52 • 105 до 1.3 • 105. Число Струхаля Sh при изменении скорости потока принимало значения от 0.007 до 0.015. Логарифмический декремент колебаний упругой подвески с подвешенным цилиндром при нулевой скорости набегающего потока 5 = 0.066. Число Скрутона, вычисленное по формуле (8), в наших экспериментах равно 2.7.
На рис. 3 и 4 представлены зависимость коэффициента затухания п от скорости набегающего потока V и зависимость частоты колебаний ф от скоростного напора д.
Эти графики демонстрируют близость экспериментальных зависимостей к линейным функциям.
Результаты определения вращательных производных т" + т®, представлены в табл. 1. Приведены также 95 %-ные доверительные интервалы.
Таблица 1. Вращательные производные т" + тЯ, цилиндра с диском и со стабилизатором
<г/£> 9/0
0.50 0.65 0.80
0.625 -2.50 ±0.06 -2.36 ±0.06 -2.30 ±0.06
0.75 -2.27 ±0.06 -2.43 ±0.06 -2.38 ±0.06
0.875 -1.93 ±0.06 -2.06 ±0.06 -2.27 ±0.06
Вращательная производная цилиндра со стабилизатором, но без соосного диска, равна —2.10±0.06. Соосный диск, диаметр которого меньше диаметра цилиндра,
7.6-
—i—1—i—1—i—1—i—1—i—1—i—'—i—1—i—1—i—1—i—1 20 40 60 80 100 120 140 160 180 200 220
q, Па
Рис. 4. Зависимость частоты колебаний up от скоростного напора q: 1,3 — цилиндр без соосного диска; 2, 4 — d = 0.875D, g = 0.65D; 1, 2 — эксперимент; 3, 4 — линейная аппроксимация.
ускоряет затухание вращательных колебаний [1], но стабилизатор ускоряет затухание колебаний гораздо успешнее.
Простая математическая модель позволяет объяснить влияние стабилизатора на затухание колебаний. Модель предполагает, что аэродинамический момент цилиндра со стабилизатором складывается из момента, действующего на цилиндр без стабилизатора, и момента, создаваемого нормальной силой, действующей на стабилизатор. Принимается, что центр давления стабилизатора находится на расстоянии Li от центра масс цилиндра. Вращательное движение цилиндра с угловой скоростью в приводит к изменению угла атаки стабилизатора на малую величину Да = в Li/v. Изменению угла атаки соответствует изменение коэффициента нормальной силы стабилизатора cy, который для стабилизатора, являющегося плоской пластиной, пропорционален изменению угла атаки cy = кДа. Необходимо учесть, что стабилизатор экранируется цилиндром с размещенным в носовой части соос-ным диском. Допустим, что экранирование можно учесть введением коэффициента П < 1 в выражение для момента. Коэффициент п принимает разные значения для цилиндра без диска и для дисков с различным диаметром. Дополнительный момент Mzs, связанный с вращательным движением стабилизатора, выражается формулой
Mzs = nqsicy Li, Cy = xeLi/v.
Влияние стабилизатора можно учесть, прибавив в правую часть уравнения движения (1) выражение для момента Mzs.
Колебания аэродинамического момента, пропорционального производной в угла наклона по времени, по фазе сдвинуты на п/2 от колебаний угла наклона цилиндра в. Работа, совершаемая дополнительным аэродинамическим моментом за время одного периода колебаний, тратится на уменьшение энергии колебательного движения, т.е. на затухание колебаний. Таким образом, влияние стабилизатора сводится к добавлению постоянной величины к вращательной производной m" + mez.
Результаты определения вращательных производных mzz представлены в табл. 2. В таблицу включены 95 %-ные доверительные интервалы.
Таблица 2. Вращательные производные mez цилиндра с диском и со стабилизатором
d/D g/D
0.50 0.65 0.80
0.625 — 1.59 ± 0.11 -1.51 ±0.11 — 1.51 ± 0.11
0.75 — 1.45 ± 0.11 -1.34 ±0.11 — 1.51 ± 0.11
0.875 — 1.72 ± 0.11 -1.50 ±0.11 — 1.50 ± 0.11
Вращательная производная mzz цилиндра со стабилизатором, но без соосного диска, равна —2.67 ± 0.11. Вращательная производная mez цилиндра со стабилизатором без соосных дисков по абсолютной величине значительно больше mz цилиндров с соосными дисками.
Проводя рассуждения, аналогичные изложенным выше для вращательных производных mz + mzz, можно получить, что влияние стабилизатора на вращательные производные mzz сводится к добавке некоторой величины к производным цилиндра без стабилизатора. При этом положительные производные становятся отрицательными. Цилиндр со стабилизатором обладает статической устойчивостью.
5. Заключение. Определены характеристики затухающих вращательных колебаний цилиндра с удлинением девять со стабилизатором в кормовой части и со-осным диском, установленным в головной части. Если для удобообтекаемых тел для описания влияния потока на затухание вращательных колебаний используются вращательные производные коэффициента момента, то для плохо обтекаемых тел характеризовать влияние потока удобно аналогами вращательных производных, которые в строгом смысле термина вращательными производными не являются. Аналоги вращательных производных для плохо обтекаемых тел зависят от диапазона амплитуды вращательных колебаний и от диапазона амплитуды угловой скорости. Вращательные производные определены для цилиндра с соосными дисками различных диаметров, расположенных на различных расстояниях между диском и передним торцом цилиндра. Предложена простая модель влияния стабилизатора на затухание колебаний.
Литература
1. Рябинин А. Н., Кауфман Д. В. Определение вращательных производных цилиндра с со-осно установленным диском в воздушном потоке. Вестник Санкт-Петербургского университета.. Математика. Механика. Астрономия 8 (66), вып. 1, 158—166 (2021). https://doi.org/10.21638 /spbu01.2021.114
2. Koenig K., Roshko A. An experimental study of geometrical effects on the drag and flow field of two bodies separated by a gap. J. Fluid Mech. 156 167—204 (1985).
3. Устройство для перевозки крупногабаритных грузов на внешней подвеске вертолета. Пат. Рос. Федерации 2002676; заяв. 10.12.1990; опубл. 15.11.1993. Бюл. №41-42.
4. Рябинин А. Н., Тюрин Б. Ф. Поведение груза, подвешенного под вертолетом. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 1, 87-91 (1993).
5. Ковалев М. А. О расчете и исследовании аэродинамических труб. Уч. зап. Ленингр. ун-та, 7, 61-86 (1939).
6. Белоцерковский С.М., Скрипач Б. К., Табачников В. Г. Крыло в нестационарном потоке газа. Москва, Наука (1971).
7. Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. 4-е изд. Москва, Наука (1974).
8. Scruton C. An experimental investigation of the aerodynamic stability of suspension bridges IABSE congress report 3, 463-473 (1948).
9. Curtis D. J., Hart J. J., Scruton C. An aerodynamic investigation for the suspended structure of the proposed Tsing Ma bridge. Eng. Struct. 7, 23-34 (1985).
10. Рябинин А.Н., Велигжанин А. А. Вращательные и поступательные колебания цилиндров малого удлинения в воздушном потоке. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 9 (67), вып. 4, 729-739 (2022).
Статья поступила в редакцию 5 апреля 2023 г.;
доработана 5 июня 2023 г.; рекомендована к печати 31 августа 2023 г.
Контактная информация:
Кауфман Даниил Витальевич — аспирант: kifsdan@mail.ru
Рябинин Анатолий Николаевич — д-р физ.-мат. наук, ст. науч. сотр., проф.; a.ryabinin@spbu.ru
Damping oscillations of a cylinder with a coaxial disc and a stabilizer
D. V. Kaufman, A. N. Ryabinin
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Kaufman D. V., Ryabinin A.N. Damping oscillations of a cylinder with a coaxial disc and a stabilizer. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11 (69), issue 1, pp. 161-170. https://doi.org/10.21638/spbu01.2024.111 (In Russian)
The damped rotational oscillations of the cylinder, which is equipped with a coaxial disk in the head part, and has a stabilizer in the tail part, are studied. The elongation of the cylinder (the ratio of length to diameter) is nine. The cylinder is mounted in the test section of the low velocities wind tunnel with a wire suspension containing steel springs. In the equilibrium position, the axis of the cylinder is horizontal and parallel to the velocity vector of the incoming flow. Under the action of the air flow, the attenuation of rotational oscillations of the cylinder changes. A semiconductor strain gauge is attached to one of the suspension springs, which measures the dependence of spring tension on time during oscillations. The voltage at the output of the strain gauge is transmitted to the PC oscilloscope. The digital signal of the oscilloscope is transmitted to the computer. After calibration of the device, the frequency and amplitude of damped rotational oscillations around a horizontal axis passing through the center of the cylinder and perpendicular to the velocity vector of the incoming flow are determined. The effect of the air flow is described by analogues of rotational derivatives, which in the case of bluff bodies depend on the amplitude of the oscillations of the angle of inclination of the body and on the amplitude of the angular velocity. A simple model of the stabilizer's effect on rotational derivatives is proposed. Keywords: wind tunnel, rotational vibrations, strain gauge, cylinder.
References
1. Ryabinin A.N., Kaufman D.V. Determination of rotational derivatives of a cylinder with a coaxially mounted disk in an air flow. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 8 (66), iss. 1, 158-166 (2021). https://doi.org/10.21638/spbu01.2021.114 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 55 (4), iss. 1, 497-503 (2022). https://doi.org/10.1134/S1063454122040197].
2. Koenig K., Roshko A. An experimental study of geometrical effects on the drag and flow field of two bodies separated by a gap. J. Fluid Mech. 156, 167-204 (1985).
3. The device for the transport of bulky cargo on external suspension of a helicopter (Patent of Russian Federation no. 2002676; appl. 10.12.1990; publ. 15.11.1993. Bull. no. 41-42. (In Russian)
4. Ryabinin A.N., Tyurin B. F. The behavior of the load suspended under a helicopter Vestnik of Saint Petersburg University. Series 1. Mathematics. Mechanics. Astronomy 1, 87—91 (1993). (In Russian)
5. Kovalev M. A. On calculation and investigation of wind tunnels. Uchenye zapiski Leningradskogo universiteta 7, 61—86 (1939). (In Russian)
6. Belotserkovkii S. M., Skripach B. K., Tabachnikov V. G. A wing in an unsteady gas stream. Moscow, Nauka Publ. (1971). (In Russian)
7. Bogoliubov N. N., Mitropolski Y.A. Asymptotic methods in the theory of non-linear oscillations New York, Gordon and Breach (1961). [Rus. ed.: Bogoliubov N. N., Mitropolski Y.A. Asimptoticheskie metody v teorii nelineinykh kolebanii. 4th. Nauka Publ. (1974). (In Russian)].
8. Scruton C. An experimental investigation of the aerodynamic stability of suspension bridges IABSE congress report 3, 463-473 (1948).
9. Curtis D. J., Hart J. J., Scruton C. An aerodynamic investigation for the suspended structure of the proposed Tsing Ma bridge. Eng. Struct. 7, 23-34 (1985).
10. Ryabinin A.N., Veligzhanin A.A. Rotational and translational oscillations of small-aspect-ratio in the air flow. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 9 (67), iss. 4, 729-739 (2022). https://doi.org/10.21638/spbu01.2022.414 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 55 (4), iss. 4, 489-496 (2022). https://doi.org/10.1134/S106345412204015X].
Received: April 5, 2023 Revised: June 5, 2023 Accepted: August 31, 2023
Authors' information:
Daniil V. Kaufman — kifadan@mail.ru Anatoly N. Ryabinin — a.ryabinin@spbu.ru