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

Задача: вывести дифференциальное уравнение собственных затухающих малых колебаний физического маятника и получить его решение.

Ключевые идеи:

• Физический маятник: тело массой m, момент inertia относительно оси вращения I, центр тяжести на расстоянии l от оси подвеса.
• Сtorия момента: при угле отклонения φ действует сила тяжести, создающая момент m g l sin φ, а при силе сопротивления вала/опорного узла возникает затухающий момент пропорциональный скорости φ'.
• Для малых колебаний sin φ ≈ φ, поэтому получаем линейное уравнение с постоянными коэффициентами.

1. Уравнение движения для затухающего физического маятника (малые углы)

• Вращательный баланс: I φ'' + b φ' + m g l sin φ = 0.
• При малых углах: sin φ ≈ φ. Получаем линейное уравнение: I φ'' + b φ' + m g l φ = 0.
• Приведем к стандартному виду: φ'' + (b/I) φ' + (m g l / I) φ = 0.

2. Анализ характера решения

• Обозначим: ω0^2 = m g l / I — естественная частота без затухания, α = b / (2 I) — коэффициент затухания.
• Характеристическое уравнение: r^2 + (b/I) r + (m g l / I) = 0.
• Корни: r_{1,2} = [ - (b/I) ± sqrt( (b/I)^2 - 4 m g l / I ) ] / 2.
• В зависимости от дисперсанта различают три случая.

3. Случаи решения

• Подзатухание (underdamped): (b/I)^2 < 4 m g l / I Пусть ωd = sqrt( m g l / I − (b/(2I))^2 ). Решение: φ(t) = e^{−(b/(2I)) t} [ C1 cos(ωd t) + C2 sin(ωd t) ]. Здесь C1, C2 задают начальные условия φ(0) и φ'(0).

• Критическое затухание: (b/I)^2 = 4 m g l / I Решение: φ(t) = (C1 + C2 t) e^{−(b/(2I)) t}.

• Переподзатухание (overdamped): (b/I)^2 > 4 m g l / I Корни r1, r2 вещественные и отрицательные: r1,2 = [ −(b/I) ± sqrt( (b/I)^2 − 4 m g l / I ) ] / 2. Решение: φ(t) = C1 e^{r1 t} + C2 e^{r2 t}.

4. Частоты и относительная логика

• Естественная беззатуханная частота: ω0 = sqrt( m g l / I ).
• Коэффициент затухания в экспоненте: α = b / (2 I).
• Для слабого затухания часто удобна форма: φ(t) = e^{−α t} [ A cos(ωd t) + B sin(ωd t) ], где ωd = sqrt( ω0^2 − α^2 ).
• Параметры (C1, C2 или A, B) задаются начальными условиями: φ(0) и φ'(0).

5. Примечания по физике и применимости

• I — момент инерции маятника относительно точки подвеса.
• l — расстояние от точки подвеса до центра масс маятника.
• b (или другой обозначенный коэффициент) — коэффициент вязкого затухания в моменте τ = −b φ'.
• Уравнение выводимо для малого угла; при больших углах нужно использовать полный sin φ, и решение будет нелинейным.
• В реальных задачах иногда используется другая формулировка: I φ'' + c φ' + k φ = 0, где k = m g l и c = b, что эквивалентно приведённой записи.

6. Пример использования (концептуальный)

• Пусть маятник имеет известные I, l, m, g и известный коэффициент затухания b.
• Вычисляете ω0 = sqrt(m g l / I) и α = b/(2I).
• Определяете дискриминант D = α^2 − ω0^2 (или напрямую сравниваете b^2 и 4 I m g l).
• В зависимости от условий выбираете подзатухание, критическое или переподзатухание и записываете соответствующую форму φ(t).
• Задаете начальные условия φ(0) и φ'(0) и находите константы.

Итак, итоговое дифференциальное уравнение для собственных затухающих малых колебаний физического маятника: I φ'' + b φ' + m g l φ = 0, которое после нормирования даёт стандартный линейный ОДУ второго порядка с решением, описанным выше, в зависимости от величины затухания.