Перераспределение энергии между поступательными и вращательными степенями свободы

Постановка задачи

Рассмотреть перераспределение энергии между вращательными и поступательными степенями свободы в системе из N тел-точек, соединенных друг с другом балками Бернулли-Эйлера.

Вывод уравнений

Рассматривается система из N тел-точек. Каждое [math]i[/math]-ое тело имеет две степени свободы - смещение вдоль вертикальной оси [math]y_{i}[/math], и угол поворота относительно вертикальной оси [math]\phi_{i}[/math]. Все тела соединены стержнями, которые описываются уравнением балки Бернулли - Эйлера. Движение каждого тела - точки описывается уравнениями:

[math] J \ddot{\phi_{i}} = -M_{i-1}(l) + M_{i}(0) [/math]
[math] m \ddot{y_{i}} = F_{i-1}(l) - F_{i}(0), [/math]
где [math]J - [/math]момент инерции тела-точки. Моменты и силы находим по определению:

[math] M = E \cdot J_{b} \cdot y''(x) [/math]
[math] F = E \cdot J_{b} \cdot y'''(x), [/math]

где [math]E - [/math] модуль юнга материала балки, [math]J_{b} - [/math] момент инерции сечения балки. Вид функции y(x) найдем из уравнения Балки - Бернулли Эйлера:

[math] E \cdot J_{b} \cdot y^{(4)} = 0 [/math]

получаем:

[math] y(x) = c_1 \frac{x^3}{6} + c_2 \frac{x^2}{2} + c_3x + c_4 [/math]
[math] \phi(x) = y'(x) = c_1 \frac{x^2}{2} + c_2x + c_3 [/math]

Для поиска коэффициентов необходимы граничные условия. Для [math]i - [/math]ого тела рассмотрим два участка: балка, соединяющая [math]i - 1[/math] и [math]i[/math] тела:

[math] y(0) = y_{i-1} [/math]
[math] y(l) = y_{i} [/math]
[math] \phi(0) = \phi_{i-1} [/math]
[math] \phi(l) = \phi_{i} [/math]

и на участке, соединяющим [math]i[/math] и[math] i+1 [/math] тела-точки:

[math] y(0) = y_{i} [/math]
[math] y(l) = y_{i+1} [/math]
[math] \phi(0) = \phi_{i} [/math]
[math] \phi(l) = \phi_{i+1} [/math]

где [math]l - [/math]длина балки.

Учитывая граничные условия и все вышеприведенные формулы, находим уравнения движения [math]i - [/math]ого тела: [math] y_{i} = \frac{EJ_{b}}{m}(\frac{12}{l^3}(y_{i+1}-2y_{i}+y_{i-1}) - \frac{6}{l^2}(\phi_{i+1}-\phi_{i-1})) [/math]
[math] \phi_{i} = \frac{EJ_{b}}{J}(\frac{6}{l^2}(y_{i+1}-y_{i}) - \frac{2}{l}(\phi_{i+1}+4\phi_{i}+\phi_{i-1})) [/math]