Dzenushko Dainis. Course project for theoretical mechanics "Double pendulum" — различия между версиями
Dainis (обсуждение | вклад) (→Solution) |
Dainis (обсуждение | вклад) (→Применение метода решения для частного случая) |
||
Строка 74: | Строка 74: | ||
We notice that using this method we can get the motion equation also for big angles, it can be linearized assuming that angles <math>\varphi,\psi</math> are small and dropping terms of second order.<br> | We notice that using this method we can get the motion equation also for big angles, it can be linearized assuming that angles <math>\varphi,\psi</math> are small and dropping terms of second order.<br> | ||
− | == | + | == Using this method for special case == |
− | + | We'll check mentionad method for special case when <math>\alpha = 0</math><br> | |
− | + | In this case problem becomes 2D.<br> | |
− | ''' | + | '''Finding rotation tensors'''<br> |
<math>\underline{\underline{P}}_1(\varphi,\underline{k}) = \underline{k}\underline{k} + (\underline{\underline{E}} - \underline{k}\underline{k})cos(\varphi) + \underline{k} \times \underline{\underline{E}}sin(\varphi)</math><br><br> | <math>\underline{\underline{P}}_1(\varphi,\underline{k}) = \underline{k}\underline{k} + (\underline{\underline{E}} - \underline{k}\underline{k})cos(\varphi) + \underline{k} \times \underline{\underline{E}}sin(\varphi)</math><br><br> | ||
<math>\underline{e}_0 = \cos(\alpha) \underline{k} + \sin(\alpha) \underline{i} = \underline{k}</math><br> | <math>\underline{e}_0 = \cos(\alpha) \underline{k} + \sin(\alpha) \underline{i} = \underline{k}</math><br> | ||
Строка 83: | Строка 83: | ||
<math>\underline{\underline{P}}_2(\psi,\underline{e})= \underline{\underline{P}}_2(\psi,\underline{k}) = \underline{k}\underline{k} + (\underline{\underline{E}} - \underline{k}\underline{k})cos(\psi) + \underline{k} \times \underline{\underline{E}}sin(\psi)</math><br><br> | <math>\underline{\underline{P}}_2(\psi,\underline{e})= \underline{\underline{P}}_2(\psi,\underline{k}) = \underline{k}\underline{k} + (\underline{\underline{E}} - \underline{k}\underline{k})cos(\psi) + \underline{k} \times \underline{\underline{E}}sin(\psi)</math><br><br> | ||
− | ''' | + | '''Finding second rod's angular vielocity'''<br> |
<math>\underline{\omega}_2 = (\dot{\varphi}+\dot{\psi})\underline{k}</math><br><br> | <math>\underline{\omega}_2 = (\dot{\varphi}+\dot{\psi})\underline{k}</math><br><br> | ||
− | ''' | + | '''Finding second rod's center of mass velocity'''<br> |
<math>\upsilon^2_c = \frac{1}{4} b^2 (\dot{\varphi}+\dot{\psi})^2 + ab\cos\psi(\dot{\varphi}+\dot{\psi})\dot{\varphi} + a^2\dot{\varphi}^2</math><br><br> | <math>\upsilon^2_c = \frac{1}{4} b^2 (\dot{\varphi}+\dot{\psi})^2 + ab\cos\psi(\dot{\varphi}+\dot{\psi})\dot{\varphi} + a^2\dot{\varphi}^2</math><br><br> | ||
− | ''' | + | '''Finding second rod's kinetic energy'''<br> |
<math>T_2 = \frac{1}{2} \left( \frac{m_2 b^2}{3}(\dot{\varphi}+\dot{\psi})^2 + m_2 ab\cos\psi(\dot{\varphi}+\dot{\psi})\dot{\varphi} + m_2 a^2\dot{\varphi}^2 \right)</math><br><br> | <math>T_2 = \frac{1}{2} \left( \frac{m_2 b^2}{3}(\dot{\varphi}+\dot{\psi})^2 + m_2 ab\cos\psi(\dot{\varphi}+\dot{\psi})\dot{\varphi} + m_2 a^2\dot{\varphi}^2 \right)</math><br><br> | ||
− | ''' | + | '''Finding second rod's potential energy'''<br> |
<math>\Pi_2 = m_2 g \left[ a \left(1-\cos\varphi \right) + \frac{b}{2}\left(2 + \sin\varphi\sin\psi - \cos\varphi\cos\psi \right) \right]</math><br><br> | <math>\Pi_2 = m_2 g \left[ a \left(1-\cos\varphi \right) + \frac{b}{2}\left(2 + \sin\varphi\sin\psi - \cos\varphi\cos\psi \right) \right]</math><br><br> | ||
− | ''' | + | '''Finding first rod's kinetic and potential energy'''<br> |
<math>T_1 = \frac{1}{2}\frac{m_1 a^2}{3}\dot{\varphi}^2</math><br> | <math>T_1 = \frac{1}{2}\frac{m_1 a^2}{3}\dot{\varphi}^2</math><br> | ||
<math>\Pi_1 = m_1 g \left( \frac{a}{2} - \frac{a}{2} \cos \varphi \right)</math><br><br> | <math>\Pi_1 = m_1 g \left( \frac{a}{2} - \frac{a}{2} \cos \varphi \right)</math><br><br> | ||
− | ''' | + | '''Finding equation of motion for this special case'''<br> |
− | + | Let's write equations for system's kinetic and potential energy:<br> | |
<math>T = \frac{1}{2}\frac{m_1 a^2}{3}\dot{\varphi}^2 + \frac{1}{2} \left( \frac{m_2 b^2}{3}(\dot{\varphi}+\dot{\psi})^2 + m_2 ab\cos\psi(\dot{\varphi}+\dot{\psi})\dot{\varphi} + m_2 a^2\dot{\varphi}^2 \right)</math><br> | <math>T = \frac{1}{2}\frac{m_1 a^2}{3}\dot{\varphi}^2 + \frac{1}{2} \left( \frac{m_2 b^2}{3}(\dot{\varphi}+\dot{\psi})^2 + m_2 ab\cos\psi(\dot{\varphi}+\dot{\psi})\dot{\varphi} + m_2 a^2\dot{\varphi}^2 \right)</math><br> | ||
<math>\Pi = m_1 g \left( \frac{a}{2} - \frac{a}{2} \cos \varphi \right) + m_2 g \left[ a \left(1-\cos\varphi \right) + \frac{b}{2}\left(2 + \sin\varphi\sin\psi - \cos\varphi\cos\psi \right) \right]</math><br> | <math>\Pi = m_1 g \left( \frac{a}{2} - \frac{a}{2} \cos \varphi \right) + m_2 g \left[ a \left(1-\cos\varphi \right) + \frac{b}{2}\left(2 + \sin\varphi\sin\psi - \cos\varphi\cos\psi \right) \right]</math><br> | ||
− | + | Now we differentiate both equations and linearise results assuming that <math>\varphi , \psi</math> are small and taking to account only infinitesimals of the first order. As the result we find equation of motion:<br> | |
<math> | <math> | ||
\begin{cases} | \begin{cases} |
Версия 17:32, 17 сентября 2012
Содержание
Subject
Specification of double pendulum oscillations
Task
Rod is connected to the roof by cylindrical joint. The second rod is connected to the end of the first rod via cylindrical joint so that when both rods are vertical the angle between joint's axis equals to
System parameters:
- Inertia tensor of first and second rods are equal. and respectively.
- length of rods are a and b, their masses and for first and second rods.
- Angle between rotation axis of joints equals to
- - angle between first rod and vertical direction
- - angle between axis of first and second rods that means the angle in second joint relative to vertical orientation of system.
Main task:
- find the equation of motion for this double pendulum
Solution
Let's find the way of solving this problem: We are going to solve the problem by using the Lagrange equation:
Selecting generalized coordinates: for generalized coordinates we are going to use angles
- In this problem we neglect non-potential interactions, so generalized non-potential forces equals to zero.
Calculating potential and kinetic energy of system:
- Potential energy of system
- Kinetic energy of system
- Kinetic energy of first rod; Where
- Inertia torque of first rod
- Potential energy of first rod
- Kinetic energy of second rod
Finding second rod's angular velocity vector:
For calculating we will find rotation tensors of both rods
Where:
- second rod's ongoing rotation axis
- second rod's initial rotation axis
But:
Now using the angular velocity addition formula we get:
So we got the equation:
Finding second rod's center of mass velocity
Finding second rod's kinetic energy
Second rod's inertia tensor:
Now we found all expressions required for second rod's kinetic energy formula:
Finding second rod's potential energy
- radius-vector of second rod's center of mass
Finding the motion equation
We need to differentiate the equations for potential and kinetic energy, the way Lagrange equation requires and use the result in it. As the result we going to get the system of two differential equations which itself is equation of motion.
We notice that using this method we can get the motion equation also for big angles, it can be linearized assuming that angles are small and dropping terms of second order.
Using this method for special case
We'll check mentionad method for special case when
In this case problem becomes 2D.
Finding rotation tensors
Finding second rod's angular vielocity
Finding second rod's center of mass velocity
Finding second rod's kinetic energy
Finding second rod's potential energy
Finding first rod's kinetic and potential energy
Finding equation of motion for this special case
Let's write equations for system's kinetic and potential energy:
Now we differentiate both equations and linearise results assuming that are small and taking to account only infinitesimals of the first order. As the result we find equation of motion:
Обсуждение результатов и выводы
В данной работе был подробно описан алгоритм решения задачи о двойном маятнике в случае когда оба шарнира циллиндрические. Затем данный метод был применен для частного случая плоской задачи.