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

Задача: вывести дифференциальное уравнение собственных затухающих малых колебаний физического маятника и получить его решение. Ключевые идеи: - Физический маятник: тело массой 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, которое после нормирования даёт стандартный линейный ОДУ второго порядка с решением, описанным выше, в зависимости от величины затухания.