Научная статья на тему 'Динамика качающейся пружины с подвижным подвесом'

Динамика качающейся пружины с подвижным подвесом Текст научной статьи по специальности «Физика»

CC BY
589
438
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕЛИНЕЙНЫЕ КОЛЕБАНИЯ / КАЧАЮЩАЯСЯПРУЖИНА / ПОДВИЖНЫЙ ПОДВЕС / АВТОПАРАМЕТРИЧЕСКИЙ РЕЗОНАНС / NONLINEAR OSCILLATIONS / SWINGING SPRING / MOVING SUPPORT / AUTOPARAMETRIC RESONANCE

Аннотация научной статьи по физике, автор научной работы — Алдошин Г. Т., Яковлев С. П.

Рассматривается аналитическая модель динамики качающейсяпружины с подвижным подвесом. При отношении частот колебаний 2:1 :1 и 2:4 :1 в системе возникают резонансы третьего и четвертого порядков. Исследование проводится методом нормальной формы Пуанкаре—Биркгофа. По сравнению с неподвижным подвесом установлено сокращение периода обмена энергией между модами и изменение амплитуд колебаний.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Dynamics of a Swinging Spring With Moving Support

We consider the analytical model of the dynamics of a swinging spring with moving support. With respect to the vibration frequencies of 2:1:1 and 2:4:1 resonances occur in the third and fourth order. The study is carried out by means of the Poincar´e—Birkhoff. normal form. It is revealed that theperiod of energy exchangebetween the modes reduces compared with a structure with fixed suspension and oscillation amplitudes change.

Текст научной работы на тему «Динамика качающейся пружины с подвижным подвесом»

УДК 534.1

Вестник СПбГУ. Сер. 1. 2012. Вып. 4

ДИНАМИКА КАЧАЮЩЕЙСЯ ПРУЖИНЫ С ПОДВИЖНЫМ ПОДВЕСОМ*

Г. Т. Алдошин1, С. П. Яковлев2

1. Балтийский государственный технический университет «ВОЕНМЕХ» им. Д.Ф.Устинова, д-р техн. наук, профессор, [email protected]

2. Балтийский государственный технический университет «ВОЕНМЕХ» им. Д.Ф.Устинова, студент, [email protected]

1. Введение. В 1928 г. советскими физиками Г. Ландсбергом и Л. Мальденшта-мом одновременно с индийским физиком Ч. В. Раманом был открыт эффект комбинационного рассеивания в спектре углекислого газа. Это необычное и загадочное явление взаимодействия различных колебательных мод между собой, сопровождающееся перекачкой энергии между модами, с квантовомеханических позиций было объяснено в 1931 г. итальянским физиком Э. Ферми [1] и получило название резонанса Ферми. Поскольку сам факт возникновения комбинационных тонов в нелинейной акустике был известен, и классическое описание движения оказывает существенную помощь при анализе квантовомеханических результатов, Л. Мальденштам предложил своим молодым сотрудникам А. Витту и Г. Горелику исследовать резонанс Ферми методами классической механики. Результаты работы были опубликованы в 1933 г. [2]. В модели молекулы СО2 атомы рассматривались как силовые центры, взаимодействующие друг с другом с силами, зависящими от расстояний, а силы моделировались невесомыми пружинами; молекула была схематизирована системой двух математических упругих маятников, скрепленных с атомом углерода. Было установлено, что характер колебаний в системе зависит от соотношения частот вертикальных и горизонтальных колебаний, и при их соотношении 1: 2 имеют место сильное взаимодействие между модами колебаний и периодическая перекачка энергии между ними, скорость и глубина которой зависят от начальных условий. Маятник совершает при этом модулированное колебание, т.е. происходит расщепление частот. Качественно выводы исследования Витта—Горелика совпали с расчетами Ферми. Сам процесс колебаний Витт и Горелик назвали параметрическим резонансом. Но обычно под параметрическими понимают колебания, вызываемые изменением параметров системы (массы, длины, жесткости и т.п.). В данном случае внешнее воздействие отсутствует, механизм возбуждения заключен в самой системе (нелинейные обратные связи), поэтому по аналогии с автоколебаниями более предпочтителен термин «автопараметрический резонанс», предложенный Минорским [3]. В теории гамильтоновых систем такое соотношение частот называют резонансом третьего рода [4]. В дальнейшем Витт и Горелик к своей задаче более не возвращались. Специалистов по нелинейным колебаниям гамильтоновых систем резонансные случаи долгое время не интересовали, что можно объяснить как трудностью проблемы, так и отсутствием практических приложений. Положение изменилось в 50-е годы прошлого века, что было вызвано необходимостью решения возникших в то время прикладных задач, в первую очередь, динамики

