Двойной маятник

редактировать
Двойной маятник состоит из двух маятников, прикрепленных встык.

In физика и математика, в области динамических систем, двойной маятник является маятником с другим маятником. ndulum, прикрепленный к его концу, представляет собой простую физическую систему, которая демонстрирует богатое динамическое поведение с сильной чувствительностью к начальным условиям. Движение двойного маятника управляется набором связанных обыкновенных дифференциальных уравнений и является хаотическим.

Содержание

  • 1 Анализ и интерпретация
    • 1.1 Лагранжиан
  • 2 Хаотическое движение
  • 3 См. Также
  • 4 Примечания
  • 5 Ссылки
  • 6 Внешние ссылки

Анализ и интерпретация

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

Двойной составной маятник Движение двойного составного маятника (из численного интегрирования уравнений движения) Траектории двойного маятника

В составном маятнике масса равна распределены по его длине. Если масса распределена равномерно, то центр масс каждой конечности находится в ее средней точке, и конечность имеет момент инерции I = 1/12 мл относительно этой точки.

Удобно использовать углы между каждой конечностью и вертикалью в качестве обобщенных координат, определяющих конфигурацию системы. Эти углы обозначены θ 1 и θ 2. Положение центра масс каждого стержня можно записать в терминах этих двух координат. Если принять начало декартовой системы координат в точке подвешивания первого маятника, то центр масс этого маятника находится в точке:

x 1 = l 2 sin ⁡ θ 1 Y 1 = - l 2 соз ⁡ θ 1 {\ Displaystyle {\ begin {align} x_ {1} = {\ frac {l} {2}} \ sin \ theta _ {1} \\ y_ {1} = - {\ frac {l} {2}} \ cos \ theta _ {1} \ end {align}}}{\ displaystyle {\ begin {align} x_ {1} = {\ frac {l} {2}} \ sin \ theta _ {1} \\ y_ {1} = - {\ frac {l} {2}} \ cos \ theta _ {1} \ end {align}}}

, а центр масс второго маятника находится в

x 2 = l ( грех ⁡ θ 1 + 1 2 грех ⁡ θ 2) Y 2 = - l (соз ⁡ θ 1 + 1 2 соз ⁡ θ 2) {\ displaystyle {\ begin {align} x_ {2} = l \ left (\ sin \ theta _ {1} + {\ tfrac {1} {2}} \ sin \ theta _ {2} \ right) \\ y_ {2} = - l \ left (\ cos \ theta _ {1} + {\ tfrac {1} {2}} \ cos \ theta _ {2} \ right) \ end {align}}}{\ displaystyle {\ begin {align} x_ {2} = l \ left (\ sin \ theta _ { 1} + {\ tfrac {1} {2}} \ sin \ theta _ {2} \ right) \\ y_ {2} = - l \ left (\ cos \ theta _ {1} + {\ tfrac { 1} {2}} \ cos \ theta _ {2} \ right) \ end {align}}}

Этой информации достаточно, чтобы записать лагранжиан.

Лагранжиан

Лагранжиан равен

L = кинетическая энергия - потенциальная энергия = 1 2 m (v 1 2 + v 2 2) + 1 2 I ( θ ˙ 1 2 + θ ˙ 2 2) - mg (y 1 + y 2) = 1 2 m (x ˙ 1 2 + y ˙ 1 2 + x ˙ 2 2 + y ˙ 2 2) + 1 2 I (θ ˙ 1 2 + θ ˙ 2 2) - мг (y 1 + y 2) {\ displaystyle {\ begin {align} L = {\ text {кинетическая энергия}} - {\ text {потенциальная энергия}} \\ = {\ tfrac {1} {2}} m \ left (v_ {1} ^ {2} + v_ {2} ^ {2} \ right) + {\ tfrac {1} {2}} I \ left ({ {\ dot {\ theta}} _ {1}} ^ {2} + {{\ dot {\ theta}} _ {2}} ^ {2} \ right) -mg \ left (y_ {1} + y_ {2} \ right) \\ = {\ tfrac {1} {2}} m \ left ({{\ dot {x}} _ {1}} ^ {2} + {{\ dot {y}} _ {1}} ^ {2} + {{\ dot {x}} _ {2}} ^ {2} + {{\ dot {y}} _ {2}} ^ {2} \ right) + { \ tfrac {1} {2}} I \ left ({{\ dot {\ theta}} _ {1}} ^ {2} + {{\ dot {\ theta}} _ {2}} ^ {2} \ right) -mg \ left (y_ {1} + y_ {2} \ right) \ end {align}}}{\ displaystyle {\ begin {align} L = {\ text {кинетическая энергия}} - {\ text {потенциальная энергия}} \\ = {\ tfrac {1} {2}} m \ left (v_ {1} ^ {2} + v_ {2} ^ {2} \ right) + {\ tfrac {1} {2}} I \ left ({{\ dot {\ theta}} _ {1}} ^ {2} + {{\ dot {\ theta}} _ {2}} ^ {2} \ right) -mg \ left (y_ {1} + y_ {2} \ right) \\ = {\ tfrac {1} {2}} m \ left ({{\ dot {x}} _ {1}} ^ {2} + {{\ dot {y}} _ {1}} ^ {2} + {{\ dot {x}} _ {2}} ^ {2} + {{\ dot {y}} _ { 2}} ^ {2} \ right) + {\ tfrac {1} {2}} I \ left ({{\ dot {\ theta}} _ {1}} ^ {2} + {{\ dot {\ theta}} _ {2}} ^ {2} \ right) -mg \ left (y_ {1} + y_ {2} \ right) \ end {выравнивается}}}

