Решение дифференциального уравнения второго порядка численным методом колебания пружинного маятника

Цель: Используя численные методы решить заданную условием задачу сводящуюся к дифференциальному уравнению либо к их системе. Исследовать фазовое движение маятника в данной задаче (нарисовать фазовую траекторию), построить осциллограмму колебаний. Научится программировать задачи с участием дифференциальных уравнений. Изучить численные методы решения и программирования таких задач.

Анализ условия, решение в общем виде:

Путём подстановки длины нити в уравнение математического маятника получаем следующую систему уравнений:

Выбор численного метода решения и анализ его сходимости:

Данная задача решается известными нам методами: Эйлера и методом Рунге-Кутта. Для начала решим её методом Эйлера. Метод Эйлера характерен своей простотой – это его основное преимущество. Однако у него малая точность (первый порядок точности). В данном методе ошибка «накапливается» потому что приращение значения функции при переходе к следующему шагу заменяется приращением кординаты касательной к кривой решения от предыдущей точки. Метод Эйлера точнее, когда вторая производная функции имеет как можно меньшее значение и «шаги» очень малы.

Алгоритм решения уравнения:

Описание алгоритма: рассмотрев выше составленные уравнения, решаем их по методу Эйлера. Изменяя время на некоторую величину dt, программа выполняет один «шаг» и отправляет все нужные нам значения в графики. Повторяя это, получаем новые значения.

Блок-схема (Для общего случая):

G, l0,a, w,ww, ww0,f, f0,t, dt, f1:extended;

Function func(t, f:extended):extended;

Procedure TForm1.Button1Click(Sender: TObject);

Видео:Механика. Л 10.1. Колебания. Вывод дифференциального уравнения пружинного маятникаСкачать

Вынуяеденные колебания. Резонанс

В реальных колебательных системах всегда имеют место потери энергии, что приводит к затуханию колебаний. Система совершает незатухающие колебания, когда на неё действует внешняя, периодически изменяющаяся сила, компенсирующая потери энергии.

