Главная Случайная страница


Категории:

ДомЗдоровьеЗоологияИнформатикаИскусствоИскусствоКомпьютерыКулинарияМаркетингМатематикаМедицинаМенеджментОбразованиеПедагогикаПитомцыПрограммированиеПроизводствоПромышленностьПсихологияРазноеРелигияСоциологияСпортСтатистикаТранспортФизикаФилософияФинансыХимияХоббиЭкологияЭкономикаЭлектроника






Вывод уравнения колебания струны

Рассмотрим натянутую струну длины l, закрепленную на концах. В положении равновесия струна направлена вдоль оси Ox. Сила натяжения T0, действующая на струну, предполагается значительной. Каждую точку струны длины l можно охарактеризовать значением ее абсциссы x и смещением этой точки в момент времени t.

Для упрощения задачи примем следующие предположения:

1. Будем рассматривать только поперечные колебания струны, предполагая, что движение происходит в одной плоскости, и что все точки струны движутся перпендикулярно оси Ox.

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

2. Будем рассматривать струну как гибкую упругую нить:

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

· понятие «нить» означает, что мы пренебрегаем формой поперечного сечения и толщиной (рассматриваем линейную плотность ρ(x)).

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

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

,

 

где – удлинение струны; – начальная длина струны.

Таким образом,

,

где k – коэффициент упругости.

Длина произвольного участка струны (рис.2.1) в любой момент времени выражается формулой:

.

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

 

Покажем также, что натяжение не зависит и от x. Найдем проекции натяжения на оси x и u (обозначим их Tx и Tu):

;

,

где α– угол касательной к кривой с осью x.

На участок действуют силы натяжения и внешние силы. Сумма проекции всех сил на ось x должна быть равна нулю (мы рассматриваем только поперечные колебания, т.е. струна не движется вдоль оси Ох). Так как внешние силы по предположению направлены вдоль оси u, то

или .

Отсюда в силу произвольности x и следует, что натяжение не зависит от x, т.е. для всех значений x и t:

.

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

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

Для вертикальной составляющей сил натяжения имеем выражение

.

Так как рассматриваем малые колебания, то

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

(2.1)

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

(2.2)

где – линейная плотность струны. Приравнивая выражения (2.1) и (2.2) и переходя к пределу при , для искомой функции получим уравнение:

или

где .

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

или

.

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

2.2. Методы решения уравнения колебания струны

2.2.1. Метод Даламбера (метод бегущих волн) для бесконечной струны

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

Итак, рассмотрим задачу Коши для уравнения колебания струны:

 

(2.3)

где – функция, задающая форму струны в начальный момент времени; – скорость точки струны в начальный момент.

Уравнение решается в явном виде с помощью замены переменных :

где

;

Аналогично,

Подставляем в уравнение

Отсюда

.

Интегрируя это равенство последовательно по каждой переменной, получим:

.

Вернемся к старым переменным:

. (2.4)

Функция описывает волну, бегущую вправо со скоростью а, а функция – волну, бегущую влево.

Функция (2.4) является общим интегралом уравнения (2.3). Теперь необходимо удовлетворить начальным условиям:

Интегрируя последнее уравнение системы, получим:

где

Или

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

Отсюда

Подставляем эти выражения в формулу (2.4) и получаем решение волнового уравнения (формула Даламбера):

. (2.5)

Рассмотрим два частных случая.

Допустим, что , т.е. струне придана начальная форма при нулевой начальной скорости. Тогда решение (2.5) принимает вид

и, следовательно, представляет собой сумму двух бегущих волн: прямой волны и обратной волны . Первая перемещается по направлению оси ОХ, а вторая – в противоположном направлении.

При решение (2.5) имеет вид

.

Пусть

,

тогда

.

И в этом случае решение начальной задачи представляет собой сумму двух волн: прямой волны и обратной волны .

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

Пример. Решить уравнение колебания бесконечной струны , удовлетворяющее условиям:

Решение: Имеем задачу свободных колебаний бесконечной струны (без краевых условий). Применяем формулу Даламбера:

2.2.2. Фазовая плоскость

Для выявления характера решения волнового уравнения (2.5) удобно воспользоваться плоскостью состояний (x, t) или «фазовой плоскостью» (рис.2.2). Прямые и называются характеристикамиуравнения (2.5). Функция вдоль характеристики сохраняет постоянное значение, функция постоянна вдоль характеристики .

Рис.2.2. Характеристический треугольник MPQ фазовой плоскости

 

Рассмотрим некоторую фиксированную точку и проведем из нее обе характеристики и , которые пересекают ось ОX в точках и . Треугольник MPQ называется характеристическим треугольником точки . Отклонение точки струны в момент времени зависит только от значений начального отклонения в вершинах P и Q треугольника MPQ и от значений начальной скорости на стороне PQ:

.

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

Решение можно представить в виде суммы

где

.

 

Наглядное представление о характере процесса распространения можно получить с помощью фазовой плоскости (x, t). Проведем характеристики через точки (–l, 0) и (l, 0). Они разбивают плоскость на шесть областей (рис.2.3).

 

Рис.2.3. Фазовая плоскость для бесконечной волны

Рассмотрим два случая.

Пусть на отрезке .

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

Области 1, 6 – колебаний нет.

Область 2: – волна движется влево.

Область 5: – волна движется вправо.

Область 4: – волны складываются.

Область 3 – колебаний нет, отклонение равно нулю.

Пусть на отрезке .

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

Области 1, 6: колебаний и отклонений нет.

Область 2: – волна бежит влево с изменением формы.

Область 5: – волна бежит вправо с изменением формы.

Область 4: – волны складываются.

Область 3: – колебаний нет, но струна не возвращается в исходное положение (если постоянная не равна нулю).

Пример. Построить профиль струны при для различных моментов времени в случае

В этом случае прямая волна движется вправо, а обратная – влево (рис.2.4).

Рис. 2.4. Профили струны для различных моментов времени в случае нулевой начальной скорости

 

Пример. Построить профиль струны для различных моментов времени в случае

Найдем выражение для :

так как в силу непрерывности при : , а при : .

В этом случае прямая волна движется вправо, а обратная – влево (рис.2.5).

Рис.2.5. Профили струны для различных моментов времени в случае нулевого начального отклонения

 

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

a=1;

l=1;

dx=.01;

x=-4*l:dx:4*l;

u=1-abs(x);

u(abs(x)>l)=0;

u_left=.5*u;

u_right=.5*u;

for t=0:.25:1.25

subplot(3,2,4*t+1);

hold on

u1=circshift(u_left,[0 -a*t/dx]);

u2=circshift(u_right,[0 a*t/dx]);

plot(x,u1,'g-','lineWidth',2);

plot(x,u2,'b-','lineWidth',2);

plot(x,u1+u2,'r-','lineWidth',3);

xlim([-4*l 4*l]);

ylim([0 1]);

grid on

xlabel('x');

ylabel('u');

title(['t=' num2str(t)]);

end

2.2.3. Метод продолжения для полубесконечной струны

Рассмотрим задачу о распространении волн на полубесконечной прямой (x ≥ 0). Следует отметить, что чаще всего имеют дело с жестким или свободным закреплением струны.

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

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

Доказательство леммы 1.

Примем за начало координат, . В этом случае условия нечетности начальных данных запишутся в виде

Функция при равна

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

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

Доказательство леммы 2.

Условие четности начальных данных имеет вид:

Заметим, что производная четной функции является функцией нечетной: Рассмотрим производную:

;

.

Так как первое слагаемое равно нулю в силу нечетности , а второе – в силу четности .

Жесткое закрепление.

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

Задача ставится следующим образом: ищем решение системы уравнений

Рассмотрим функции и (x), являющиеся нечетными продолжениями функций и . Тогда функция

определена для всех . В силу леммы 1 .

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

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

Свободное закрепление.

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

Делаем четное продолжение функций и . Получим решение уравнения колебаний в виде функции

,

определенной для всех . В силу леммы 2 .

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

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

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

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

2.2.4. Метод продолжения для конечной струны (начальная и конечная точки жестко закреплены)

Рассмотрим краевую задачу для ограниченного отрезка (0, l). Будем искать решение уравнения

,

удовлетворяющее граничным условиям

и начальным условиям

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

,

где (x) и (x) – функции, подлежащие определению. Начальные условия

определяют значения (x) и (x) в интервале (0, l).

Для того чтобы удовлетворить нулевым граничным условиям, наложим на функции (x) и (x) требования нечетности относительно точек x = 0, x = l:

Сопоставляя эти равенства, получим:

и аналогично для Ψ(x), т.е. (x) и (x) – периодические функции с периодом 2l.

Нетрудно видеть, что условия нечетности относительно начала координат и условия периодичности определяют продолжение (x) и (x) на всей прямой . Подставляя их, получаем решение задачи.

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

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

Тогда по формуле Даламбера:

=

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

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

,

Тогда по формуле Даламбера:

=

Упражнения.

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

2.1.

Ответ: .

2.2.

Ответ: .

2.3.

Ответ: .

 

2.4.

Ответ: .

2.5.

Ответ: .

2.6.

Ответ: .

2.7.

Ответ: .

2.8.

Ответ: .

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

2.9.

Ответ:

2.10.

Ответ:

2.11.

Ответ: .

2.12.

Ответ:

2.13.


Ответ:

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

2.14.

Нарисовать профиль струны для моментов времени

Ответ:

2.15.

Нарисовать профиль струны для моментов времени

Ответ:

2.16.

Нарисовать профиль струны для моментов времени


Ответ:

2.17.

Нарисовать профиль струны для моментов времени


Ответ:

2.18.

Нарисовать профиль струны для моментов времени


Ответ:

2.19.

Нарисовать профиль струны для моментов времени

Ответ:

2.3. Метод Фурье (метод стоячих волн) или метод разделения переменных

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

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

где X (x)функция только переменного x; T (t)функция только переменного t.

Подставляя предполагаемую форму решения в уравнение, получим:

или, после деления на множитель ,

.

Правая часть этого равенства является функцией только переменного t, а левая – только х. Фиксируя, например, некоторое значение х и меняя t (или наоборот), получим, что правая и левая части равенства при изменении своих аргументов сохраняют постоянное значение:

,

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

Очевидно, что нас интересуют нетривиальные решения.

Граничные условия дают:

Отсюда следует

.

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

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

Рассмотрим отдельно случаи, когда параметр отрицателен, равен нулю или положителен.

1. При задача не имеет нетривиальных решений. Действительно, общее решение уравнения имеет вид

,

Граничные условия дают:

; .

Отсюда и, следовательно, .

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

.

Граничные условия дают:

; .

Отсюда и, следовательно, .

3. При общее решение уравнения может быть записано в виде

.

Граничные условия дают:

;

Нетривиальное решение получаем только в случае или . Отсюда

.

Этим собственным значениям соответствуют собственные функции

,

где – произвольная постоянная.

Пусть , тогда собственными функциями являются

.

Аналогично решаем уравнение относительно :

,

где и – произвольные постоянные.

Следовательно, функции

являются частными решениями данного уравнения. В силу линейности и однородности уравнения сумма частных решений также удовлетворяет этому уравнению и граничным условиям. Получаем общее решение:

.

Начальные условия позволяют определить и :

;

.

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

; .

Подставив эти коэффициенты в общее решение, мы удовлетворим краевым условиям и получим решение уравнения.

Простейшие задачи Штурма – Лиувилля для уравнения .

Вид условия Собственные значения и функции для
 
 
 
 

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

Решение. Общее решение имеет вид

.

Из начальных условий определим и :

Тогда .

.

Отсюда .

Подставив эти коэффициенты в общее решение, получим решение уравнения:

.

Можно построить в среде MATLAB поверхность решения данного волнового уравнения (рис.2.6).

Рис.2.6. Поверхность решения уравнения колебания ограниченной струны

[x,t]=meshgrid(0:.1:5);

u=1

Последнее изменение этой страницы: 2017-07-07

lectmania.ru. Все права принадлежат авторам данных материалов. В случае нарушения авторского права напишите нам сюда...