Первый член - это линейная кинетическая энергия центра массы тел, а второй член - это кинетическая энергия вращения вокруг центра масс каждого стержня. Последний член - это потенциальная энергия тел в однородном гравитационном поле. точечная нотация указывает производную по времени рассматриваемой переменной.

Подстановка координат выше и изменение уравнения дает

L = 1 6 ml 2 (θ ˙ 2 2 + 4 θ ˙ 1 2 + 3 θ ˙ 1 θ ˙ 2 cos ⁡ (θ 1 - θ 2)) + 1 2 mgl (3 cos ⁡ θ 1 + cos ⁡ θ 2). {\ displaystyle L = {\ tfrac {1} {6}} ml ^ {2} \ left ({{\ dot {\ theta}} _ {2}} ^ {2} +4 {{\ dot {\ theta) }} _ {1}} ^ {2} +3 {{\ dot {\ theta}} _ {1}} {{\ dot {\ theta}} _ {2}} \ cos (\ theta _ {1} - \ theta _ {2}) \ right) + {\ tfrac {1} {2}} mgl \ left (3 \ cos \ theta _ {1} + \ cos \ theta _ {2} \ right).}{\ displaystyle L = {\ tfrac {1} {6}} ml ^ {2} \ left ({{\ dot {\ theta }} _ {2}} ^ {2} +4 {{\ dot {\ theta}} _ {1}} ^ {2} +3 {{\ dot {\ theta}} _ {1}} {{\ d ot {\ theta}} _ {2}} \ cos (\ theta _ {1} - \ theta _ {2}) \ right) + {\ tfrac {1} {2}} mgl \ left (3 \ cos \ theta _ {1} + \ cos \ theta _ {2} \ right).}

Есть только одна сохраняющаяся величина (энергия) и нет сохраняющихся импульсов. Два обобщенных импульса могут быть записаны как

p θ 1 = ∂ L ∂ θ ˙ 1 = 1 6 ml 2 (8 θ ˙ 1 + 3 θ ˙ 2 cos ⁡ (θ 1 - θ 2)) p θ 2 = ∂ L ∂ θ ˙ 2 = 1 6 ml 2 (2 θ ˙ 2 + 3 θ ˙ 1 cos ⁡ (θ 1 - θ 2)). {\ displaystyle {\ begin {align} p _ {\ theta _ {1}} = {\ frac {\ partial L} {\ partial {{\ dot {\ theta}} _ {1}}}} = {\ tfrac {1} {6}} ml ^ {2} \ left (8 {{\ dot {\ theta}} _ {1}} + 3 {{\ dot {\ theta}} _ {2}} \ cos ( \ theta _ {1} - \ theta _ {2}) \ right) \\ p _ {\ theta _ {2}} = {\ frac {\ partial L} {\ partial {{\ dot {\ theta}} _ {2}}}} = {\ tfrac {1} {6}} ml ^ {2} \ left (2 {{\ dot {\ theta}} _ {2}} + 3 {{\ dot {\ theta }} _ {1}} \ cos (\ theta _ {1} - \ theta _ {2}) \ right). \ End {align}}}{\ displaystyle {\ begin {align} p _ {\ theta _ {1}} = {\ frac {\ частичное L} {\ partial {{\ dot {\ theta}} _ {1}}}} = {\ tfrac {1} {6}} ml ^ {2} \ left (8 {{\ dot {\ theta} } _ {1}} + 3 {{\ dot {\ theta}} _ {2}} \ cos (\ theta _ {1} - \ theta _ {2}) \ right) \\ p _ {\ theta _ { 2}} = {\ frac {\ partial L} {\ partial {{\ dot {\ theta}} _ {2}}}} = {\ tfrac {1} {6}} ml ^ {2} \ left (2 {{\ dot {\ theta}} _ {2}} + 3 {{\ dot {\ theta}} _ {1}} \ cos (\ theta _ {1} - \ theta _ {2}) \ справа). \ end {align}}}