*Доклад на Международной конференции по механике «Шестые Поляховские чтения» 31 января— 3 февраля 2012 г., Санкт-Петербург, Россия.

© Г. Т. Алдошин, С. П. Яковлев, 2012

летательных аппаратов и водного транспорта. Уравнения движения таких объектов представляют собой сложную нелинейную систему пространственных уравнений, и для возможности их решения проводилась линеаризация уравнений с разделением на продольное и боковое движение. Но с накоплением опыта эксплуатации сверхзвуковых самолетов и быстроходных судов выяснилось, что в том случае, когда колебания по тангажу происходят с частотой в два раза большей, чем по крену, энергия продольных колебаний передается боковым и возникает автопараметрический резонанс [5]. В настоящее время задаче нелинейных колебаний автономных систем с двумя степенями свободы посвящено много работ (библиографию см. [4, 6-7]). Обрела второе дыхание и задача Витта—Горелика, получившая название «качающаяся пружина» [8-16]. Обзор полученных результатов содержится в [16].

Настоящая статья рассматривает ранее не изучавшиеся колебания качающейся пружины с подвижным подвесом и представляет уточнение модели Витта—Горелика колебаний молекулы СО2.

2. Постановка задачи. Основные уравнения. Материальная точка массы т подвешена на невесомой пружине жесткости с к подвесу массы т\, который может смещаться по горизонтальной оси и удерживается между упорами пружинами жесткости С1. Начало координат выбрано в положении устойчивого равновесия системы: маятник находится в нижнем вертикальном положении, подвес занимает среднее положение между упорами (рис. 1, а). Система имеет три степени свободы: Х1 —отклонение маятника от положения равновесия, у1 — отклонение от подвеса, у2 — смещение подвеса. На рисунке 1, б приведена механическая модель линейной молекулы СО2: атомы С, О с массами тс, то расположены по одной линии, атом С — в центре между атомами О; силы взаимодействия между атомами моделируются линейными безмассовыми пружинами жесткости кч при растяжении-сжатии и жесткости к^ при изгибе. Линейная молекула СО2 имеет четыре колебательных степени свободы; соответствующие формы колебаний и их частоты показаны на рис. 1, б. Но так как изгибные колебания могут происходить как в плоскости рисунка, так и в плоскости, ей перпендикулярной, изгибная частота ^ дважды вырождена (изгибные колебания имеют две степени свободы, но одну частоту). Два других колебания происходят в осевом направлении: симметричные колебания VI и антисимметричные ^3. Так как колебания атомов молекулы совершаются относительно ее центра масс, точка подвеса (атом С) смещается, что не учитывалось в модели Витта—Горелика.

Дальнейшее рассмотрение задачи удобно проводить в безразмерных координатах, выбрав в качестве линейного масштаба равновесную длину маятника Ь = ¿о + тд/с (/о —длина ненапряженной пружины) и масштаба времени \/Ь/д. Тогда

XI У1 У2 т Гд, Х1 = Т> = У2 = Т' г = у Т1- (1)

За масштаб энергии принимаем величину тдЬ; кинетическую энергию Т и потенциальную энергию П запишем как

Т = ^2 + (у1+у2)2+Иу2}, П = 1а2[(1-а)2-(1-а)2 + 2ру22]-х1, (2)

где 1 = \/(Ь + Ж1)2 + у\/Ь—длина пружины, /л = т\/т, а2 = к2/к2, к2 = с/т, а = ¡о/Ь, в = с1/с (точка над переменными означает дифференцирование,

Рис. 1. Расчетные схемы: а) качающаяся пружина с подвижным подвесом; б) механическая модель молекулы СО2.

черта над переменными будет опущена; значение потенциальной энергии в положении статического равновесия принято нулевым). Гамильтониан системы Н = Т + П после несложных преобразований можно записать в виде

Я = + (г/1 + т? + т + а2(12- 1) - 2(<т2 -!)(/-!)+ 2/За2у22] - Х1. (3)