Колебания, возникающие в системе под действием периодической внешней силы, называются вынуяеденными. Переменная внешняя сила, приложенная к колебательной системе и вызывающая вынужденные колебания, называется вынуждающей силой F (О-

В колебательном контуре компенсируются потери энергии, если в контур включён источник электрической энергии, основными характеристиками которого являются электродвижущая сила (ЭДС) и внутреннее сопротивление.

Источник электрической энергии в электротехнике называют источником ЭДС (источником напряжения). В механических системах потери энергии компенсируются работой внешних сил по преодолению сил трения.

Видео:14. Дифференциальные уравнения второго порядка, допускающие понижение порядкаСкачать

Вывод дифференциального уравнения вынуяеденных колебаний пружинного маятника

Вынужденные колебания пружинного маятника происходят под действием вынуждающей силы F _ehm.(t), изменяющейся со временем t по гармоническому закону с циклической частотой со

где F вын.(0 — мгновенное значение (в момент времени г) вынуждающей силы,

F т.вын — амплитуда колебаний вынуждающей силы.

На колебательную систему, кроме вынуждающей силы F _вын, одновременно действуют упругая сила F_vnp и сила трения F_mp, равные

где к — коэффициент упругости пружины, г — коэффициент пропорциональности.

Уравнение движения тела массой пружинного маятника, по второму закону Ньютона, имеет вид

Пусть маятник движется вдоль оси о х (рис. 190). Запишем уравнение (19.81) в проекции на ось ох

разделим уравнение на т — массу тела и введём обозначения

тогда получим дифференциальное уравнение вынужденных колебаний пружинного маятника

Вывод дифференциального уравнения вынужденных колебаний в электрическом колебательном контуре

В реальном электрическом колебательном контуре электрическое сопротивление R отлично от нуля (рис.191). Поэтому, если возникают свободные колебания в контуре, то они постепенно затухают.

Формула закона Ома для участка цепи 1 -R-L-2 имеет вид

где (g>i — q>2) — разность потенциалов обкладок конденсатора,

равная qjj— (р-> = iL,

С q — электрический заряд на обкладках конденсатора,

С — ёмкость конденсатора,

J — сила электрического тока, J = С 1Я_,

Е (t) — электродвижущая сила (ЭДС) внешнего источника напряжения

Запишем уравнение (19.83) в виде

затем разделим на индуктивность катушки L и введём обозначения

где р — коэффициент затухания свободных колебаний,

а>₀ — циклическая частота свободных затухающих колебаний, когда

сопротивление R контура равно нулю (R = 0).

Подставим эти обозначения в (19.84)

Формула (19.85) представляет собой дифференциальное уравнение вынуяеденных колебаний в электрическом колебательном контуре.

Уравнения (19.85) и (19.82) аналогичны друг другу, поэтому заменим их уравнением вида

электрического колебательного контура S = q, 2ft = ,со₀ = 1 .

Уравнение (19.86) является линейным неоднородным (с правой частью, отличной от нуля) дифференциальным уравнением второго порядка. Общее решение этого уравнения представляет собой сумму решения однородного уравнения (S _одн) и частного решения (S _част) неоднородного уравнения

Однородное уравнение для колеблющейся величины S имеет вид

Для пружинного маятника S = х, 2/3 = —,со₀ 2 = —, F₀ =f _т._вЬ1Н, а Д ля

где частота со равна

Уравнение (19.88) описывает собственные затухающие колебания величины S. Они прекращаются через некоторый интервал времени (множитель е убывает со временем по экспоненте, поэтому им можно пренебречь).

Величина S совершает незатухающие (установившиеся во времени) вынужденные колебания, описываемые частным решением уравнения (19.87)

где А — амплитуда колебаний, А = S _тах,

со — циклическая частота колебаний, равная циклической частоте вынуждающей силы,

А ср — сдвиг по фазе между вынуждающей силой и вынужденными колебаниями.

На рис. 192 показан график установившихся вынужденных колебаний величины S. Характер вынужденных колебаний определяется, как величиной внешней силы, так и свойствами колебательной системы. Собственные затухающие колебания имеют место в течение интервала времени t = t _раск , когда происходит раскачка колебательной системы.

Время раскачки t _раск зависит от коэффициента затухания Д чем меньше величина Д тем больше это время. Когда прекращаются собственные затухающие колебания, система совершает вынужденные (незатухающие) установившиеся колебания.

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

Для вывода формул, определяющих амплитуду А установившихся вынужденных колебаний и сдвиг по фазе А C L_§_ в уравнение (19.86)

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

Первое слагаемое на векторной диаграмме (рис. 193) изображается вектором, длина которого равна А со ”. Он направлен в отрицательном направлении горизонтальной оси (противоположно вектору, представляющему третье слагаемое). Третье слагаемое, стоящее в левой части уравнения, представлено на векторной диаграмме вектором, длиной А со₀ 2 , направленным по горизонтальной оси в положительном направлении (вправо).

Второе слагаемое на векторной диаграмме изображается вектором длиной, равной 2А /> со. Этот вектор повёрнут от

горизонтальной оси на угол 71 против часовой стрелки.

На рис. 194,195 представлен результат сложения трёх векторов, когда со со₀.

Векторная сумма векторов с длинами, равными 2А [3 со и А со₀ 2 должна быть равна вектору длиной F₀ (рис. 194).

Результат сложения векторов зависит от соотношения частот со

и С0о. Суммарный вектор длиной F₀ является диагональю

прямоугольника со сторонами 2А /3 со и А (со₀ со ). Он составляет с

горизонтальной осью угол А со_а. В формуле стоит квадрат выражения (со 2 ,, — со 2 ) 2 .

Согласно (19.90), амплитуда А вынужденных колебаний зависит от соотношения собственной частоты со₀ и частоты со вынуждающей силы, а так же от коэффициента затухания /1. Заметим, что увеличение амплитуды колебаний ограничено, независимо от соотношения частот со₀ и со. Это объясняется тем, что с ростом амплитуды увеличивается скорость колебательного движения, а вместе с ней и сила трения (сила сопротивления). Колебательная система интенсивно расходует энергию. Тогда амплитуда автоматически принимает значение, при котором прибыль энергии извне компенсирует расходы на преодоление трения (на излучение и др.).

Сдвиг по фазе А ср, как следует из векторной диаграммы, определяется формулой

Из формул (19.90) и (19.91) следует, что амплитуда А вынужденных колебаний и сдвиг по фазе А ср между вынужденными колебаниями и вынуждающей силой определяются параметрами, характеризующими вынуждающую силу, такими, как амплитуда А, частота со, а так же параметрами колебательной системы, а именно, собственной частотой колебаний со₀, коэффициентом затухания /1.

В механической колебательной системе (пружинный маятник) величины А и А ср зависят от массы т колеблющего тела, а у электрического колебательного контура — индуктивности L катушки. Амплитуда А и сдвиг по фазе А ср установившихся вынужденных колебаний не зависят от времени t и начальных условий.

Вынужденные колебания отстают по фазе от вынуждающей силы, вызывающей их, на величину А ср, зависящую от частоты со вынуждающей силы, которая периодически изменяется во времени. Объясним, почему это происходит на примере механической колебательной системы (пружинного маятника).

Считаем, что в начальный момент времени (t = 0) он покоится. Вынуждающая сила F (?) равна нулю. Затем внешняя сила F (?) смещает тело массы т из положения равновесия, совершая положительную работу.

Через четверть периода (t = Т_) величина внешней силы F (7) и

смещения х достигают максимального значения. Затем, внешняя сила убывает со временем, а тело под действием возвращающей силы (силы упругости F у„р) движется к положению равновесия.

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

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

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

Следует заметить, что параметры, характеризующие вынужденные колебания, происходящие под действием внешней периодической силы, определяются добротностью Q.

Видео:Линейное неоднородное дифференциальное уравнение второго порядка с постоянными коэффициентамиСкачать

Решение дифференциального уравнения второго порядка численным методом колебания пружинного маятника

Простейшими из колебаний являются гармонические. Это колебания, при которых колеблющаяся величина изменяется со временем по закону синуса или косинуса.

Рассмотрим пружинный маятник (Рис. 1.7.1).

Рис. 1.7.1. Пружинный маятник

В состоянии покоя сила тяжести уравновешивается упругой силой:

(1.7.1)

Если сместить шарик от положения равновесия на расстояние х, то удлинение пружины станет равным Δl₀ + х. Тогда результирующая сила примет значение:

(1.7.2)

Учитывая условие равновесия (1.7.1), получим:

(1.7.3)

Знак «минус» показывает, что смещение и сила имеют противоположные направления.

Упругая сила f обладает следующими свойствами:

1. Она пропорциональна смещению шарика из положения равновесия;
2. Она всегда направлена к положению равновесия.

Для того, чтобы сообщить системе смещение х, нужно совершить против упругой силы работу:

(1.7.4)

Эта работа идет на создание запаса потенциальной энергии системы:

(1.7.5)

Под действием упругой силы шарик будет двигаться к положению равновесия со все возрастающей скоростью . Поэтому потенциальная энергия системы будет убывать, зато возрастает кинетическая энергия (массой пружины пренебрегаем). Придя в положение равновесия, шарик будет продолжать двигаться по инерции. Это — замедленное движение и прекратится тогда, когда кинетическая энергия полностью перейдет в потенциальную. Затем такой же процесс будет протекать при движении шарика в обратном направлении. Если трение в системе отсутствует, шарик будет колебаться неограниченно долго.

Уравнение второго закона Ньютона в этом случае имеет вид:

(1.7.6)

Преобразуем уравнение так:

(1.7.7)

Вводя обозначение , получим линейное однородное дифференциальное уравнение второго порядка:

(1.7.8)

Прямой подстановкой легко убедиться, что общее решение уравнения (1.7.8) имеет вид:

(1.7.9)

где а — амплитуда и φ — начальная фаза колебания — постоянные величины. Следовательно, колебание пружинного маятника является гармоническим (Рис. 1.7.2).

Рис. 1.7.2. Гармоническое колебание

Вследствие периодичности косинуса различные состояния колебательной системы повторяются через определенный промежуток времени (период колебаний) Т, за который фаза колебания получает приращение 2π. Рассчитать период можно с помощью равенства:

(1.7.10)

(1.7.11)

Число колебаний в единицу времени называется частотой:

(1.7.12)

За единицу частоты принимается частота такого колебания, период которого равен 1 с. Такую единицу называют 1 Гц.

Из (1.7.11) следует, что:

(1.7.13)

Следовательно, ω₀ — это число колебаний, совершаемое за 2π секунд. Величину ω₀ называют круговой или циклической частотой. Используя (1.7.12) и (1.7.13), запишем:

(1.7.14)

Дифференцируя (1.7.9) по времени, получим выражение для скорости шарика:

(1.7.15)

Из (1.7.15) следует, что скорость также изменяется по гармоническому закону и опережает смещение по фазе на ½π. Дифференцируя (1.7.15), получим ускорение:

(1.7.16)

1.7.2. Математический маятник

Математическим маятником называют идеализированную систему, состоящую из нерастяжимой невесомой нити, на которой подвешено тело, вся масса которого сосредоточена в одной точке.

Отклонение маятника от положения равновесия характеризуют углом φ, образованным нитью с вертикалью (Рис. 1.7.3).

Рис. 1.7.3. Математический маятник

При отклонении маятника от положения равновесия возникает вращательный момент, который стремится вернуть маятник в положение равновесия:

(1.7.17)

Напишем для маятника уравнение динамики вращательного движения, учитывая, что момент его инерции равен ml 2 :

(1.7.18)

Это уравнение можно привести к виду:

(1.7.19)

Ограничиваясь случаем малых колебаний sinφ ≈ φ и вводя обозначение:

(1.7.20)

уравнение (1.7.19) может быть представлено так:

(1.7.21)

что совпадает по форме с уравнением колебаний пружинного маятника. Следовательно, его решением будет гармоническое колебание:

(1.7.22)

Из (1.7.20) следует, что циклическая частота колебаний математического маятника зависит от его длины и ускорения свободного падения. Используя формулу для периода колебаний (1.7.11) и (1.7.20), получим известное соотношение:

(1.7.23)

1.7.3. Физический маятник

Физическим маятником называется твердое тело, способное совершать колебания вокруг неподвижной точки, не совпадающей с центром инерции. В положении равновесия центр инерции маятника С находится под точкой подвеса О на одной с ней вертикали (Рис. 1.7.4).

Рис. 1.7.4. Физический маятник

При отклонении маятника от положения равновесия на угол φ возникает вращательный момент, который стремится вернуть маятник в положение равновесия:

(1.7.24)

где m — масса маятника, l — расстояние между точкой подвеса и центром инерции маятника.

Напишем для маятника уравнение динамики вращательного движения, учитывая, что момент его инерции равен I:

(1.7.25)

Для малых колебаний sinφ ≈ φ. Тогда, вводя обозначение:

(1.7.26)

(1.7.27)

что также совпадает по форме с уравнением колебаний пружинного маятника. Из уравнений (1.7.27) и (1.7.26) следует, что при малых отклонениях физического маятника от положения равновесия он совершает гармоническое колебание, частота которого зависит от массы маятника, момента инерции и расстояния между осью вращения и центром инерции. С помощью (1.7.26) можно вычислить период колебаний:

(1.7.28)

Сравнивая формулы (1.7.28) и (1.7.23) получим, что математический маятник с длиной:

(1.7.29)

будет иметь такой же период колебаний, что и рассмотренный физический маятник. Величину (1.7.29) называют приведенной длиной физического маятника. Следовательно, приведенная длина физического маятника — это длина такого математического маятника, период колебаний которого равен периоду колебаний данного физического маятника.

Точка на прямой, соединяющей точку подвеса с центром инерции, лежащая на расстоянии приведенной длины от оси вращения, называется центром качания физического маятника. По теореме Штайнера момент инерции физического маятника равен:

(1.7.30)

где I₀ — момент инерции относительно центра инерции. Подставляя (1.7.30) в (1.7.29), получим:

(1.7.31)

Следовательно, приведенная длина всегда больше расстояния между точкой подвеса и центром инерции маятника, так что точка подвеса и центр качания лежат по разные стороны от центра инерции.

1.7.4. Энергия гармонических колебаний

При гармоническом колебании происходит периодическое взаимное превращение кинетической энергии колеблющегося тела Е_к и потенциальной энергии Е_п, обусловленной действием квазиупругой силы. Из этих энергий слагается полная энергия Е колебательной системы:

(1.7.32)

Распишем последнее выражение

(1.7.33)

Но к = mω 2 , поэтому получим выражение для полной энергии колеблющегося тела

(1.7.34)

Таким образом полная энергия гармонического колебания постоянна и пропорциональна квадрату амплитуды и квадрату круговой частоты колебания.

1.7.5. Затухающие колебания .

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

Если в системе кроме квазиупругой силы действуют силы сопротивления среды (силы трения), то второй закон Ньютона можно записать так:

.	(1.7.34.а)

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

,	(1.7.34.б)

где r – коэффициент трения, характеризующий свойства среды оказывать сопротивление движению. Подставим (1.7.34б) в (1.7.34а):

,	(1.7.34.в)

где β – коэффициент затухания; ω 0 – круговая частота собственных колебаний системы.

Решение уравнения(1.7.34.в) существенно зависит от знака разности: , где ω – круговая частота затухающих колебаний. При круговая частота ω является действительной величиной и решение (1.7.34.в) будет следующим:

.

(1.7.35)

График этой функции показан на рис.1.7.5 сплошной кривой 1, а штриховой линией 2 изображено изменение амплитуды:

.	(1.7.35.а)

Период затухающих колебаний зависит от коэффициента трения и определяется формулой

.	(1.7.35.б)

При очень малом трении период затухающего колебания близок к периоду незатухающего свободного колебания (1.7.35.б)

Рис.1.7.5. Затухающее колебание	Рис.1.7.6. Апериодический процесс

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

;

Следовательно, коэффициент затухания и логарифмический декремент затухания связаны достаточно простой зависимостью:

При сильном затухании из формулы (1.7.37) видно, что период колебания является мнимой величиной. Движение в этом случае уже называется апериодическим . График апериодического движения в виде показан на рис. 1.7.6. Незатухающие и затухающие колебания называют собственными или свободными . Они возникают вследствие начального смещения или начальной скорости и совершаются при отсутствии внешнего воздействия за счет первоначально накопленной энергии.

1.7.6. Вынужденные колебания. Резонанс .

Вынужденными колебаниями называются такие, которые возникают в системе при участии внешней силы, изменяющейся по периодическому закону.

Предположим, что на материальную точку кроме квазиупругой силы и силы трения действует внешняя вынуждающая сила

,

где F 0 – амплитуда; ω – круговая частота колебаний вынуждающей силы. Составим дифференциальное уравнение (второй закон Ньютона):

,

,

(1.7.38)

где .

Решение дифференциального уравнения (3.19) является суммой двух колебаний: затухающих и незатухающих с амплитудой

,

(1.7.39)

Амплитуда вынужденного колебания (1.7.39) прямо пропорциональна амплитуде вынуждающей силы и имеет сложную зависимость от коэффициента затухания среды и круговых частот собственного и вынужденного колебания. Если ω 0 и β для системы заданы, то амплитуда вынужденных колебаний имеет максимальное значение при некоторой определенной частоте вынуждающей силы, называемой резонансной .

Само явление – достижение максимальной амплитуды для заданных ω 0 и β – называют резонансом.

Рис. 1.7.7. Резонанс

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

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

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

1.7.7. Автоколебания

Существуют и такие колебательные системы, которые сами регулируют периодическое восполнение растраченной энергии и поэтому могут колебаться длительное время.

Незатухающие колебания, существующие в какой-либо системе при отсутствии переменного внешнего воздействия, называются автоколебаниями , а сами системы – автоколебательными.

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

Рис. 1.7.8. Блок-схема автоколебаний

Во многих случаях автоколебательные системы можно представить тремя основными элементами (рис.1.7.8): 1) собственно колебательная система; 2) источник энергии; 3) регулятор поступления энергии в собственно колебательную систему. Колебательная система каналом обратной связи (рис. 6) воздействует на регулятор, информирую регулятор о состоянии этой системы.

Классическим примером механической автоколебательной системы являются часы, в которых маятник или баланс являются колебательной системой, пружина или поднятая гиря – источником энергии, а анкер – регулятором поступления энергии от источника в колебательную систему.

Многие биологические системы (сердце, легкие и др.) являются автоколебательными. Характерный пример электромагнитной автоколебательной системы – генераторы автоколебательных колебаний.

1.7.8. Сложение колебаний одного направления

Рассмотрим сложение двух гармонических колебаний одинакового направления и одинаковой частоты:

x 1 =a 1 cos(ω 0 t + α 1 ), x 2 =a 2 cos(ω 0 t + α 2 ).

Гармоническое колебание можно задать с помощью вектора, длина которого равна амплитуде колебаний, а направление образует с некоторой осью угол, равный начальной фазе колебаний. Если этот вектор вращается с угловой скоростью ω 0 , то его проекция на выбранную ось будет изменяться по гармоническому закону. Исходя из этого, выберем некоторую ось Х и представим колебания с помощью векторов а 1 и а 2 (рис.1.7.9).

Рис.1.7.9

Вектор а является суммой векторов а 1 и а 2 . Проекция вектора а на ось Х равна сумме проекций векторов а 1 и а 2 :

Следовательно, вектор а представляет собой результирующее колебание. Этот вектор вращается с той же угловой скоростью, что и векторы а 1 и а 2 . Таким образом, результирующее движение представляет собой гармоническое колебание с частотой ω 0 , амплитудой а и начальной фазой α. Используя теорему косинусов, находим значение амплитуды результирующего колебания:

(1.7.40)

Из рис.1.7.6 следует, что

.

Схемы, в которых колебания изображаются графически в виде векторов на плоскости, называются векторными диаграммами.

Из формулы 1.7.40 следует. Что если разность фаз обоих колебаний равна нулю, амплитуда результирующего колебания равна сумме амплитуд складываемых колебаний. Если разность фаз складываемых колебаний равна , то амплитуда результирующего колебания равна . Если частоты складываемых колебаний не одинаковы, то векторы, соответствующие этим колебаниям будут вращаться с разной скоростью. В этом случае результирующий вектор пульсирует по величине и вращается с непостоянной скоростью. Следовательно, в результате сложения получается не гармоническое колебание, а сложный колебательный процесс.

1.7.9. Биения

Рассмотрим сложение двух гармонических колебаний одинакового направления мало отличающихся по частоте. Пусть частота одного из них равна ω , а второго ω+∆ω, причем ∆ω 1 =a cos ωt, x 2 =a cos(ω+∆ω)t.

Сложив эти выражения и используя формулу для суммы косинусов, получаем:

(1.7.41)

(во втором множителе пренебрегаем членом по сравнению с ω). График функции (1.7.41) изображен на рис. 1.7.10.

Рис.1.7.10

Колебания (1.7.41) можно рассматривать как гармоническое колебание частотой ω, амплитуда которого изменяется по закону . Эта функция является периодической с частотой в два раза превышающей частоту выражения, стоящего под знаком модуля, т.е. с частотой ∆ω. Таким образом, частота пульсаций амплитуды, называемая частотой биений, равна разности частот складываемых колебаний.

1.7.10. Сложение взаимно перпендикулярных колебаний (фигуры Лиссажу)

Если материальная точка совершает колебания как вдоль оси х, так и вдоль оси у, то она будет двигаться по некоторой криволинейной траектории. Пусть частота колебаний одинакова и начальная фаза первого колебания равна нулю, тогда уравнения колебаний запишем в виде:

где α – разность фаз обоих колебаний.

Выражение (1.7.42) представляет заданное в параметрическом виде уравнение траектории, по которой движется точка, участвующая в обоих колебаниях. Если исключить из уравнений (1.7.42) параметр t, то получим уравнение траектории в обычном виде:

(1.7.43)

Уравнение (1.7.43) представляет собой уравнение эллипса, оси которого ориентированы произвольно относительно координатных осей х и у. Ориентация эллипса и величина его полуосей зависят от амплитуд а и b и разности фаз α. Рассмотрим некоторые частные случаи:

α=mπ (m=0, ±1, ±2, …). В этом случае эллипс вырождается в отрезок прямой

,

(1.7.44)

где знак плюс соответствует нулю и четным значениям m (рис 1.7.8.а), а знак минус – нечетным значениям m (рис.1.7.8.б). Результирующее колебание является гармоническим с частотой ω, амплитудой , совершающимся вдоль прямой (1.7.44), составляющей с осью х угол (рис.1.7.11).

Рис.1.7.11.а

Рис.1.7.11. б

α=(2m+1)

(m=0, ±1, ±2, …). В этом случае уравнение имеет вид

Это уравнение эллипса, оси которого совпадают с осями координат, а его полуоси равны амплитудам (рис. 1.7.12). Если амплитуды равны, то эллипс становится окружностью.

Рис.1.7.12

Если частоты взаимно перпендикулярных колебаний отличаются на малую величину ∆ω, их можно рассматривать как колебания одинаковой частоты, но с медленно изменяющейся разностью фаз. В этом случае уравнения колебаний можно записать

x=a cos ωt, y=b cos[ωt+(∆ωt+α)]

и выражение ∆ωt+α рассматривать как разность фаз, медленно изменяющуюся со временем по линейному закону. Результирующее движение в этом случае происходит по медленно изменяющейся кривой, которая будет последовательно принимать форму, отвечающую всем значениям разности фаз от -π до+π.

Если частоты взаимно перпендикулярных колебаний не одинаковы, то траектория результирующего движения имеет вид довольно сложных кривых, называемых фигурами Лиссажу . Пусть, например, частоты складываемых колебаний относятся как 1 : 2 и разность фаз π/2. Тогда уравнения колебаний имеют вид

x=a cos ωt, y=b cos[2ωt+π/2].

За то время, пока вдоль оси х точка успевает переместиться из одного крайнего положения в другое, вдоль оси у, выйдя из нулевого положения, она успевает достигнуть одного крайнего положения, затем другого и вернуться. Вид кривой показан на рис. 1.7.13. Кривая при таком же соотношении частот, но разности фаз равной нулю показана на рис.1.7.14. Отношение частот складываемых колебаний обратно отношению числа точек пересечения фигур Лиссажу с прямыми, параллельными осям координат. Следовательно, по виду фигур Лиссажу можно определить соотношение частот складываемых колебаний или неизвестную частоту. Если одна из частот известна.

Рис.1.7.13

Рис.1.7.14

Чем ближе к единице рациональная дробь, выражающая отношение частот колебаний, тем сложнее получающиеся фигуры Лиссажу.

1.7.11. Распространение волн в упругой среде

Если в каком-либо месте упругой (твёрдой жидкой или газообразной) среды возбудить колебания её частиц, то вследствие взаимодействия между частицами это колебание будет распространяться в среде от частицы к частице с некоторой скоростью υ. процесс распространения колебаний в пространстве называется волной .

Частицы среды, в которой распространяется волна, не вовлекаются волной в поступательное движение, они лишь совершают колебания около своих положений равновесия.

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

На рис. 1.7.12 показано движение частиц при распространении в среде поперечной волны. Номерами 1,2 и т. д. обозначены частицы отстающие друг от друга на расстояние, равное (¼ υT), т.е. на расстояние, проходимое волной за четверть периода колебаний, совершаемых частицами. В момент, времени принятый за нулевой, волна, распространяясь вдоль оси слева направо, достигла частицы 1, вследствие чего частица начала смещаться из положения равновесия вверх, увлекая за собой следующие частицы. Спустя четверть периода частица 1 достигает крайнего верхнего положения равновесия частица 2. По пришествие ещё четверти периода первая часть будет проходить положение равновесия, двигаясь в направлении сверху вниз, вторая частица достигнет крайнего верхнего положения, а третья частица начнёт смещаться вверх из положения равновесия. В момент времени равный T, первая частица закончит полный цикл колебания и будет находиться в таком же состоянии движения, как чальный момент. Волна к моменту времени T, пройдя путь (υT), достигнет частицы 5.

На Рис. 1.7.13 показано движение частиц при распространении в среде продольной волны. Все рассуждения, касающиеся поведения частиц в поперечной волне, могут быть отнесены и к данному случаю с заменой смещений вверх и вниз смещениями вправо и влево.

Из рисунка видно, что при распространении продольной волны в среде создаются чередующиеся сгущения и разряжения частиц (места сгущения обведены на рисунке пунктиром), перемещающиеся в направлении распространения волны со скоростью υ.

Рис. 1.7.15

Рис. 1.7.16

На рис. 1.7.15 и 1.7.16 показаны колебания частиц, положения, равновесия которых лежат на оси x. В действительности колеблются не только частицы, расположенные вдоль оси x, а совокупность частиц, заключённых в некотором объёме. Распространяясь от источников колебаний, волновой процесс охватывает всё новые и новые части пространства, геометрическое место точек, до которых доходят колебания к моменту времени t, называется фронтом волны (или волновым фронтом). Фронт волны представляет собой ту поверхность, которая отделяет часть пространства, уже вовлеченную в волновой процесс, от области, в которой колебания ещё не возникли.

Геометрическое место точек, колеблющихся в одинаковой фазе, называется волновой поверхностью . Волновую поверхность можно провести через любую точку пространства, охваченного волновым процессом. Следовательно, волновых поверхностей существует бесконечное множество, в то время как волновой фронт каждый момент времени только один. Волновые поверхности остаются не подвижными (они проходят через положения равновесия частиц, колеблющихся в одной фазе ). Волновой фронт всё время перемещается.

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

Рис. 1.7.17

Пусть плоская волна распространяется вдоль оси x . Тогда все точки сферы, положения, равновесия которых имеет одинаковую координату x (но различие значения координат y и z), колеблются в одинаковой фазе.

На Рис. 1.7.17 изображена кривая, которая даёт смещение ξ из положения равновесия точек с различными x в некоторый момент времени. Не следует воспринимать этот рисунок как зримое изображение волны. На рисунке показан график функций ξ ( x, t) для некоторого фиксированного момента времени t. Такой график можно строить как для продольной так и для поперечной волны.

Расстояние λ, на короткое распространяется волна за время, равное периоду колебаний частиц среды, называется длиной волны . Очевидно, что

где υ – скорость волны, T – период колебаний. Длину волны можно определить также как расстояние между ближайшими точками среды, колеблющимися с разностью фаз, равной 2π (см. рис. 1.7.14)

Заменив в соотношении(1.7.45) T через 1/ν (ν – частота колебаний), получим

К этой формуле можно придти также из следующих соображений. За одну секунду источник волн совершает ν колебаний, порождая в среде при каждом колебании один «гребень» и одну «впадину» волны. К тому моменту, когда источник будет завершать ν — е колебание, первый «гребень» успеет пройти путь υ. Следовательно, ν «гребней» и «впадин» волны должны уложиться в длине υ.

1.7.12. Уравнение плоской волны

Уравнением волны называется выражение, которое дает смещение колеблющейся частицы как функцию ее координат x, y, z и времени t :

(имеются в виду координаты равновесного положения частицы). Эта функция должна быть периодической относительно времени t , и относительно координат x, y, z. . Периодичность по времени вытекает из того, что точки, отстоящие друг от друга на расстоянии λ , колеблются одинаковым образом.

Найдем вид функции ξ в случае плоской волны, предполагая, что колебания носят гармонический характер. Для упрощения направим оси координат так, чтобы ось x совпадала с направлением распространения волны. Тогда волновые поверхности будут перпендикулярными к оси x и, поскольку все точки волновой поверхности колеблются одинаково, смещение ξ будет зависеть только от x и t :

Рис.1.7.18

Пусть колебания точек, лежащих в плоскости x = 0 (рис. 1.7.18), имеют вид

Найдем вид колебания точек в плоскости, соответствующей произвольному значению x . Для того, чтобы пройти путь от плоскости x =0 до этой плоскости, волне требуется время ( υ – cкорость распространения волны). Следовательно, колебания частиц, лежащих в плоскости x , будут отставать по времени на τ от колебаний частиц в плоскости x = 0 , т.е. будут иметь вид

Итак, уравнение плоской волны (продольной, и поперечной), распространяющейся в направлении оси x , выглядит следующим образом:

(1.7.47)

Величина а представляет собой амплитуду волны. Начальная фаза волны α определяется выбором начала отсчета x и t . При рассмотрении одной волны начало отсчета времени и координаты обычно выбирают так, чтобы α была равной нулю. При совместном рассмотрении нескольких волн сделать так, чтобы для всех них начальные фазы равнялись нулю, как правило, не удается.

Зафиксируем какое – либо значение фазы, стоящей в уравнении (1.7.47), положив

(1.7.48)

Это выражение определяет связь между временем t и тем местом x , в котором фаза имеет зафиксированное значение. Вытекающее из него значение dx/dt дает скорость, с которой перемещается данное значение фазы. Продифференцировав выражение (1.7.48), получим

.

(1.7.49)

Таким образом, скорость распространения волны υ уравнении (1.7.47) есть скорость перемещения фазы, в связи с чем, ее называют фазовой скоростью.

Согласно (1.7.49) dx/dt> 0, следовательно, уравнение (1.7.47) описывает волну, распространяющуюся в сторону возрастания x .

Волна, распространяющаяся в противоположном направлении, описывается уравнением

(1.7.50)

Действительно, приравняв константе фазу волны (1.7.50) и продифференцировав получившееся равенство, придем к соотношению

,

из которого следует, что волна (1.7.50) распространяется в сторону убывания x .

Уравнению плоской волны можно придать симметричный относительно x и t вид. Для этого введем величину

,

(1.7.51)

которая называется волновым числом. Умножив числитель и знаменатель последнего выражения на частоту ν, и вспомнив, что , можно представить волновое число в виде

.

(1.7.52)

Раскрыв в уравнении волны

круглые скобки и используя волновое число, придем к следующему уравнению плоской волны, распространяющейся вдоль оси :

(1.7.53)

Уравнение волны, распространяющейся в сторону убывания x :

При выводе формулы (1.7.53) мы предполагали, что амплитуда колебаний не зависит от x . Для плоской волны это наблюдается в том случае, когда энергия волны не поглощается средой. При распространении в поглощающей энергию среде интенсивность волны с удалением от источника колебаний постепенно уменьшается – наблюдается затухание волны. Опыт показывает, что в однородной среде такое затухание происходит по экспоненциальному закону:

Соответственно уравнение плоской волны, с учетом затухания , имеет следующий вид:

(1.7.54)

(a 0 – амплитуда в точках плоскости x = 0).

© ФГОУ ВПО Красноярский государственный аграрный университет, 2013

Видео:15. Линейные однородные дифференциальные уравнения второго порядка с постоянными коэффициентамиСкачать

Свободные колебания пружинного маятника. Общие сведения

Цель работы. Ознакомиться с основными характеристиками незатухающих и затухающих свободных механических колебаний.

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

Приборы и принадлежности. Штатив со шкалой, пружина, набор грузов различной массы, сосуд с водой, секундомер.

1. Свободные колебания пружинного маятника. Общие сведения

Колебаниями называются процессы, в которых периодически изменяется одна или несколько физических величин, описывающих эти процессы. Колебания могут быть описаны различными периодическими функциями времени. Простейшими колебаниями являются гармонические колебания – такие колебания, при которых колеблющаяся величина (например, смещение груза на пружине) изменяется со временем по закону косинуса или синуса. Колебания, возникающие после действия на систему внешней кратковременной силы, называются свободными.

Рассмотрим одну из простейших колебательных систем – пружинный маятник, представляющий собой груз массой m, подвешенный на абсолютно упругой пружине с коэффициентом жесткости k
(рис. 1). Пусть l0 – длина пружины без подвешенного к ней груза. При подвешивании груза под действием силы тяжести пружина растянется на x1 так, что маятник будет находиться в положении равновесия вследствие равенства модулей силы тяжести mg и упругой силы Fупр: mg = kx1, стремящейся вернуть груз в положение равновесия (полагается, что деформации пружины идеально упругие и подчиняются закону Гука).

Если груз вывести из положения равновесия, отклонив на величину x, то сила упругости возрастает: Fупр = – kx2= – k(x1 + x). Дойдя до положения равновесия, груз будет обладать отличной от нуля скоростью и пройдет положение равновесия по инерции. По мере дальнейшего движения будет увеличиваться отклонение от положения равновесия, что приведет к возрастанию силы упругости, и процесс повторится в обратном направлении. Таким образом, колебательное движение системы обусловлено двумя причинами: 1) стремлением тела вернуться в положении равновесия и 2) инерцией, не позволяющей телу мгновенно остановиться в положении равновесия. В отсутствии сил трения колебания продолжались бы сколь угодно долго. Наличие силы трения приводит к тому, что часть энергии колебаний переходит во внутреннюю энергию и колебания постепенно затухают. Такие колебания называются затухающими.