Эти выражения могут быть инвертированы в получить

θ ˙ 1 = 6 ml 2 2 p θ 1-3 cos ⁡ (θ 1 - θ 2) p θ 2 16-9 cos 2 ⁡ (θ 1 - θ 2) θ ˙ 2 = 6 ml 2 8 p θ 2 - 3 cos ⁡ (θ 1 - θ 2) p θ 1 16 - 9 cos 2 ⁡ (θ 1 - θ 2). {\ displaystyle {\ begin {align} {{\ dot {\ theta}} _ {1}} = {\ frac {6} {ml ^ {2}}} {\ frac {2p _ {\ theta _ {1 }} - 3 \ cos (\ theta _ {1} - \ theta _ {2}) p _ {\ theta _ {2}}} {16-9 \ cos ^ {2} (\ theta _ {1} - \ theta _ {2})}} \\ {{\ dot {\ theta}} _ {2}} = {\ frac {6} {ml ^ {2}}} {\ frac {8p _ {\ theta _ { 2}} - 3 \ cos (\ theta _ {1} - \ theta _ {2}) p _ {\ theta _ {1}}} {16-9 \ cos ^ {2} (\ theta _ {1} - \ theta _ {2})}}. \ end {align}}}{\ displaystyle {\ begin {align} {{\ dot {\ theta}} _ {1}} = {\ frac {6} {ml ^ {2}}} {\ frac {2p _ {\ theta _ {1}} - 3 \ cos (\ theta _ {1} - \ theta _ {2}) p _ {\ theta _ {2}}} {16- 9 \ cos ^ {2} (\ theta _ {1} - \ theta _ {2})}} \\ {{\ dot {\ theta}} _ {2}} = {\ frac {6} {ml ^ {2}}} {\ frac {8p _ {\ theta _ {2}} - 3 \ cos (\ theta _ {1} - \ theta _ {2}) p _ {\ theta _ {1}}} {16 -9 \ соз ^ {2} (\ тета _ {1} - \ тета _ {2})}}. \ Конец {выровнено}}}

Остальные уравнения движения записываются как

p ˙ θ 1 = ∂ L ∂ θ 1 = - 1 2 ml 2 (θ ˙ 1 θ ˙ 2 sin ⁡ (θ 1 - θ 2) + 3 gl sin ⁡ θ 1) p ˙ θ 2 = ∂ L ∂ θ 2 = - 1 2 ml 2 (- θ ˙ 1 θ ˙ 2 sin ⁡ (θ 1 - θ 2) + gl sin ⁡ θ 2). {\ displaystyle {\ begin {align} {{\ dot {p}} _ {\ theta _ {1}}} = {\ frac {\ partial L} {\ partial \ theta _ {1}}} = - {\ tfrac {1} {2}} ml ^ {2} \ left ({{\ dot {\ theta}} _ {1}} {{\ dot {\ theta}} _ {2}} \ sin (\ theta _ {1} - \ theta _ {2}) + 3 {\ frac {g} {l}} \ sin \ theta _ {1} \ right) \\ {{\ dot {p}} _ {\ theta _ {2}}} = {\ frac {\ partial L} {\ partial \ theta _ {2}}} = - {\ tfrac {1} {2}} ml ^ {2} \ left (- {{ \ dot {\ theta}} _ {1}} {{\ dot {\ theta}} _ {2}} \ sin (\ theta _ {1} - \ theta _ {2}) + {\ frac {g} {l}} \ sin \ theta _ {2} \ right). \ end {align}}}{\ displaystyle {\ begin {align} {{\ dot {p}} _ {\ theta _ {1}}} = {\ frac {\ partial L} {\ partial \ theta _ {1}}} = - {\ tfrac {1} {2}} ml ^ {2} \ left ({{\ dot {\ theta}}} _ {1}} {{\ dot {\ theta}} _ {2}} \ sin (\ theta _ {1} - \ theta _ {2}) + 3 {\ frac {g} {l}} \ sin \ theta _ {1} \ right) \\ {{\ dot {p}} _ {\ theta _ {2}}} = {\ frac {\ partial L} {\ partial \ theta _ {2}}} = - {\ tfrac {1} {2}} ml ^ {2} \ left (- {{\ dot {\ theta}} _ {1}} {{\ dot {\ theta}} _ {2}} \ sin (\ theta _ {1} - \ theta _ {2}) + {\ frac {g} {l}} \ sin \ theta _ {2} \ right). \ end {align}}}