Уравнения Лагранжа получаются такими: х1 + а2х 1 + (а2 - 1) (1 -

1 + Х\

I

1-1 Т

У1 + (1 + м)У2 +2ва2 у2 =0.

У1 + VI + VI ~ (сг2 - 1)—^—2/1 = 0,

(4)

Линеаризация уравнений (4) в окрестности положения равновесия приводит к системе линейных уравнений

Х1 + а2Х1 = 0,

У1 + У2 + У1 =0, (5)

У1 + (1 + м)У2 +2ва2 у2 =0.

В отличие от колебаний пружины с постоянной точкой подвеса движения по отдельным координатам не разделяются, и уравнения (5) описывают колебательный процесс с частотами главных колебаний:

IБ ± ^Б2 - 8/1/3<т2

з = уг-, й' = 2/За2 + (1 + /г).

(6) 47

0

Уиттекер [17] был, по-видимому, первым, кто обратил внимание на то, что отбрасывание членов высших порядков малости при линеаризации может не только сказаться на точности решения, но и исказить физический характер явления, так что приближенное решение даже качественно может не соответствовать точному.

В теории устойчивости система (5), являющаяся линейной частью системы (4), соответствует критическому случаю, в котором движение определяется структурой нелинейных членов. Еще с ранних работ Н. М. Ляпунова известно, что движение автономной системы имеет тот или иной вид в зависимости от того, существует ли между частотами Qi линейной системы (5) целочисленная линейная зависимость

п1П1 + П2^2 + П3П3 = 0. (7)

Если сумма модулей целых чисел

Ы + Ы + Ы = к, (8)

то в системе имеет место внутренний резонанс к-го порядка.

При решении конкретной механической задачи в резонансных случаях частоты линейной системы зависят от параметров вполне определенным образом. В таблице приведены значения параметров системы, при которых реализуются резонансы третьего и четвертого порядка, т.е. отношения частот ^2, О3 близки или равны 2:1:1 и 2:4:1 соответственно.

Порядок резонанса к 3 : 1 - 2 + ( — 1) • 1 + ( — 1) • 1 = 0 4 : (-1) -2 + 1 • 4 + (-2) -1 = 0

Жесткости пружин с,с±, Н/м с = 20; С1 = 10 с = 20; С1 = 27

Массы 771,7711, кг пг = 1; ГП1 = 4, 5 т=1;т1= 0,82

Равновесная длина Ь, м 1,962 1,962

Парциальные безразмерные частоты

Частота колебаний

груза на пружине а 2 2

Частота боковых

колебаний маятника 1 1

Частота колебаний

подвеса 2(3/{¡л,1)<т 0,853 2,436

Главные безразмерные частоты Пх = 2; П2 = 1, 237; П3 = 0, 762 Пх = 2; П2 = 3, 806; П3 = 0, 954

Начальные условия жх(0) = 0, 1; ш(0) = 0,02; 2/2(0) = 3,75 • 10~3 Х!(0) = т(0) = у2(0) = 0 жх(0) = 0, 1; 2/1(0) = 0,02; 2/2(0) = -8,32 • КГ3 ¿1(0) = 2/1(0) = 2/2(0) =0

Результаты численного решения при резонансе порядка к = 3 для начальных условий, заданных в таблице, приведены на рис. 2, а. Там же (рис. 2, б), для сравнения, приведено решение линейной системы (5). Их сравнение наглядно иллюстрирует нелинейные эффекты: зависимость частоты от амплитуды и перекачку энергии между модами колебаний.

3. Анализ задачи методом инвариантной нормализации гамильтониана. Известно, что в нерезонансном случае решение задачи может быть получено асимптотическими методами. В классической и небесной механике соизмеримость частот вызывает появление «вековых членов», так что ряды в теории возмущений расходятся.

Для исследования рассматриваемой задачи используем метод приведения к нормальной форме Пуанкаре—Биркгофа, получивший в последнее десятилетие широкое

Рис. 2. Численное решение для резонанса к = 3: а) решение исходной системы (4); б) решение линеаризованной системы (5).

применение в разнообразных нелинейных задачах. При этом гамильтонова система упрощается настолько, что становится возможным получение решения до членов третьего порядка и выше. Ниже используется алгоритм Журавлева—Петрова [14], который пригоден для всех случаев: резонансный, нерезонансный, автономный, неавтономный.