Незатухающие свободные колебания

Сначала рассмотрим колебания пружинного маятника, на который не действуют силы трения – незатухающие свободные колебания. Согласно второму закону Ньютона c учетом знаков проекций на ось X

(1)

Из условия равновесия смещение, вызываемое силой тяжести: . Подставляя в уравнение (1), получим: . Разделив правую и левую часть этого уравнения на m и принимая, что a = d2x/dt2, получим дифференциальное уравнение

. (2)

Это уравнение называется дифференциальным уравнением гармонических колебаний пружинного маятника. Из этого уравнения следует, что после прекращения внешнего воздействия, приводящего к первоначальному отклонению системы от положения равновесия, движение груза обусловлено только действием упругой силы (сила тяжести вызывает постоянное смещение).

Общее решение однородного дифференциального уравнения второго порядка (2) имеет вид

. (3)

Данное уравнение называется уравнением гармонических колебаний. Наибольшее отклонение груза от положения равновесия А0 называется амплитудой колебаний. Величина , стоящая в аргументе косинуса, называется фазой колебания. Постоянная φ0 представляет собой значение фазы в начальный момент времени (t = 0) и называется начальной фазой колебаний. Величина

(4)

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

. (5)

Рассмотрим свободные колебания пружинного маятника при наличии силы трения (затухающие колебания). В простейшем и вместе с тем наиболее часто встречающемся случае сила трения пропорциональна скорости υ движения:

