Рассмотрим консервативную систему с голономными и стационарными связями, имеющую одну степень свободы. В этом случае выражение для кинетической энергии имеет вид
T=12a(q)˙q2 , (63)
где всегда положительный обобщенный коэффициент инерции может зависеть от обобщенной координаты q.
Для получения линейной математической модели разложим a(q) в ряд Маклорена по степеням qи отбросим все слагаемые, начиная со второго порядка малости (т.е. рассмотрим только малые колебания около положения устойчивого равновесия).
Тогда a(q)≈a(0)=a и выражение (63) будет иметь вид T=12a˙q2 .
Рассуждая аналогично, для потенциальной энергии можно получить выражение П=12cq2, где обобщенный коэффициент жесткости c(q)≈c(0)=c.
При таких колебаниях механической системы (линейной изначально либо линеаризованной, как это сделано нами выше) дифференциальное уравнение движения может быть получено, например, способом Лагранжа (49).
Тогда a¨q=−cq; (64.а)
либо ¨q+k2q=0 ; где k=√ca . (64.б)
Дополняя это уравнение начальными условиями t=0;q(0)=q0;˙q(0)=˙q0, будем иметь полную математическую модель для отыскания закона движения q(t).
Виды решения этого дифференциального уравнения общеизвестны:
q=C1eikt+C2e−ikt=Asin(kt)+Bcos(kt)=A0sin(kt+α), (65)
здесь C1;C2, либо A;B, либо A0;α — постоянные интегрирования.
При заданных начальных условиях t=0;q(0)=q0;˙q(0)=˙q0 второй вид (65) будет
q=q0sin(kt)+˙q0kcos(kt), (66)
для третьего вида амплитуда колебаний A0=√q20+(˙q0k)2,
начальная фаза колебаний α=arctgα=q0k˙q0.
Как видно, величина kимеет смысл частоты свободных колебаний (число колебаний за 2π единиц времени); она связана с периодом τ этих колебаний соотношением τ=2πk.
ПРИМЕР 20. Воздушный взрыв дал кораблю начальную скорость в вертикальном направлении V0. Зная площадь его ватерлинии S, водоизмещение D и присоединенную массу жидкости в вертикальном направлении λ, найти для корабля период и амплитуду вертикальных колебаний. Развалом бортов и сопротивлением воды пренебречь.
РЕШЕНИЕ. Обозначим вертикальное перемещение корабля y. Запишем дифференциальное уравнение вертикальных колебаний корабля, спроецировав на вертикаль силы, на него действующие (см.рис.23).
Тогда (Dg+λ)¨y=−Fарх+D=−γSy,
где γ — удельный вес воды.
Приведем уравнение к виду (64.б)
¨y+k2y=0 , здесь k2=γSgD+gλ .
Период вертикальных колебаний будет τ=2πk=2π√D+gλγSg .
Запишем начальные условия: t=0;y0=0;˙y0=V0. Подставив их в (66), получим:
y=V0kcos(kt), где k=√γSgD+gλ .