bilim.techО проекте

Гармонический осциллятор

Одно уравнение описывает подвеску самосвала, дверной доводчик и вибрацию фундамента. Меняй массу, жесткость и демпфирование и смотри, как колебания затухают.

01 Что это такое

Дифференциальное уравнение задает закон через скорость и ускорение, а его решение дает всю траекторию.

Груз на пружине тянут вниз и отпускают. Ньютон не говорит, где груз будет через секунду. Он говорит, какое у груза ускорение при данном отклонении и скорости: пружина тянет к равновесию, амортизатор тормозит. Уравнение m·x″ + c·x′ + k·x = 0 связывает отклонение x с его производными. Решить его значит найти функцию x(t), которая выполняет этот закон в каждый момент. Подвеска самосвала, дверной доводчик и фундамент насоса описываются одним и тем же уравнением, меняются только три числа.

время t, сотклонение xζ = 0,15: качаетсяζ = 1: быстро и без качанийζ = 3: ползет медленноx₀
Одна пружина и один груз, разное трение. Чем больше коэффициент демпфирования ζ, тем меньше качаний, но слишком большое трение снова замедляет возврат.
m
Масса, кг. Инерция: чем тяжелее груз, тем труднее его разогнать и остановить.
k
Жесткость пружины, Н/м. Сила, с которой пружина возвращает груз при отклонении на 1 м.
c
Демпфирование, Н·с/м. Сила вязкого трения при скорости 1 м/с. Ее дает амортизатор: поршень в масле.
02 Опыт
Оттяни груз и отпусти. Меняй трение и смотри, как гаснут колебания. ТОЧНОЕ РЕШЕНИЕ УРАВНЕНИЯ
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

Как груз на пружине возвращается в равновесие?

Рис. 5.3
ГОТОВЫЕ СЛУЧАИ
Собственная частота ω₀
Коэффициент демпфирования ζ
Затухание до 5% амплитуды

ⓘ График строится по точной формуле решения для каждого из трех режимов, без пошагового счета. Быстрые колебания показаны замедленно, медленные ускоренно: темп указан над графиком.

03 Формулы по буквам

Уравнение, частота и демпфирование

Первая формула это сам закон движения. Вторая говорит, как часто груз качается без трения. Третья решает, будут ли колебания вообще. Под каждой формулой числа из опыта.

Какой закон управляет грузом?

m · x″ + c · x′ + k · x = 0

Это второй закон Ньютона, где все силы перенесены в одну сторону. Каждое слагаемое это сила в ньютонах: инерция груза, вязкое трение в амортизаторе и упругость пружины. В любой момент они уравновешивают друг друга.

m · x″
инерция: масса на ускорение, вторая производная x по времени
c · x′
вязкое трение: растет со скоростью, первой производной x
k · x
пружина: растет с отклонением от равновесия

Как часто груз качается без трения?

ω₀ = √(k / m)

Жесткая пружина возвращает груз быстрее, тяжелый груз разгоняется медленнее. Отсюда собственная частота. Период T = 2π / ω₀ это время одного полного качания, частота f = 1 / T это число качаний в секунду.

ω₀
собственная угловая частота, рад/с
T
период, с: T = 2π / ω₀
f
частота, Гц: f = 1 / T

Будут колебания или нет?

ζ = c2√(k · m)

ζ сравнивает трение с тем, что нужно для возврата без качаний. При ζ < 1 груз качается и затухает: x(t) ≈ x₀ · e−ζω₀t · cos(ωdt), где ωd = ω₀√(1 − ζ²). При ζ = 1 возврат самый быстрый без перехода через ноль. При ζ > 1 груз медленно ползет к равновесию.

ζ
коэффициент демпфирования, безразмерный
ωd
частота затухающих колебаний, чуть ниже ω₀
e−ζω₀t
огибающая: во сколько раз уменьшился размах к моменту t

Откуда берутся уравнение, синусоида и экспонента

На груз действуют две силы. Пружина тянет к равновесию с силой −k·x, амортизатор сопротивляется движению с силой −c·x′. По второму закону Ньютона m·x″ = −k·x − c·x′. Переносим все в левую часть и получаем m·x″ + c·x′ + k·x = 0. Вес груза в уравнение не входит, если x отсчитывать от положения равновесия на сжатой пружине.

Без трения остается x″ = −(k/m)·x: нужна функция, вторая производная которой равна ей самой со знаком минус, умноженной на число. Так ведут себя синус и косинус: (cos ω₀t)″ = −ω₀²·cos ω₀t. Отсюда ω₀² = k/m, и груз качается по синусоиде.

С трением ищем решение в виде x = ert. Каждая производная просто умножает его на r, и уравнение превращается в квадратное: m·r² + c·r + k = 0. Если дискриминант отрицательный (это ζ < 1), корни комплексные: r = −ζω₀ ± iωd. Действительная часть −ζω₀ дает экспоненту-огибающую, мнимая часть ωd дает синусоиду внутри нее. При ζ > 1 оба корня действительные и отрицательные: две убывающие экспоненты без колебаний.

x₀ · e^(−ζω₀t)−x₀ · e^(−ζω₀t)Tx₀t
Синусоида с частотой ωd зажата между экспонентами. За каждый период размах уменьшается в одно и то же число раз.
04 Где и зачем это применяют

Везде, где есть масса, упругость и трение

Допущения модели+

Пружина линейная: сила строго пропорциональна отклонению, как в законе Гука. Трение вязкое: сила пропорциональна скорости, как в масляном амортизаторе. Сухое трение, люфты и упоры в модель не входят. Груз движется по одной оси, внешних сил нет, начальная скорость ноль. Настоящая подвеска сложнее: у нее нелинейные пружины, разное сопротивление на сжатие и отбой и несколько масс, но первую оценку частоты и затухания дает именно это уравнение.

Еще в направлении «Математика»

Все темы направления «Математика» →