где r – постоянная, называемая коэффициентом сопротивления. Знак минус показывает, что сила трения и скорость имеют противоположные направления. Уравнение второго закона Ньютона в проекции на ось Х при наличии упругой силы и силы трения

Данное дифференциальное уравнение с учетом υ = dx/dt можно записать

, (8)

где – коэффициент затухания; – циклическая частота свободных незатухающих колебаний данной колебательной системы, т. е. при отсутствии потерь энергии (β = 0). Уравнение (8) называют дифференциальным уравнением затухающих колебаний.

Чтобы получить зависимость смещения x от времени t, необходимо решить дифференциальное уравнение (8). В случае малых затуханий () решение уравнения можно записать следующим образом:

, (9)

где А0 и φ0 – начальная амплитуда и начальная фаза колебаний;
– циклическая частота затухающих колебаний при ω >> ω ≈ ω0.

Движение груза в этом случае можно рассматривать как гармоническое колебание с частотой ω и переменной амплитудой, меняющейся по закону:

. (10)

На графике функции (9), рис. 2, пунктирными линиями показано изменение амплитуды (10) затухающих колебаний.

Рис. 2. Зависимость смещения х груза от времени t при наличии силы трения

Для количественной характеристики степени затухания колебаний вводят величину, равную отношению амплитуд, отличающихся на период, и называемую декрементом затухания:

