1612134389-38f1c0b45e1f51152856fa5b86345f0b (829489), страница 14
Текст из файла (страница 14)
В данном разделе мы рассмотрим только этот последний случай.Рассмотрим влияние трения на собственные колебания одномерногогармонического осциллятора, описываемого координатой x. Учтем силутрения прямо во втором законе Ньютона. Из опыта известно, что эта силанаправлена противоположно скорости, и при малом абсолютном значениискорости пропорциональна ей.
Это позволяет записать силу трения в видеfтр = −αẋ, где α — положительная константа, характеризующая интенсивность трения. Соответствующее уравнение движенияmẍ = −αẋ − kxразделим на m и введем обозначенияk,ω02 = mλ= α .2mЗдесь ω0 — частота свободных колебаний в отсутствие трения, а величинаλ называется коэффициентом затухания. Уравнение свободных колебанийзапишется в виде(22.1a)ẍ + 2λẋ + ω02 x = 0.Полагая x = Re ert ,получаем для r характеристическое уравнение. Это квадратное уравнение,решения которого+r1,2 = −λ ± λ2 − ω02 .Общее решение есть линейная суперпозиция двух решений:x = Re C1 er1 t + C2 er2 t .(22.2)При сильном трении (λ > ω0 ) оба характеристических показателя отрицательны и общее решение x, убывая без колебаний, стремится к нулю.Такое движение называется апериодическим затуханием.При λ = ω0 движение также оказывается апериодическим:x = (C1 + C2 t) e−λt .Глава III. КОЛЕБАНИЯ92При более слабом трении (λ < ω0 ) характеристические показателикомплексно сопряжены+r1,2 = −λ ± iω,ω = ω02 − λ2 ,и решение имеет видx = ae−λt cos(ωt + ϕ),(22.3)где a и ϕ — вещественные постоянные.
Такое движение представляет собойзатухающие колебания.В важном предельном случае λ ω0 амплитуда колебаний за время T = 2π/ω почти не меняется. Средние за это время квадраты координаты и скорости можно вычислять, считая множитель e−λt постоянным.Эти средние квадраты пропорциональны e−2λt и средняя энергия колебаний убывает по этому же закону:E(t) = E0 e−2λt .Рассмотрим влияние трения на движение осциллятора под действиемвнешней гармонической силы f cos (γt + ϕ).
Уравнение движения запишется в видеf(22.1b)ẍ + 2λẋ + ω02 x = m cos (γt + ϕ).Решение этого неоднородного линейного уравнения представляют собойсумму двух слагаемых — описанного выше свободного движения и вынужденного (установившегося) колебания. Последнее решение удобно искатьв комплексной форме. Для этого запишемcos (γt + ϕ) = Re ei(γt+ϕ) .Отыскивая решение в видеx = Re Bei(γt+ϕ) ,найдем комплексную амплитудуB = beiδ =f.m(ω02 − γ 2 + 2iλγ)Мнимая часть B всегда отрицательна, а действительная (при возрастании γот нуля до бесконечности) меняет знак с положительного на отрицательныйпри γ = ω0 , значит аргумент δ меняется в интервале 0 > δ > −π.§ 22.
Колебания при наличии силы трения93Для вещественной амплитуды b и аргумента δ имеемb=f,2m (ω0 − γ 2 )2 + 4λ2 γ 2ctg δ =γ 2 − ω02.2λγ(22.4)Легко найти, что при заданной амплитуде силы f амплитуда вынужденныхколебаний b максимальна на частоте+(22.5)γm = ω02 − 2λ2и равнаbmax = b(γm ) =f.2mλ ω02 − λ2(22.6)Графики функций b и δ в зависимости от частоты вынуждающей силы γи при различных значениях λ представлены на рис. 32. При малых λ картина колебаний близка к той, что имела место при отсутствии трения (см.рис. 31). А именно, фаза δ близка к нулю (т. е.
вынужденные колебания совершаются почти в фазе с действующей силой) при γ < ω0 и фаза δ близкак (−π) (т. е. вынужденные колебания совершаются почти в противофазес действующей силой) при γ > ω0 , а переход от одного режима к другомупроисходит в узком интервале частот γ вблизи ω0 .Рис. 32. Зависимость амплитуды вынужденных колебаний b и фазы δ от частотывынуждающей силы γ для различных значений коэффициента затухания λСвободные колебания осциллятора с трением затухают, со временемостаются только вынужденные колебания видаx = b cos(γt + ϕ + δ).Глава III.
КОЛЕБАНИЯ94Так как фаза δ отрицательна, аргумент косинуса этого решения имеет значение, которое аргумент косинуса вынуждающей силы имел в более ранниймомент времени; говорят, что вынужденное колебание отстает по фазе отвынуждающей силы.Рассмотрим окрестность резонанса в практически интересном случаеслабого трения λ ω0 . Обозначив расстройку частоты(22.7) = γ − ω0и заменив в (4) приближенноγ 2 − ω02 = (γ + ω0 )(γ − ω0 ) ≈ 2ω0 ,2λγ ≈ 2λω0 ,(22.8)получаемb=f,√2mω0 2 + λ2ctg δ = .λ(22.9)Максимальная амплитуда вынужденных колебаний обратно пропорциональна λ:f.(22.10)bmax =2mω0 λКвадрат амплитуды колебаний при = ±λ в два раза меньше по сравнению с максимальным значением.
Изменение фазы δ также происходит,в основном, в этом же интервале частот. Таким образом, характерная ширина области резонанса равна λ.Полученное решение именно в этом интервале частот существенно отличается от решения для резонанса без трения, даваемого формулой (21.4). Если расстройка частоты больше коэффициента затухания (т. е.если || λ), то влияние трения становится несущественным и решениядля резонанса с трением и без трения практически совпадают.§ 23. Колебания при наличии гироскопических силВ § 19 и 20 мы рассматривали свободные колебания систем, движущихся под действием потенциальных сил.
Линейные колебания, соответствующие малым отклонениям от положения равновесия, в этом случаепроисходят вблизи минимума потенциальной энергии U , а отдельное нормальное колебание (19.24) происходит вдоль линии: x = A cos (ωt + ϕ).Рассмотрим теперь линейные колебания систем, в которых помимо потенциальных сил действуют непотенциальные гироскопические силы Fg , линейные по скорости частицы и ортогональные этой скорости.§ 23.
Колебания при наличии гироскопических сил9523.1. Гироскопические силыМы стартуем от функции Лагранжа для частицыL(r, v) = L0 + Lg ,L0 = m v2 − U (r),2Lg = v · C(r),(23.1)где C(r) — известная функция координат r. Слагаемое L0 соответствуетдвижению частицы в потенциальном поле U (r), а добавка Lg приводитк появлению гироскопической силы Fg . Действительно, уравнения Лагранжаd ∂L ∂L−=0dt ∂v∂rимеют видmr̈ = −∇U (r) + Fg ,Fg = [v, [∇, C(r)]].(23.2)Так как Fg v = 0, работа гироскопической силы равна нулю, а энергияE = m v2 + U (r)2сохраняется.Приведем примеры гироскопических сил.При движении частицы с зарядом e в постоянном магнитном поле B(r)(см.
§ 10) функцияC(r) = ec A(r),где A(r) — векторный потенциал, а гироскопическая сила совпадает с силойЛоренца:eFg = [v, B(r)].cПри движении частицы в системе отсчета, вращающейся с постояннойугловой скоростью Ω (см. § 17.2), функцияC(r) = m [Ω, r],а гироскопическая сила совпадает с кориолисовой силой:Fg = 2m[v, Ω].Линейные колебания при наличии гироскопических сил имеют интересные особенности.
В частности, малые колебания могут происходить нетолько вблизи минимума потенциальной энергии U , но также и вблизи максимума U , а траектория таких колебаний с определенной частотой не обязательно имеет вид прямой (19.24). Рассмотрим три поучительных примера.Глава III. КОЛЕБАНИЯ9623.2. Заряженная частица в потенциальном и магнитном поляхПусть потенциальная энергияU (x, y) =11kx x2 + ky y 2 ,22а постоянное однородное магнитное поле направлено по оси z, т.
е.B = (0, 0, Bz ).Мы будем рассматривать движение частицы лишь в плоскости xy. Потенциальная энергия U (x, y) имеет экстремум в точке x = y = 0. Обозначим2> 0. Уравнения движения частицы имеютωB = eBz /(mc) и kB = mωBвидmẍ + kx x − mωB ẏ = 0,mÿ + ky y + mωB ẋ = 0.Будем искать решения этих уравнений в виде колебанийx = Re(Aeiωt ),y = Re(Beiωt ),где комплексная амплитуда A = aeiϕ . Чтобы система однородных уравнений(kx − mω 2 )A − imωB ωB = 0,imωB ωA + (ky − mω 2 )B = 0имела нетривиальное решение, должно выполняться уравнениеm2 ω 4 − mω 2 (kB + kx + ky ) + kx ky = 0.Легко проверить, что корни этого уравнения+2ω1,2= 1 kB + kx + ky ± (kB + kx + ky )2 − 4kx ky2mположительны при выполнении одного из двух условий:1) либо приkx > 0, ky > 0(23.3)(при этом потенциальная энергия имеет минимум в точке x = y = 0);2) либо при(23.4a)kx < 0, ky < 0§ 23.
Колебания при наличии гироскопических сил97и достаточно большом магнитном полеkB > −kx − ky + 2 kx ky(23.4b)(при этом потенциальная энергия имеет максимум в точке x = y = 0).В обоих этих случаях решения представляют собой колебания вблизи начала координат:x = aα cos(ωα t + ϕα ),y = bα sin(ωα t + ϕα ),гдеbα =mωB ωαaα ,ky − mωα2α = 1, 2.Рассмотрим для определенности вариант kx > ky > 0. В этом случаепервое из найденных колебаний представляет собой движение по часовой стрелке по эллипсу с большой осью, направленной вдоль оси x, а второе — в обратном направлении по эллипсу с большой осью, лежащей вдольоси y (подробнее об этом случае можно прочитать, например, в задаче 6.36из [3]).
Свободное движение осциллятора представляет собой суперпозицию найденных колебаний. Эти колебания можно назвать нормальными,обобщая тем самым понятие нормального колебания: движения в направлениях осей x и y происходят с одной и той же частотой, но со сдвигом фаз.Привести функцию Лагранжа к диагональному виду с помощью линейного преобразования только координат невозможно (переход к нормальнымкоординатам связан в этом случае с каноническим преобразованием — см.задачи 11.7 и 11.9 из [3]).На рис. 33 показана траектория движения частицы для второго варианта.