Вводим нормальные координаты

х1 ,

у

У 2 ~ 73 2/1 72-73

722/1 ~ У 2 72-73

(9)

где 72,3 = П з/(2ва2 — (м + 1)^2 з) —коэффициенты распределения амплитуд между второй и третьей степенями свободы.

Представим гамильтониан в виде ряда в окрестности равновесия по степеням малого параметра е:

Н = Но + еЛ + е2Л + ...,

Но = + сг2.г'2 + а2(г>2 + С%2) + а3(го2 + П2г2)},

21

F1 = -(a2-l){x(y + z)2},

1, 2 Г (у + г)4

4

— х2(у + г)2

(10)

где а2,з = 1 + 272,3 + (м + 1)72 з —главные квазиинерционные коэффициенты.

Искомые решения х, у, г и частоты колебаний также представим в виде рядов по степеням е:

х = хо + ех1 + .. ., П1 = Пю + еПц + . ...

у = уо + еу1 + .. ., П2 = П20 + е^21 +

г = го + ег1 + . .., Пз = П30 + еПз1 +

(11)

х

г

j*$'D')+#&&9%&$#')$(I%(H#((#&'$&'+'(+H'

Тогда после подстановки (11) в (4) получим систему рекуррентных уравнений

Хо + хо = 0, Уо + ^Уо = 0, х1 + П2иХ1 = -а 1(у0 + г0)2, У1 + Щ1У1 = ---х0(Уо + ¿о),

a2

¿o + O^zo = 0, ¿1 + Q^zi = ----х0(у0 + z0).

2 Л (12)

0

а2- 1

аз

Поскольку система (12) автономна, за начальный момент времени t = 0 можно выбрать такой, когда x(0) = a, y(0) = b. Тогда начальными условиями для системы (12) будут

xo(0) = a, xi(0)=0, ..., xi(0) = 0, ( )

yo(0)= b, yi(0) = 0, ..., yi(0)=0. ( 3)

Все остальные начальные возмущения приняты нулевыми. Анализируя уравнения первого приближения после подстановки в них решения порождающей системы, можно показать, что при значениях главных частот ^i = 2, П2 = О3 = 1 наблюдается внутренний резонанс порядка k = 3 между всеми парциальными системами, который и будем исследовать методом инвариантной нормализации [14].

Получив решение системы уравнений для невозмущенной части гамильтониана Ho, подставим его в рекуррентную формулу нормальной формы и ограничимся первым приближением:

_ _ q

F = F1 = -[-X(V + W)2 + U (V + W)(Y + Z)+ X(Y + Z)2]. (14)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

8

Здесь X, U, Y, V, Z, W — медленно меняющиеся амплитуды координат и импульсов соответственно. Нормальную форму (14) удобно анализировать в комплексных переменных Биркгофа z\ = l/\/2U + а/2iX, z2 = V + iY, z3 = W + iZ\

F = F = 2(F=-^[Z1(Z2 + Z3)2 -Z^Zi+ZS)2}. (15)

16

Система уравнений первого приближения

dZx 16 v л>

Z2 = ^- = -^Z1(Z2 + Z3) (16)

aZ 2 о

имеет интегралы

Ho = const ^ 21Zi|2 + a2|Z2|2 + аз^з|2 = C2,

F = const, (17)

G2 = (Z2 - Z3)(Z2 - Z3) = const = (y0 - ¿o)2 + (y0 - z0)2.

С использованием последних можно составить уравнение для переменной £ = |Z2 + Z3|2:

§ = O-0-JL,. (18)

Первый интеграл уравнения (18) имеет форму интеграла энергии

Здесь П = 9/8(2/3£3 — 2С'2£2) —потенциал. Таким образом, перешли к аналогии потенциальной ямы, исследуя которую, можно получить зависимость периода перекачки энергии между различными модами колебаний от начальных условий при любом конечном отклонении энергии Е от минимума. Наиболее ярко данный процесс проявляется, когда уо ^ хо ^ хо (рис. 3). Такая неустойчивость системы при малом возмущении сепаратрисы получила название процесса срыва вертикальных колебаний [9].

О 20 40 60 80 100

Рис. ,3. Графики y(t),z(t) (результаты численного решения точных уравнений) и их огибающие Y(t),Z(t) при срыве вертикальных колебаний.

Аналитически найденные выражения для огибающих имеют вид

X(t) = x0(|th(r - 7-о/2)| + |th(r - 3tq/2)I - 1), Y(t) + Z(t) = 2v/2.T0(sch(r - rn/2) + sch(r - 3r0/2)), (2Q)

r = (3/2xn)i, то = 8, 478, T= — = 56,488.

3xo

Качественно картина явления остается той же, как и в случае закрепленного подвеса. Количественные отличия выражаются в сокращении периода обмена энергией между различными степенями свободы системы Т примерно в два раза и в том, что максимум амплитуды суммарного бокового колебания в 2а/2 (не в 2, как в случае неподвижного подвеса [14, 15]) раза превышает максимум амплитуды вертикального, в том числе и когда y(0) = 0 или z(0) = 0. Оба эти отличия вызваны, по всей видимости, подвижностью подвеса, и, как следствие, наличием в системе переносного движения.

В следующей статье будет дано описание нелинейных колебаний молекулы углекислого газа с использованием более точной механической модели (рис. 1, б).

Литература

1. Fermi E. "Über der Ramanefect des kohlendioxyds // Z. Phys. 1931, 71. N2. S. 250-259.

2. Витт А., Горелик Г. Колебания упругого маятника как пример колебаний двух параметрически связанных линейных систем // ЖТФ. 1933. Т. III. Вып. 2-3. С. 294-307.

3. Minorsky N. Nonlinear Oscillations. New York: D. V. Nastrland Company Inc., 1962.

4. Маркеев А. П. Устойчивость гамильтоновых систем // Нелинейная механика / под ред. В. М. Матросова, А. В. Карапетяна. М.: Физматлит, 2001. С. 114-130.

5. Алдошин Г. Т. К вопросу о линеаризации уравнений Лагранжа // Пятые Поляхов-ские чтения. Избранные труды Международной научной конференции по механике. Санкт-Петербург: 3-6 февраля 2009 г. СПб., 2009. С. 39-46.

6. Куницын А. Л., Маркеев А. П. Устойчивость в резонансных случаях // Итоги науки и техники. Серия Общая механика. М.: ВИНИТИ, 1979. Т. 4. С. 58-139.

7. Холостова О. В. О нелинейных колебаниях гамильтоновой системы с двумя степенями свободы при резонансе третьего порядка // ПММ. 2010. Вып. 5. С. 789-811.

8. Найфе А. Х. Методы теории возмущений. М.: Мир, 1976. 535 с. (Nayfeh A. H. Perturbation Methods. J. Wiley, 1973.)

9. Старжинский В. М. Прикладные методы нелинейных колебаний. М.: Наука, 1977. 235 с.

10. Богаевский В. Н., Повзнер А. Я. Алгебраические методы в теории возмущений. М.: Наука, 1987. 255 с.

11. Цельман Ф. Х. О «перекачке» энергии между связанными осцилляторами в случае резонанса третьего порядка // ПММ. 1970. Т. 14. Вып. 5. С. 957-962.

12. Маркеев А. П., Чеховская Т. А. О колебаниях упругого маятника // Изв. РАН. МТТ. 1998. №2. С. 18-26.

13. Зарипов М. Н., Петров А. Г. Нелинейные колебания качающейся пружины // Докл. Акад. наук. 2004. Т. 399, №3. С. 347-352.

14. Петров А. Г. Нелинейные колебания качающейся пружины при резонансе // Изв. РАН. МТТ. 2006. №5. С. 18-28.

15. Петров А. Г., Фомичев А. В. О нелинейных трехмерных колебаниях тяжелой материальной точки на пружине при резонансе // Изв. РАН. МТТ. 2008. №5. С. 15-26.

16. Петров А. Г., Шундерюк М. М. О нелинейных колебаниях тяжелой материальной точки на пружине // Изв. РАН. МТТ. 2010. №2. С. 15-26.

17. Уиттекер Е. Т. Аналитическая динамика. М.; Л.: ОНТИ, 1937. 500 с. — (Whit-taker E. T. A Treatise on the Analytical Dynamics of Particles and Rigid Bodies. Cambridge, Cambridge Univ. Press, 1927. 430 p.)

Статья поступила в редакцию 26 июня 2012 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.