. (11)

Часто используют натуральный логарифм этой величины. Такой параметр называется логарифмическим декрементом затухания:

. (12)

Если за время t‘ амплитуда уменьшается в n раз, то из уравнения (10) следует, что

. (13)

Отсюда для логарифмического декремента получаем выражение

. (14)

Если за время t‘ амплитуда уменьшается в е раз (е = 2,71 – основание натурального логарифма), то система успеет совершить число колебаний

. (15)

Следовательно, логарифмический декремент затухания – величина, обратная числу колебаний, совершаемых за то время, за которое амплитуда уменьшается в е раз. Чем больше θ, тем быстрее происходит затухание колебаний.

2. Методика эксперимента и экспериментальная установка

Рис. 3. Схема установки

Установка состоит из штатива 1 с измерительной шкалой 2. К штативу на пружине 3 подвешиваются грузы 4 различной массы. При изучении затухающих колебаний в задании 2 для усиления затухания используется кольцо 5, которое помещается в прозрачный сосуд 6 с водой.

В задании 1 (выполняется без сосуда с водой и кольца) в первом приближении затуханием колебаний можно пренебречь и считать гармоническими. Как следует из формулы (5) для гармонических колебаний зависимость T 2 = f (m) – линейная, из которой можно определить коэффициент жесткости пружины k по формуле

