1.5.1. Свободные колебания.В случае, когда действующая на материальную точку позиционная сила стремиться вернуть ее в исходное положение, движение точки будет носить колебательный характер. Такую силу принято называть восстанавливающей.
Пусть материальная точка может двигаться по известной траектории под действием восстанавливающей силы, структура которой имеет вид ; здесь - проекция силы на касательную к траектории, - известный коэффициент пропорциональности, а - криволинейная координата, являющаяся отклонением точки от начального положения. Тогда дифференциальное уравнение движения материальной точки вдоль траектории будет
.
Перепишем уравнение в виде
, где . (1.29)
Это линейное однородное дифференциальное уравнение второго порядка.
Предположим, что его решение должно иметь вид . (1.30)
Для проверки правильности предположения найдем выражение для и подставим его и (1.30) в (1.29). Вынося общий множитель за скобку, имеем
.
Полученное произведение равно нулю, если . Это означает, что предположение (1.30) справедливо, а решение уравнения (1.29) имеет вид
, где . (1.31)
Если известны, например, начальные условия движения ( - начальное положение, - начальная скорость), для определения постоянных величин и начальные условия следует подставить в выражениях для и . Полученные уравнения
позволяют найти и . (1.32)
Таким образом, под действием восстанавливающей силы материальная точка совершает движение по синусоидальному закону, т.е. гармоническое колебательное движение. Эти колебания называют свободными; наибольшее отклонение - амплитудой, аргумент - фазой, а - начальной фазой колебаний. Величина называется круговой (или циклической) частотой колебаний; она определяет число колебаний, совершаемых точкой за секунд (здесь - период колебаний). Заметим, что частота зависит только от массы точки и коэффициента пропорциональности (т.е. инерционно – жесткостных характеристик), и не зависит от начальных условий движения. В связи с этим частоту свободных колебаний иногда называют собственной частотой механической системы.
Заметим, что если начало отсчета координаты находится не в положении статического равновесия, дифференциальное уравнение движения оказывается неоднородным, а его решение будет суммой из полученного выше решения однородного дифференциального уравнения и частного решения. Физический смысл последнего – расстояние от начала отсчета до положения статического равновесия, около которого и происходят свободные колебания точки.
В частном случае, когда траектория точки – прямая линия, касательная и траектория совпадают; а оказывается полным ускорением точки.
ПРИМЕР 1.7. В результате воздушного взрыва корабль получил вертикальную скорость . Найти период и амплитуду вертикальной качки корабля, если известны площадь его ватерлинии , водоизмещение и присоединенная масса жидкости в вертикальном направлении . Развалом бортов и сопротивлением воды пренебречь.
РЕШЕНИЕ. Обозначим вертикальное перемещение корабля (см.рис.1.11).
Запишем дифференциальное уравнение вертикальных колебаний корабля, спроецировав на вертикаль силы, на него действующие (сила веса и архимедова сила поддержания): ;
где - удельный вес воды, - объем подводной части корабля в равновесном состоянии.
При записи уравнения учтено, что в равновесном положении водоизмещение корабля и сила поддержания равны и противоположно направлены.
Приведем полученное уравнение к виду (1.29)
, здесь .
Тогда период вертикальных колебаний будет .
Система уравнений (1.32) позволяет найти амплитуду и начальную фазу вертикальных колебаний как
.
1.5.2. Свободные колебания при сухом трении.Обсудим влияние сухого трения в ходе решения задачи о движении груза массы по горизонтальной шероховатой поверхности (см. рис.1.12), коэффициент трения скольжения для которой задан; груз присоединен к стене пружиной длины и жесткости .
Нанесем на рисунок силы, действующие на груз при его движении в положительном направлении оси . Запишем дифференциальное уравнение движения груза
и приведем его к виду . (1.33)
Это неоднородное дифференциальное уравнение второго порядка, решение которого есть сумма решений однородного уравнения (1.29) и частного решения. Последнее находится достаточно просто, если предположить, что его величина постоянна и равна . После подстановки в неоднородное уравнение, получим .
Тогда общее решение неоднородного уравнения запишется как .
Заметим, что полученное решение справедливо только при совпадении направлений перемещения и скорости точки (при либо ). При движении материальной точки в противоположном направлении (например, ), сила трения изменит направление на противоположное ( на рис.1.12), что приведет к изменению знака правой части дифференциального уравнения и, как следствие, знака частного решения.
Окончательно имеем:
при ;
при . (1.34)
Зависимость перемещения груза от времени изображена на рис.1.13.
В настоящем пособии ограничимся обсуждением приведенного результата (с подробным решением задачи можно ознакомиться, например, в [1]). Хотя движение, описываемое формулами (1.34), строго говоря, не является периодическим (так как с течением времени максимальные отклонения точки от положения равновесия уменьшаются), в механике такое движение принято называть затухающими колебаниями.
Заметим, что, во-первых, анализируемое движение должно рассматриваться как последовательность этапов, в пределах каждого из которых величина скорости точки сохраняет свой знак; время, положение и скорость в конце предыдущего этапа оказываются начальными условиями для последующего этапа. В результате период затухающих колебаний равен периоду свободных колебаний, а амплитуды колебаний образуют арифметическую прогрессию, т.е. . (1.35)
Во-вторых, движение точки заканчивается в тот момент, когда ее максимальное отклонение (при этом скорость точки равна нулю) оказывается в зоне, называемой зоной нечувствительности. Ширина зоны определяется величиной силы сухого трения (в этой зоне сила упругости пружины оказывается меньше силы сухого трения, т.е. , и точка не может продолжить свое движение).
1.5.3. Свободные колебания при линейно – вязком сопротивлении.Рассмотрим движение материальной точки под действием восстанавливающей силы и силы сопротивления , пропорциональной скорости точки. Составим дифференциальное уравнение движения
и приведем его к виду
, где . (1.36)
Предположим, что его решение должно иметь вид
. (1.37)
Для проверки правильности предположения найдем выражения для и подставим их и (1.37) в (1.36). Вынося общий множитель за скобку, имеем
.
Полученное произведение равно нулю, если . Это означает, что сделанное предположение (1.37) справедливо. Однако характер движения будет зависеть от соотношения величин и .
1. Пусть . Тогда - действительные отрицательные числа и решение будет иметь вид
. (1.38)
Характер движения в этом случае представлен на рис.1.14; движение является апериодическим, эспоненциальным, что говорит о большом сопротивлении.
2. Пусть . Тогда - действительные отрицательные числа и решение будет иметь вид . (1.39)
Характер движения аналогичен представленному на рис.1.14.
3. Пусть . Тогда - комплексные числа и решение может быть приведено к виду
. (1.40)
Характер движения - затухающие колебания, - представлен на рис.1.15.
Период затухающих колебаний
(1.41)
несколько больше периода свободных колебаний, а последовательность амплитуд представляет в этом случае убывающую геометрическую прогрессию, знаменатель которой называется декрементом колебаний. Его величина будет
. (1.42)
В некоторых задачах используется модуль его натурального логарифма; он называется логарифмическим декрементом колебаний .
В тех случаях, когда сила сопротивления пропорциональна скорости в степени выше первой (например, ), дифференциальное уравнение нелинейно, а его решение не может быть записано в виде аналитической зависимости. По результатам численного решения такого уравнения строится график (при большом сопротивлении график похож на рис.1.14; при малом – на рис.1.15, однако интервалы времени между наибольшими отклонениями и отношения соседних отклонений не будут постоянны).
Иногда на практике приходится иметь дело с обратной задачей, когда по записи процесса колебательного движения необходимо сделать вывод о структуре силы сопротивления.
Алгоритм действий следующий: во-первых, измеряются интервалы времени между наибольшими отклонениями. Если интервалы одинаковы, то сила сопротивления является либо силой сухого трения, либо силой линейно – вязкого сопротивления. Во-вторых, измеряются наибольшие отклонения от средней линии; если они образуют арифметическую прогрессию, преобладает сухое трение, а если геометрическую, то линейно-вязкое. В противном случае действует сила сопротивления, пропорциональная скорости точки в степени выше первой.
ПРИМЕР 1.8. (задача 32.64 из [2]). Для определения вязкости жидкости Кулон употреблял следующий метод: подвесив на пружине тонкую пластинку А, он заставлял ее колебаться сначала в воздухе, а затем в той жидкости, вязкость которой надлежало определить, и находил продолжительность одного размаха: - в первом случае и - во втором. Сила трения между пластинкой и жидкостью может быть выражена формулой , где - площадь пластинки, - ее скорость, - коэффициент вязкости. Пренебрегая трением между пластинкой и воздухом, определить коэффициент по найденным из опыта величинам и , если масса пластинки .
РЕШЕНИЕ. Изобразим на рис.1.16 пластинку в сосуде с жидкостью и нанесем силы, на нее действующие.
Если начало координаты поместить в положение статического равновесия, то проекция уравнения движения пластинки на ось будет совпадать с уравнением (1.36).
При колебаниях в воздухе сила будет отсутствовать, а дифференциальное уравнение движения будет совпадать с уравнением (1.29). Эти обстоятельства позволяют воспользоваться полученными выше формулами.
Возведем в квадрат выражение для периода затухающих колебаний и учтем, что частота свободных колебаний может быть выражена через период как . Тогда получим выражение
, где .
После несложных преобразований получим .
1.5.4. Вынужденные колебания материальной точки при действии гармонического возбуждения. В некоторых случаях кроме восстанавливающей силы может дополнительно действовать возмущающая сила (является заданной функцией времени) либо кинематическое возбуждение (движение точки подвеса по заданному закону); последнее, согласно (1.25), может быть учтено добавлением силы инерции переносного движения.
Рассмотрим наиболее распространенный случай, когда возмущающая сила изменяется по гармоническому закону, т.е. .
Дифференциальное уравнение в этом случае будет
.
Оно может быть переписано в виде
, где . (1.43)
Это неоднородное дифференциальное уравнение второго порядка, решение которого есть сумма решений однородного уравнения (1.29) и частного решения. В случае, когда , последнее находится достаточно просто, если предположить, что его вид , где - неизвестная постоянная. Подставим выражений для и в неоднородное уравнение и сравним коэффициенты при синусе в левой и правой частях равенства. Получим, что .
Тогда общее решение неоднородного уравнения запишется как . (1.44)
Как уже говорилось выше, постоянные интегрирования и могут быть определены из начальных условий. Таким образом, искомое движение является суммой гармонических колебаний, происходящих с собственной частотой и частотой вынуждающей силы .
Остановимся подробнее на втором слагаемом, описывающем не зависящие от начальных условий вынужденные колебания (в технических приложениях эти колебания называют вибрацией).
Введем величину , представляющую смещение точки при действии на нее постоянной силы . Величина отношения амплитуды вибрации к этому смещению показывает, во сколько раз динамическое отклонение превосходит ; эта величина называется коэффициентом динамичности . Зависимость модуля коэффициента динамичности от отношения называется амплитудно-частотной характеристикой, а зависимость отставания смещения от вынуждающей силы – фазово-частотной характеристикой. Эти характеристики приведены на рис.1.17.а и 1.17.б соответственно.
Анализ амплитудно-частотной характеристики показывает, что при стремлении частоты вынуждающей силы к нулю, амплитуда вибрации стремится к величине статического отклонения; при бесконечно большой частоте вынуждающей силы амплитуда вибрации стремится к нулю; при совпадении частоты вибрации и собственной частоты амплитуда вибрации может стать бесконечно большой (явление резонанса).
Анализ фазово-частотной характеристики показывает, что смещение и вынуждающая сила до резонанса совпадают по направлению; после резонанса их направления противоположны.
Явление резонанса. Рассмотрим более подробно ситуацию, когда . Частное решение уравнения (1.43) в этом случае следует искать в виде . Подставив это предположение в (1.43), получаем два соотношения
.
Отсюда , . Окончательно получим
. (1.45)
График этой зависимости (при нулевых начальных условиях) приведен на рис.1.18.
Как видно, при происходит неограниченное возрастание амплитуды вибрации, причем рост амплитуды линейно зависит от времени. Из этого факта следует важный вывод: если рабочая частота механизма находится за частотой резонанса, последнюю следует проходить быстро (чтобы амплитуда вибрации не успела достигнуть опасных значений, вызывающих разрушение механизма).
Явление биений. В случае, когда частота близка к собственной частоте , благодаря наложению колебаний наступает своеобразное явление, называемое биением.
Опуская выкладки (при желании с ними можно ознакомиться, например, в [1]), для вибрации получаем зависимость
. (1.46)
График этого движения представлен на рис.1.19.
Показанные здесь биения представляют собой колебания, происходящие с частотой вынуждающей силы, причем амплитуда этих колебаний медленно меняется, следуя также периодическому закону. Чем ближе величина к собственной частоте , тем амплитуда биений больше; при биения переходят в резонанс. После прохода резонанса вновь возникают биения, исчезающие с ростом частоты вынуждающего воздействия.
Кинематическое возбуждение. Как было сказано выше, кинематическое возбуждение колебаний (например, гармонические колебания конца пружины жесткости , на другом конце которой прикреплен груз массы ) эквивалентно действию на груз вынуждающей силы, равной силе инерции переносного движения .
Выражение для коэффициента динамичности для этого случая примет вид
. (1.47)
Амплитудно-частотная и фазово-частотная характеристики движения в этом случае приведены на рис.1.20а и 1.20.б соответственно.
Очевидно, что при вынужденные колебания отсутствуют; с ростом частоты кинематического возбуждения амплитуда вибрации растет, пока не достигнет области резонанса; при дальнейшем возрастании частоты кинематического возбуждения амплитуда вынужденных колебаний стремиться к амплитуде кинематического возбуждения, при этом кинематическое возбуждение и вибрация происходят в противофазе (т.е. ).
ПРИМЕР 1.9. Виброметр используется для определения амплитуд вертикальных колебаний одной из частей машины. В подвижной системе прибора демпфер отсутствует. Собственная частота колебаний виброметра – 6 Гц. В каком диапазоне должна находиться частота вибрации, чтобы погрешность измерения ее амплитуды не превышала 5%?
РЕШЕНИЕ. Очевидно, что в рассматриваемом примере имеет место кинематическое возбуждение колебаний. Из рис.1.20.а и формулы (1.47) видно, что при амплитуда вибрации стремиться к амплитуде кинематического возбуждения.
Очевидно, что в таком случае показания виброметра будут наиболее точно регистрировать колебания объекта, на котором он установлен; при этом кинематическое возбуждение и относительное движение груза будут в противофазе (т.е. ).
Воспользуемся формулой (1.47) для расчета амплитуды вибрации и выразим из нее отношение . Потребуем, чтобы величина этого отношения отличалась от 1 (их точное совпадение) не более, чем на 5%, т.е. была бы равна 1.05.
Тогда на границе искомого диапазона должно выполняться равенство
.
Выполнив расчет, получим, что
.
Из графика 1.20.а видно, что и в дорезонансной области существуют частоты вибрации, где соотношение
выполняется, однако конструировать виброметр для работы в столь узком диапазоне частот не имеет смысла.
Действие полигармонического или периодического возбуждения.В случае, когда вынуждающая сила представляет собой полигармоничекое воздействие (например, ; - целое число), в силу линейности задачи решение может быть получено, как сумма из решений, в каждом из которых рассматривается действие только одного моногармонического воздействия, т.е. .
В случае, когда вынуждающая сила оказывается периодической функцией, можно попробовать разложить ее в тригонометрический ряд Фурье. Если в результате разложения можно ограничиться учетом действия нескольких членов ряда (амплитудные значения остальных членов ряда существенно меньше оставленных, а частоты не совпадают с частой свободных колебаний, т.е. резонансные явления отсутствуют), воздействие, с достаточной степенью точности, полагают полигармоническим. В противном случае можно воспользоваться определением периодического решения в замкнутой форме. В ограниченных рамках настоящего пособия этот материал не приводится, однако ознакомиться с ним можно, например, в [1].
1.5.5. Влияние сопротивления на вынужденные колебания. Рассмотрим движение материальной точки вдоль оси под действием линейной восстанавливающей силы , силы линейно – вязкого сопротивления и гармонической вынуждающей силы . Составим дифференциальное уравнение движения
.
Оно может быть переписано в виде
, (1.48)
где .
Это неоднородное дифференциальное уравнение второго порядка, решение которого есть сумма решений однородного уравнения (1.40) и частного решения. В случае, когда , последнее будем искать в виде , где - неизвестные постоянные. Подставим выражений для и в неоднородное уравнение и сравним коэффициенты при синусе и косинусе в левой и правой частях равенства. Получим систему из двух алгебраических уравнений, решив которую найдем неизвестные величины
. (1.49)
Тогда общее решение неоднородного уравнения запишется как (1.50)
где , а и - постоянные величины, определяемые, например, из начальных условий.
Через некоторый промежуток времени (из-за сопротивления) остается только последнее слагаемое (вибрация); при этом перемещение точки сдвинуто по фазе относительно вынуждающей силы на угол .
На рис.1.20 приведены амплитудно-частотная и фазово-частотная характеристики этого движения.
Отметим несколько особенностей: наличие сопротивления делает амплитуду вибрации при резонансе конечной; с ростом сопротивления величина амплитуды уменьшается (уже при значении явно выраженный максимум отсутствует), а частота, при которой наблюдается максимум, смещается в сторону начала координат; сдвиг фаз между перемещением и вынуждающей силой на каждой частоте вынуждающей силы имеет (при заданном сопротивлении) свою величину.
Вопросы и задачи для самоконтроля
1. Какие модели используются в классической механике для пространства, времени, объектов и их взаимодействия?
2. Что такое структура силы? Приведите примеры наиболее часто встречающихся структур сил.
3. Сформулируйте законы классической механики. В каких системах отсчета они справедливы? Какая система отсчета называется инерциальной?
4. Запишите проекции основного закона классической механики на оси декартовой координатной системы (на оси естественной координатной системы, на оси цилиндрической координатной системы).
5. Достаточно ли знания структуры силы, действующей на точку известной массы, чтобы найти закон ее движения?
6. Можно ли судить о структуре силы, действующей на точку заданной массы, если известны ее траектория и закон движения по ней?
7. Чем отличаются дифференциальные уравнения движения материальной точки в инерциальной и неинерциальной системах отсчета?
8. Запишите условие покоя материальной точки в неинерциальной системе отсчета.
9. Почему частоту свободных колебаний иногда называют собственной частотой?
10. Как выглядит график зависимости координаты (отклонение от положения статического равновесия) от времени при наличии сухого трения? Какую прогрессию образуют амплитуды этих колебаний? Чем обусловлено наличие зоны нечувствительности?
11. Какие типы движения могут иметь место, если кроме восстанавливающей силы на материальную точку действует сила линейно – вязкого сопротивления?
12. Расскажите о явлениях биений и резонансе.
13. Как влияет линейно – вязкое сопротивление на амплитуду вибрации?
14. Что такое кинематическое возбуждение колебаний? В чем отличие амплитудно – частотных характеристик при кинематическом и силовом возбуждениях?
15. Решите следующие задачи из [2] : 26.5; 26.13; 27.5; 27.9; 27.33; 27.18; 27.21; 32.26; 32.62; 32.81; 33.2; 33.14.