summaryrefslogtreecommitdiff
path: root/physics/double_pendulum
diff options
context:
space:
mode:
authorAndrew <saintruler@gmail.com>2020-04-22 22:22:29 +0400
committerAndrew <saintruler@gmail.com>2020-04-22 22:22:29 +0400
commit5e94dedcc44c77327b37f34792bc1a47a5be60a7 (patch)
treea5734c6e84215a34568abea6363e29fde31145de /physics/double_pendulum
parent693d8b180c690601ce44d640352b0d44151edb76 (diff)
Добавил лабораторные по физике
Diffstat (limited to 'physics/double_pendulum')
-rw-r--r--physics/double_pendulum/images/Double-Pendulum.pngbin0 -> 35512 bytes
-rw-r--r--physics/double_pendulum/images/Double_pendulum_flips_graph.pngbin0 -> 232088 bytes
-rw-r--r--physics/double_pendulum/images/Taipei_101_Tuned_Mass_Damper.pngbin0 -> 653012 bytes
-rw-r--r--physics/double_pendulum/images/mayak.pngbin0 -> 118820 bytes
-rw-r--r--physics/double_pendulum/images/pendulum.pngbin0 -> 25299 bytes
-rw-r--r--physics/double_pendulum/images/taipei_101.jpgbin0 -> 11978 bytes
-rw-r--r--physics/double_pendulum/images/two_states.pngbin0 -> 39778 bytes
-rw-r--r--physics/double_pendulum/images/vectors.pngbin0 -> 28013 bytes
-rw-r--r--physics/double_pendulum/pendulum.nav47
-rw-r--r--physics/double_pendulum/pendulum.snm8
-rw-r--r--physics/double_pendulum/pendulum.tex476
11 files changed, 531 insertions, 0 deletions
diff --git a/physics/double_pendulum/images/Double-Pendulum.png b/physics/double_pendulum/images/Double-Pendulum.png
new file mode 100644
index 0000000..da43df1
--- /dev/null
+++ b/physics/double_pendulum/images/Double-Pendulum.png
Binary files differ
diff --git a/physics/double_pendulum/images/Double_pendulum_flips_graph.png b/physics/double_pendulum/images/Double_pendulum_flips_graph.png
new file mode 100644
index 0000000..537631e
--- /dev/null
+++ b/physics/double_pendulum/images/Double_pendulum_flips_graph.png
Binary files differ
diff --git a/physics/double_pendulum/images/Taipei_101_Tuned_Mass_Damper.png b/physics/double_pendulum/images/Taipei_101_Tuned_Mass_Damper.png
new file mode 100644
index 0000000..7184000
--- /dev/null
+++ b/physics/double_pendulum/images/Taipei_101_Tuned_Mass_Damper.png
Binary files differ
diff --git a/physics/double_pendulum/images/mayak.png b/physics/double_pendulum/images/mayak.png
new file mode 100644
index 0000000..c365d75
--- /dev/null
+++ b/physics/double_pendulum/images/mayak.png
Binary files differ
diff --git a/physics/double_pendulum/images/pendulum.png b/physics/double_pendulum/images/pendulum.png
new file mode 100644
index 0000000..004647f
--- /dev/null
+++ b/physics/double_pendulum/images/pendulum.png
Binary files differ
diff --git a/physics/double_pendulum/images/taipei_101.jpg b/physics/double_pendulum/images/taipei_101.jpg
new file mode 100644
index 0000000..5396e0f
--- /dev/null
+++ b/physics/double_pendulum/images/taipei_101.jpg
Binary files differ
diff --git a/physics/double_pendulum/images/two_states.png b/physics/double_pendulum/images/two_states.png
new file mode 100644
index 0000000..39b4fc6
--- /dev/null
+++ b/physics/double_pendulum/images/two_states.png
Binary files differ
diff --git a/physics/double_pendulum/images/vectors.png b/physics/double_pendulum/images/vectors.png
new file mode 100644
index 0000000..5bc3a65
--- /dev/null
+++ b/physics/double_pendulum/images/vectors.png
Binary files differ
diff --git a/physics/double_pendulum/pendulum.nav b/physics/double_pendulum/pendulum.nav
new file mode 100644
index 0000000..b90e3f7
--- /dev/null
+++ b/physics/double_pendulum/pendulum.nav
@@ -0,0 +1,47 @@
+\headcommand {\slideentry {0}{0}{1}{1/1}{}{0}}
+\headcommand {\beamer@framepages {1}{1}}
+\headcommand {\slideentry {0}{0}{2}{2/2}{}{0}}
+\headcommand {\beamer@framepages {2}{2}}
+\headcommand {\slideentry {0}{0}{3}{3/3}{}{0}}
+\headcommand {\beamer@framepages {3}{3}}
+\headcommand {\slideentry {0}{0}{4}{4/4}{}{0}}
+\headcommand {\beamer@framepages {4}{4}}
+\headcommand {\slideentry {0}{0}{5}{5/5}{}{0}}
+\headcommand {\beamer@framepages {5}{5}}
+\headcommand {\slideentry {0}{0}{6}{6/6}{}{0}}
+\headcommand {\beamer@framepages {6}{6}}
+\headcommand {\slideentry {0}{0}{7}{7/7}{}{0}}
+\headcommand {\beamer@framepages {7}{7}}
+\headcommand {\slideentry {0}{0}{8}{8/8}{}{0}}
+\headcommand {\beamer@framepages {8}{8}}
+\headcommand {\slideentry {0}{0}{9}{9/9}{}{0}}
+\headcommand {\beamer@framepages {9}{9}}
+\headcommand {\slideentry {0}{0}{10}{10/10}{}{0}}
+\headcommand {\beamer@framepages {10}{10}}
+\headcommand {\slideentry {0}{0}{11}{11/11}{}{0}}
+\headcommand {\beamer@framepages {11}{11}}
+\headcommand {\slideentry {0}{0}{12}{12/12}{}{0}}
+\headcommand {\beamer@framepages {12}{12}}
+\headcommand {\slideentry {0}{0}{13}{13/13}{}{0}}
+\headcommand {\beamer@framepages {13}{13}}
+\headcommand {\slideentry {0}{0}{14}{14/14}{}{0}}
+\headcommand {\beamer@framepages {14}{14}}
+\headcommand {\slideentry {0}{0}{15}{15/15}{}{0}}
+\headcommand {\beamer@framepages {15}{15}}
+\headcommand {\slideentry {0}{0}{16}{16/16}{}{0}}
+\headcommand {\beamer@framepages {16}{16}}
+\headcommand {\slideentry {0}{0}{17}{17/17}{}{0}}
+\headcommand {\beamer@framepages {17}{17}}
+\headcommand {\slideentry {0}{0}{18}{18/18}{}{0}}
+\headcommand {\beamer@framepages {18}{18}}
+\headcommand {\slideentry {0}{0}{19}{19/19}{}{0}}
+\headcommand {\beamer@framepages {19}{19}}
+\headcommand {\slideentry {0}{0}{20}{20/20}{}{0}}
+\headcommand {\beamer@framepages {20}{20}}
+\headcommand {\slideentry {0}{0}{21}{21/21}{}{0}}
+\headcommand {\beamer@framepages {21}{21}}
+\headcommand {\beamer@partpages {1}{21}}
+\headcommand {\beamer@subsectionpages {1}{21}}
+\headcommand {\beamer@sectionpages {1}{21}}
+\headcommand {\beamer@documentpages {21}}
+\headcommand {\gdef \inserttotalframenumber {21}}
diff --git a/physics/double_pendulum/pendulum.snm b/physics/double_pendulum/pendulum.snm
new file mode 100644
index 0000000..2caa33d
--- /dev/null
+++ b/physics/double_pendulum/pendulum.snm
@@ -0,0 +1,8 @@
+\beamer@slide {}{2}
+\beamer@slide {}{3}
+\beamer@slide {}{5}
+\beamer@slide {}{13}
+\beamer@slide {}{16}
+\beamer@slide {}{17}
+\beamer@slide {}{18}
+\beamer@slide {}{20}
diff --git a/physics/double_pendulum/pendulum.tex b/physics/double_pendulum/pendulum.tex
new file mode 100644
index 0000000..f0e53e0
--- /dev/null
+++ b/physics/double_pendulum/pendulum.tex
@@ -0,0 +1,476 @@
+\documentclass{beamer}
+
+\usepackage[T2A]{fontenc}
+\usepackage[utf8]{inputenc}
+\usepackage{graphicx}
+\graphicspath{ {./images/} }
+\usepackage{tempora}
+
+\usepackage{amsmath}
+\usepackage[english,russian]{babel}
+
+\usetheme{Madrid}
+
+\title{Двойной математический маятник}
+\author{Гущин Андрей, 131 группа, факультет КНиИТ}
+
+\begin{document}
+
+\maketitle
+
+\begin{frame}
+ \frametitle{Определение}
+
+\begin{columns}
+ \column{0.6\textwidth}
+
+ \textbf{Двойным маятником} называются два скрепленных
+ математических маятника, двигающихся в одной
+ плоскости, причем точка привеса первого маятника
+ неподвижна, а точка привеса второго маятника
+ совпадает с тяжелой материальной точкой первого маятника.
+
+ \column{0.4\textwidth}
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[height=0.5\textheight]{Double-Pendulum.png}
+ \caption{Схема двойного маятника}
+ \label{}
+ \end{figure}
+
+\end{columns}
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Математическая модель}
+
+\begin{columns}
+ \column{0.6\textwidth}
+ Обозначим как $O_1$ и $O_2$ точки привеса, $l_1$ и $l_2$ длины,
+ $Q_1$ и $Q_2$ веса первого и второго маятника соответственно.
+
+ В качестве обобщенных координат этой системы выберем углы
+ $\varphi_1$ и $\varphi_2$, которые составляют соответственно
+ $l_1$ и $l_2$ с вертикалью.
+
+ \column{0.4\textwidth}
+ \begin{figure}[H]
+ \centering
+ \includegraphics{pendulum.png}
+ \caption{}
+ \label{}
+ \end{figure}
+
+\end{columns}
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Кинетическая энергия системы}
+
+ Подсчитаем кинетическую энергию системы $T$, которая
+ состоит из суммы кинетических энергий первого маятника $T_1$
+ и второго маятника $T_2$.
+
+ \[ T_1 = \frac{Q_1 l_1}{2g} \dot{\varphi_1}^2 \]
+
+ \[ T_2 = \frac{Q_2}{2g} v^2 \]
+ где $v$ -- скорость точки $Q_2$. Вектор скорости $v$ можно рассматривать
+ как сумму вектора вразательной скорости $v'$ точки $Q_2$ вокруг точки
+ $O_2$ и вектора $v_{02}$ скорости точки $O_2$:
+ \[ v = v' + v_{02} \]
+
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{}
+
+ Вектор $v'$ перпендикулярен к $l_2$ и вектор $v_{02}$ перпендикулярен
+ к $O_1 O_2$ $\implies$ из треугольника, образованного векторами
+ $v'$, $v_{02}$ и $v$, видно, что угол между $v'$ и $v_{02}$ будет
+ равен $\pi - (\varphi_1 - \varphi_2)$, откуда
+ \[ v^2 = v'^2 + v_{02}^2 + 2v' v_{02} \cos (\varphi_1 - \varphi_2) \]
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[width=0.6\textwidth] {vectors.png}
+ \caption{}
+ \label{}
+ \end{figure}
+
+ Но $v' = l_1 \dot{\varphi_2}$, $v_{02} = l_2 \dot{\varphi_1}$ и для
+ малых колебаний можно принять $\cos (\varphi_1 - \varphi_2) = 1$, тогда
+ \[ v^2 = l_1^2 \varphi_2^2 + l_1^2 \varphi_1^2 + 2 l_1 l_2 \dot{\varphi_1} \dot{\varphi_2} \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{}
+
+ Из данных вычислений получаем
+ \[ T_2 = \frac{Q_2}{2g} l_1^2 \varphi_2^2 + l_1^2 \varphi_1^2 + 2 l_1 l_2 \dot{\varphi_1} \dot{\varphi_2}\]
+
+ Кинетическую энергию системы теперь запишем в виде:
+ \[
+ T = \frac{Q_1 l_1^2 + Q_2 l_1^2}{2g} \dot{\varphi_1}^2 +
+ \frac{Q_2 l_2^2}{2g} \dot{\varphi_2}^2 +
+ \frac{Q_2 l_1 l_2}{g} \dot{\varphi_1} \dot{\varphi_2}
+ \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Потенциальная энергия системы}
+
+ Так как рассматриваемая система находится под действием сил тяжести,
+ то потенциальная энергия системы будет равна
+ \[ P = Q_1 h_1 + Q_2 h_2 \]
+ где $h_1$ и $h_2$ -- высота точек $m_1$ и $m_2$, под некоторым
+ произвольно выбранным уровнем, который расположен на расстоянии
+ $l_1 + l_2$ от точки $O_1$, тогда
+ \[ h_1 = l_1 + l_2 - l_1 \cos \varphi_1 \]
+ \[ h_2 = l_1 + l_2 - l_1 \cos \varphi_1 - l_2 \cos \varphi_2 \]
+
+ Разлагая косинусы в ряд и ограничиваяст в этих разложениях вторыми
+ степенями малых углов отклонения, получим:
+ \[
+ P = Q_1 \frac{l_1}{2} \varphi_1^2 +
+ Q_2 \left( \frac{l_1}{2} \varphi_1^2 + \frac{l_2}{2} \varphi_2^2 \right)
+ \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Уравнения движения и их интегрирование}
+
+ Зная кинетическую и потенциальную энергию системы, запишем
+ уравнения движения системы (уравнения Лагранжа второго рода)
+ в виде:
+ \[
+ \frac{Q_1 + Q_2}{2g} l_1^2 \varphi_1 +
+ \frac{Q_2 l_1 l_2}{g} \dot{\varphi_2} +
+ (Q_1 + Q_2) l_1 \varphi_1 = 0
+ \]
+
+ \[
+ \frac{Q_2 l_2^2}{g} \ddot{\varphi_2} +
+ \frac{Q_2 l_1 l_2}{g} \varphi_1 +
+ Q_2 l_1 \varphi_2 = 0
+ \]
+
+ Коэффициенты этих уравнения обозначаются в виде
+
+ \[
+ a_{11} = \frac{Q_1 + Q_2}{g} l_1^2, \,
+ a_{12} = \frac{Q_1 l_1 l_2}{g}, \,
+ a_{22} = \frac{Q_2 l_2^2}{g},
+ \]
+ \[
+ c_{11} = (Q_1 + Q_2) l_1, \,
+ c_{22} = Q_2 l_2, \,
+ (c_{12} = 0)
+ \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Уравнение частот}
+
+ Запишем уравнение частот:
+ \[
+ \left( \frac{g}{l_1} - k^2 \right) \left( \frac{g}{l_2} - k^2 \right) -
+ k^4 \frac{Q_2^2}{(Q_1 + Q_2) Q_2} = 0
+ \]
+
+ Введем обозначения:
+ \[
+ \frac{g}{l_1} = n_1^2, \,
+ \frac{g}{l_2} = n_2^2, \,
+ \frac{Q_2^2}{(Q_1 + Q_2) Q_2} = \chi^2
+ \]
+
+ Тогда уравнение частот примет вид:
+ \[ (n_1^2 - k^2) (n_2^2 - k^2) - \chi^2 k^4 = 0 \]
+ \[ (1 - \chi^2) k^2 - (n_1^2 - n_2^2) k^2 + n_1^2 n_2^2 = 0 \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{}
+
+ Корнями уравнения частот являются:
+ \[
+ k_{1,2}^2 = \frac{1}{2 (1 - \chi^2)}
+ \left( n_1^2 + n_2^2 \pm \sqrt{(n_2^2 - n_1^2)^2 + 4 \chi^2 n_1^2 n_2^2} \right)
+ \]
+
+ Выражения, определяющие $k_1^2$ и $k_2^2$, будут положительны.\
+ Значения $k_1$ и $k_2$ определяют собственные частоты системы.
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Уравнения главных колебаний}
+
+ Уравнениями главных колебаний системы в рассматриваемом случае являются:
+
+ \[ \varphi_1^{(1)} = C_1 (c_{11} - k_1^2 a_{11}) \sin (k_1 t + \alpha_1) \]
+ \[ \varphi_2^{(1)} = C_1 k_1^2 a_{12} \sin (k_1 t + \alpha_1) \]
+
+ \[ \varphi_1^{(2)} = C_2 (c_{11} - k_2^2 a_{11}) \sin (k_2 t + \alpha_2) \]
+ \[ \varphi_2^{(2)} = C_2 k_2^2 a_{12} \sin (k_2 t + \alpha_2) \]
+
+ где $C_1$, $C_2$, $\alpha_1$, $\alpha_2$ -- произвольные постоянные,
+ подлежащие определению из начальных условий.
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{}
+
+ Отношение амплитуд главных колебаний будет:
+
+ \[
+ \beta_1 = \frac{a_{12}}{a_{11}} \frac{k_1^2}{\frac{c_{11}}{a_{11}} - k_1^2}, \,
+ \beta_2 = \frac{a_{12}}{a_{11}} \frac{k_2^2}{\frac{c_{11}}{a_{11}} - k_2^2}
+ \]
+
+ Но так как $\frac{c_{11}}{a_{11}} = n_1^2$, то
+
+ \[
+ \beta_1 = \frac{a_{12}}{a_{11}} \frac{k_1^2}{n_1^2 - k_1^2}, \,
+ \beta_2 = \frac{a_{12}}{a_{11}} \frac{k_1^2}{n_1^2 - k_2^2}
+ \]
+
+ Из этих равенств следует, что
+ \[ n_1^2 > k_1^2, \, n_2^2 > k_1^2 \]
+ \[ k_2^2 > n_1^2, \, k_2^2 > n_2^2 \]
+
+ Поэтому \[ \beta_1 > 0, \, \beta_2 < 0 \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{}
+
+\begin{columns}
+
+ \column{0.5\textwidth}
+
+ Имеем, что при главном колебании низшей частоты $k_1$ знаки $\varphi_1$ и
+ $\varphi_2$ одинаковы, а при колебаниях высшей частоты $k_2$ знаки
+ $\varphi_1$ и $\varphi_2$ различны. Это означает, что в первом главном
+ колебании прямые $l_1$ и $l_2$ отклоняются в одну сторону от вертикали и
+ отношение углов отклонения при этом остается постоянным
+ $\varphi_1^{(1)} = \beta_1 \varphi_2^{(1)}$.
+
+ \column{0.5\textwidth}
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[height=0.35\textheight]{two_states.png}
+ \caption{}
+ \label{}
+ \end{figure}
+
+\end{columns}
+
+ Во втором главном колебании отклонение указанных прямых происходит по
+ разные стороны от вертикали и отношение углов также остается неизменным
+ $\varphi_1^2 = \beta_2 \varphi_2^2$.
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{}
+
+ Общее решение исходной системы дифференциальных уравнений будет вида:
+ \[
+ \varphi_1 = C_1 (c_{22} - k_1^2 a_{22}) \sin (k_1 t + \alpha_1) +
+ C_2 (c_{22} - k_2 a_{22}) \sin (k_2 t + \alpha_2)
+ \]
+ \[
+ \varphi_2 = C_1 k_1^2 a_{12} \sin (k_1 t + \alpha_1) +
+ C_2 k_2^2 a_{12} \sin (k_2 t + \alpha_2)
+ \]
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Главные координаты}
+
+ Введем новые переменные:
+ \[
+ \Theta_1 = C_1 \sin (k_1 t + \alpha_1), \,
+ \Theta_2 = C_2 \sin (k_2 t + \alpha_2)
+ \]
+
+ Тогда
+ \[ \varphi_1 = (c_{22} - k_1^2 a_{22}) \Theta_1 + (c_{22} - k_2^2 a_{22}) \Theta_2 \]
+ \[ \varphi_2 = k_1^2 a_{12} \Theta_1 - k_2^2 a_{12} \Theta_2 \]
+
+ Из этого следует, что новые переменные полностью описывают движение
+ рассматриваемой механической системы, следовательно, их можно выбрать в
+ качестве обобщенных координат системы.
+
+ Приведённые ранее выражения для $\beta_1$ и $\beta_2$ указывают, что
+ $\Theta_1$ и $\Theta_2$ являются главными координатами системы.
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Применение}
+
+ Конструкции, похожие на двойной маятник зачастую применяют в местах,
+ где необходимо уменьшить амплитуду механических вибраций. Такие вибрации
+ могут причинять дискомфорт, урон или даже полный отказ системы.
+
+\begin{columns}
+
+ \column{0.6\textwidth}
+
+ Как пример можно привести небоскрёбы, в которых применяются
+ \textbf{инерционные демпферы}. Как правило, демпферы представляют собой
+ огромные бетонные или стальные блоки, установленные в небоскребах или
+ других конструкциях и перемещаемые для компенсации резонансной
+ частоты колебаний конструкции с помощью пружин, жидкости или маятников.
+
+ \column{0.4\textwidth}
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[height=0.35\textheight]{taipei_101.jpg}
+ \caption{}
+ \label{}
+ \end{figure}
+
+
+
+\end{columns}
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Применение}
+
+\begin{columns}
+
+ \column{0.4\textwidth}
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[height=0.4\textheight]{Taipei_101_Tuned_Mass_Damper.png}
+ \caption{Демпфер в небоскребе Тайбэй 101}
+ \label{}
+ \end{figure}
+
+ \column{0.6\textwidth}
+
+ В таких конструкциях в качествет основного \textit{обратного} маятника
+ выступает само здание, а дополнительная масса прикреплена для завершения
+ двойного маятника.
+
+ Когда верх здания приходит в движение, гигантский шар, раскачивается
+ подобно гигантскому маятнику. Он ударяет по масленым амортизаторам,
+ которые рассеивают энергию колебаний. Таким образом, когда здание
+ отклоняется в одну сторону, маятник двигается в другую, сокращая,
+ таким образом, раскачивание небоскрёба.
+
+\end{columns}
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Применение}
+
+\begin{columns}
+
+ \column{0.6\textwidth}
+
+ Похожую конструкцию имеют плавающие маяки, где в качестве главного маяка
+ выступает сам поплавок, а вторым маятником является фонарь.
+
+ \column{0.4\textwidth}
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[height=0.5\textheight]{mayak.png}
+ \caption{}
+ \label{}
+ \end{figure}
+
+\end{columns}
+
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Применение}
+
+ Еще один пример двойного маятника мы имеем в колоколе. В 1876 году в Кёльне
+ имел место случай, на первый взгляд очень странный, -- не удавалось
+ заставить звонить большой колокол, только что подвешенный тогда на башне
+ собора: когда пытались звонить, язык совершал относительно колокола
+ столь малые колебания, что не удавалось произвести удара, хотя язык и был
+ достаточно длинен для того, чтобы достать до стенок колокола. На основе
+ предыдущих рассуждений было установлени, что для этого колокола
+ $l_1 - l = 65.3$см и $\lambda = 66.7$см, так что при ничтожности массы
+ языка по сравнению с массой колокола приближенной равенство их фоз
+ было обнаружено из совпадения этих значений.
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Применение}
+
+\begin{columns}
+
+ \column{0.7\textwidth}
+
+ Двойной маятник подвергается \textbf{хаотическому движению},
+ и его движение очень чувствительно к начальным значениям. На картинке
+ показано количество времени, которое должно пройти перед тем, как маятник
+ перевернется как функция от начальной позиции. Цвет каждого пикселя
+ показывает перевернется ли маятник в течение
+ \begin{itemize}
+ \item $10 \sqrt{l / g}$ секунд (зеленый)
+ \item $100 \sqrt{l / g}$ секунд (красный)
+ \item $1000 \sqrt{l / g}$ секунд (фиолетовый)
+ \item $10000 \sqrt{l / g}$ секунд (синий)
+ \end{itemize}
+
+ Начальные значения, которые не ведут к перевороту в течение $10000 \sqrt{l / g}$
+ нарисованы белым цветом.
+
+ \column{0.3\textwidth}
+
+ \begin{figure}[H]
+ \centering
+ \includegraphics[height=0.4\textheight]{Double_pendulum_flips_graph.png}
+ \caption{}
+ \label{}
+ \end{figure}
+
+\end{columns}
+
+\end{frame}
+
+\begin{frame}
+ \frametitle{Источники}
+
+ \begin{itemize}
+ \item Голубева, О.В. Теоретическая механика / О.В. Голубева. - Москва: ``Высшая школа'', 1968. - 487с.
+ \item Леви-Чивита, T. Курс теоретической механики / Т. Леви-Чивита, У. Амальди - Москва: ``Москва'', 1951. - 556с.
+ \item Аганова, А.Ю. Инерционный демпфер сердце тейбей 101 / А.Ю. Аганова, Н.Д. Комарова // Инновационная наука. - 2015.
+ \item https://en.wikipedia.org/wiki/Double\_pendulum
+ \item https://en.wikipedia.org/wiki/Tuned\_mass\_damper
+ \end{itemize}
+
+\end{frame}
+
+
+
+\end{document} \ No newline at end of file