, (16)

где – угловой коэффициент наклона прямой T 2 от m.

Задание 1. Определение зависимости периода собственных колебаний пружинного маятника от массы груза.

1. Определить период колебаний пружинного маятника при различных значениях массы груза m. Для этого с помощью секундомера для каждого значения m трижды измерить время t полных n колебаний (n ≥10) и по среднему значению времени вычислить период . Результаты занести в табл. 1.

2. По результатам измерений построить график зависимости квадрата периода T2 от массы m. Из углового коэффициента графика определить жесткость пружины k по формуле (16).

Результаты измерений для определения периода собственных колебаний

, с

, с

🎦 Видео

МЗЭ 2022 Численное решение дифференциальных уравнений. Неявный метод Эйлера. Ложкин С.А.Скачать

5 Численное решение дифференциальных уравнений Part 1Скачать

Численное решение задачи Коши методом ЭйлераСкачать

Дифференциальные уравнения для самых маленькихСкачать

Дифференциальные уравнения, 1 урок, Дифференциальные уравнения. Основные понятияСкачать

19. Метод вариации произвольных постоянных. Линейные неоднородные диф уравнения 2-го порядкаСкачать

Дифференциальные уравнения. Задача Коши. Метод Эйлера.Скачать

16. Линейные неоднородные дифференциальные уравнения 2-го порядка с постоянными коэффициентамиСкачать

Линейное неоднородное дифференциальное уравнение с постоянными коэффициентами 4y''-y=x^3-24x #1Скачать

13. Как решить дифференциальное уравнение первого порядка?Скачать

py127 Численное решение обыкновенных дифуровСкачать

Решение системы дифференциальных уравнений методом ЭйлераСкачать

Дифференциальные уравнения, 8 урок, Линейные дифференциальные уравнения с const коэф-ами 2 порядкаСкачать