Эти последние четыре уравнения представляют собой явные формулы для временной эволюции системы с учетом ее текущего состояния. Невозможно пойти дальше и интегрировать эти уравнения аналитически, чтобы получить формулы для θ 1 и θ 2 как функций времени. Однако можно выполнить это интегрирование численно, используя метод Рунге Кутта или аналогичные методы.

Хаотическое движение

График времени, в течение которого маятник перевернулся, как функция начальных условий Длительная выдержка двойного маятника, демонстрирующего хаотическое движение (отслеживается с помощью светодиода )

двойной маятник совершает хаотическое движение и показывает чувствительную зависимость от начальных условий. Изображение справа показывает количество времени, прошедшее до того, как маятник перевернется, в зависимости от исходного положения при отпускании в состоянии покоя. Здесь начальное значение θ 1 изменяется вдоль направления x от −3 до 3. Начальное значение θ 2 изменяется вдоль направления y, от От −3 до 3. Цвет каждого пикселя указывает, переворачивается ли маятник в пределах:

  • 10√ ⁄ g (зеленый)
  • 100√ ⁄ g (красный)
  • 1000√ ⁄ g (фиолетовый) или
  • 10000√ ⁄ g (синий).
Три двойные маятники с почти одинаковыми начальными условиями расходятся со временем, показывая хаотический характер системы.

Начальные условия, которые не приводят к перевороту в пределах 10000√ ⁄ g нанесены белым цветом.

Граница центральной белой области частично определяется законом сохранения энергии следующей кривой:

3 cos ⁡ θ 1 + cos ⁡ θ 2 = 2. {\ displaystyle 3 \ cos \ theta _ {1} + \ cos \ theta _ {2} = 2.}{\ displaystyle 3 \ cos \ theta _ {1 } + \ cos \ theta _ {2} = 2.}

В пределах области, определенной этой кривой, то есть если

3 cos ⁡ θ 1 + cos ⁡ θ 2>2, {\ displaystyle 3 \ cos \ theta _ {1} + \ cos \ theta _ {2}>2,}{\displaystyle 3\cos \theta _{1}+\cos \theta _{2}>2,}

тогда любой маятник энергетически невозможно перевернуть. Вне этой области маятник может перевернуться, но это сложный вопрос Чтобы определить, когда он перевернется. Аналогичное поведение наблюдается для двойного маятника, состоящего из двух точечных масс, а не двух стержней с распределенной массой.

Отсутствие собственной частоты возбуждения привело к использование систем двойного маятника в конструкциях сейсмостойкости в зданиях, где само здание перевернутый маятник, и вторичная масса соединена, чтобы завершить двойной маятник.

См. Также

Примечания

Ссылки

  • Мейрович, Леонард (1986). Элементы вибрационного анализа (2-е изд.). McGraw-Hill Наука / Инженерия / Математика. ISBN 0-07-041342-8.
  • Эрик У. Вайстейн, Двойной маятник (2005), ScienceWorld (содержит подробные сведения о задействованных сложных уравнениях) и "Двойной маятник »Роба Морриса, Wolfram Demonstrations Project, 2007 (анимация этих уравнений).
  • Питер Линч, Двойной маятник, (2001). (Моделирование Java-апплета.)
  • Северо-Западный университет, Двойной маятник, (Моделирование Java-апплета.)
  • Теоретическая группа астрофизики высоких энергий в UBC, Двойной маятник, (2005).

Внешние ссылки

Последняя правка сделана 2021-05-18 14:06:44
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте